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BACKGROUND 


SUMMARY 


The  purpose  of  the  FY-71  effort  was  to  initiate  a program 
for  improving  vehicular  communications  in  the  urban  environment 
and  to  provide  UMTA  with  technical  support  in  this  and  related 
areas.  A general  program  outline  is  given  in  the  PPA  UM03-1 
dated  July  24,  1970. 

Requirements  imposed  on  land-mobile  communication  services 
have  become  more  stringent  in  recent  years.  First,  the  number 
of  users  continues  to  increase,  while  available  spectrum  for 
this  service  remains  severely  restricted  both  for  technical 
reasons  and  for  reasons  of  public  policy;  second,  the  emerging 
interest  in  data  communications  and  in  automatic  vehicle  mon- 
itoring ( AVM)  systems  implies  demands  for  increased  communi- 
cation capacity  utilizing  wideband  signalling  methods. 

Timely  provision  of  this  increased  capacity  requires  that 
we  understand  in  detail  the  urban  communication  environment. 

The  measurement  of  these  channel  properties,  therefore,  is  not 
academic;  in  fact,  it  is  integral  to  both  design  and  evaluation 
of  communication  systems  capable  of  carrying  out  the  difficult 
task  of  reliable  urban  communication  at  reasonable  cost. 

The  analytic  techniques  and  special  purpose  equipment  being 
developed  for  UMTA  by  TSC  have  immediate,  practical  application 
to  an  important,  continuing  UMTA  demonstration  program,  namely, 
the  AVM  tests  presently  planned  for  Philadelphia.  Two  main 
areas  of  application  are  envisaged:  (1)  spectrum  occupancy  and 

interference,  and  (2)  design  assistance  for  waveforms  and 
modems . 

1 . Spectrum  occupancy  and  interference 


The  AVM  program  can  be  assisted  in  several  task  areas  using 

the  TSC  equipment. 

a.  Quantitative  measurements  of  RFI  and  EMI  generated  by 
the  test  systems  can  be  made  and  their  effect  on  the 
EM  environment  in  their  bands  of  operation  can  be 
calculated.  These  measurements  can  help  set  filtering 
and  shielding  requirements  for  transmitters  and 
receivers . 

b.  These  data  will  provide  essential  information  to  be 
presented  to  the  FCC  when  requesting  approval  for  the 
use  of  certain  bands  or  modulation  formats. 

c.  The  TSC  equipment  will  be  available  to  carry  out  other 
spectrum  occupancy  measurements  in  response  to  FCC 
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requests  for  further  data. 


2 . Design  assistance 


The  channel  characterization  and  noise  measurement  equip- 
ment can  be  used  as  a design  tool  in  the  formulation  of 
range  tone  location  systems  and  in  the  design  and  test  of 
digital  modems. 

The  TSC  program  consists  of  two  parts  described  in  the  sections 
Experimental  Phase  and  Analytical  Phase. 

EXPERIMENTAL  PHASE 

First,  a mobile  experimental  station  is  being  developed 
with  the  capability  to  measure  and  record  both  the  multipath 
structure  of  any  particular  urban  channel  and  the  associated 
noise  environment.  This  will  be  accomplished  by  outfitting  a 
van  to  make  noise  measurements  and  also  to  be  the  mobile  receiv- 
ing site  for  suitable  designed  probing  signals  which  will  be 
transmitted  from  fixed  locations. 

A thorough  review  of  the  literature  has  supported  our  view 
that  the  present  measurements  and  standards  for  urban  noise  are 
inadequate.  Consequently,  we  are  preparing  to  make  broad-band 
noise  measurements  over  a wide  dynamic  range  in  order  to  preserve 
the  basic  noise  structure.  This  data  will  be  recorded  on  mag- 
netic tape  and  will  be  available  both  for  analysis  and  simula- 
tion. The  results  of  the  noise  and  multipath  measurements  will 
then  be  used  in  formulating  a clearer  description  of  these  pro- 
cesses. This  capability  will  also  be  available  to  make  compre- 
hensive measurements  on  the  electromagnetic  noise  produced  by 
electric  vehicles.  The  effects  of  suppression  techniques  can 
also  be  measured  and  compared. 

The  probing  signals  will  be  specifically  designed  to  decom- 
pose the  multipath  structure  into  its  fundamental  components. 

This  is  a standard  technique  in  the  study  of  fading  dispersive 
channels  which  is  being  carried  over  directly  to  the  urban  chan- 
nel. If  done  properly,  each  of  these  fundamental  components 
can  be  recorded  on  tape  and  used  as  an  input  to  a "playback" 
channel  simulator.  Therefore,  the  emphasis  of  the  first  part 
of  this  program  is  to  assemble  a channel  simulator,  designed  to 
test  hardware  for  the  particular  environment  in  which  a demon- 
stration is  to  take  place.  We  estimate  that  it  will  be  possible 
to  record  an  environment  and  be  prepared  to  evaluate  equipment 
in  a period  of  about  four  months.  Hopefully,  this  will  provide 
a quick  reaction  to  programs  and  changes  in  programs  as  they 
develop . 
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ANALYTICAL  PHASE 


The  second  part  of  this  program  is  directed  toward  analyz- 
ing the  noise  and  multipath  characteristics  and  developing  more 
reliable  configurations  for  command  and  control.  In  particular, 
since  UMTA  will  most  likely  be  interested  in  the  flow  of  digital 
data  between  terminals  it  will  be  necessary  to  study  specifica- 
tions and  requirements  for  digital  data  flow  over  the  existing 
voice  links.  At  present,  land  mobile  bands  are  allocated  for 
voice  transmission  and  only  recently  have  proposals  been  made 
for  digital  data  flow  in  the  same  channel.  Allocations  of  chan- 
nels for  digital  data  transmission  are  available  on  an  experi- 
mental basis.  The  intended  use  of  these  channels  is  exclusively 
for  Automatic  Vehicle  Monitoring  (AVM)  systems,  not  for  general 
digital  transmission  of  voice  and  other  message  traffic.  In  ad- 
dition, most  of  the  available  hardware  is  not  suited  for  this 
use  and  is  spectrally  inefficient. 

Thus,  effort  is  to  be  directed  toward  the  eventual  construc- 
tion of  a communication  system  simulator  into  which  the  channel 
simulator  could  be  inserted.  This  would  allow  detailed  design 
and  specification  of  command  and  control  systems  for  the  urban 
environment . 

This  approach  to  the  command  portion  of  the  problem  has  the 
following  merits. 

As  the  frequency  assignments  change,  it  is  relatively  easy 
to  retake  the  measurements  without  altering  the  basic  system. 
Since  the  assemblage  is  modular  in  nature  other  channels  of  in- 
terest can  be  inserted  without  losing  the  capability  of  the  com- 
munication simulator.  Furthermore,  the  system  can  be  directly 
employed  in  analyzing  and  designing  various  forms  of  AVM  systems 
which  use  the  land  mobile  bands.  The  great  flexibility  of  a 
communication  simulator  will  allow  the  designer  to  tailor  his 
system  to  a close  replica  of  the  channel  and  to  quickly  evaluate 
competing  techniques. 

In  addition,  we  will  be  in  a position  to  suggest  modifica- 
tions of  existing  equipment  that  will  yield  the  minimum  impact 
on  the  cost  of  demonstration  programs. 

RELATED  WORK 

During  the  course  of  the  program,  we  have  given  support  to 
the  demonstration  programs  office  in  the  assessment  of  competing 
AVM  systems.  This  has  been  directed  primarily  at  proposal 
evaluation  and  technical  recommendations.  The  level  of  future 
TSC  participation  in  the  AVM  program  will  be  determined  by  the 
demonstration  programs  office  and  will  be  responsive  to  their 
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needs . 


We  are  also  undertaking  preliminary  investigation  into  the 
general  area  of  traffic  flow  and  traffic  control  in  the  urban 
area.  This  work  is  being  directed  toward  assessing  the  status 
of  present  traffic  control  systems  and  determining  what  the 
state  of  the  technology  is  in  this  area.  Coordination  is  also 
being  made  with  other  groups  within  TSC  who  are  working  in  this 
and  related  areas. 
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TECHNICAL  DISCUSSION 


INTRODUCTION 

GENERAL 

A program  was  initiated  at  TSC  to  study  the  interference  on 
the  communication  channel  in  the  urban  and  suburban  environments 
and  to  establish  the  capability  for  comparative  evaluation  of 
competing  components,  subsystems  and  systems.  In  this  first 
year,  major  effort  was  concentrated  on: 

1.  Defining  the  characteristics  of  interference 

2.  Establishing  interference  measuring  facility 

3.  Constructing  a channel  simulator  facility. 

In  the  following  sections,  we  will  discuss  the  general  prob 
lem,  the  noise  program,  and  the  channel  characterization  program 

GENERAL  PROBLEM 

Reviewing  the  past  experimental  and  theoretical  work,  we 
believe  there  are  essentially  three  major  types  of  interference 
in  the  urban  and  suburban  environment;  namely,  additive  noise  -- 
Gaussian  and  impulsive,  Co-channel  interference,  and  fading  dis- 
persive interference  — reflection  and  scattering  off  the  struc- 
tures, terrain,  etc.  The  experimental  program  was  designed  for 
the  following  purposes: 

1.  To  study  experimental  techniques  for  the  measurement  of 
wideband  man-made  noise  in  the  VHF-UHF  region  and  to 
establish  a capability  for  the  routine  performance  of 
such  measurements. 

2.  To  identify  the  major  sources  of  radiation  that  will 
interfere  with  the  performance  of  any  land  transporta- 
tion communication  system  being  planned  in  the  frequen- 
cy band  mentioned  above. 

3.  To  determine  the  information  needed  for  the  design  and 
evaluation  of  any  communication  or  command  and  control 
system  of  mass  transit  or  land,  air,  water  or  highway 
communication.  In  combination  with  channel  probing 
measurements  this  information  will  provide  a complete 
characterization  of  the  urban  and  suburban  channel. 

4.  To  determine  the  radiation  characteristics  of  existing 
and  proposed  transportation  systems  such  as  electric 
cars,  trolley  cars,  electric  trains,  etc. 
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Figure  1 shows  the  major  tasks  and  problem  areas  and  the 
logical  flow  of  the  data  in  the  measurement  program. 

NOISE  MEASUREMENT  PROGRAM 


A communication  system  planner  would  like  to  know  the  de- 
tailed characteristics  of  the  background  noise  in  order  to  evalu 
ate  its  effects  on  system  performance.  This  information  include 
the  time  and  space  variation  of  the  noise  as  well  as,  for  predic 
tion  purposes,  its  correlation  with  factors  such  as  population 
density,  vehicular  traffic,  power  consumption  and  topographical 
characteristics . 

As  has  been  mentioned,  the  requirements  are  more  stringent 
for  digital  communication  systems,  which  require  a more  thorough 
specification  of  the  statistical  characteristics  of  the  noise 
than  do  analog  systems.  By  and  large  this  information  is  not 
currently  available  and  this  lack  of  design  data  hampers  design- 
ers of  wideband  systems,  i.e.,  for  our  purposes  defined  as  sys- 
tems whose  bandwidth  exceeds  10  kHz.  Therefore,  the  purpose  of 
the  experiment  under  consideration  is  to  measure  broadband,  wide 
dynamic  range,  impulsive-type  noise  in  selected  regions  of  the 
VHF  and  UHF  bands.  The  bandwidth  of  the  measuring  system  has 
been  tentatively  selected  to  be  1 MHz.  The  predominant  noise  is 
expected  to  be  generated  by  man-made  sources  such  as  automobile 
engines,  industrial  equipment,  etc. 

One  of  the  first  tasks  is,  of  course,  the  determination  of 
the  statistical  quantities  describing  the  noise  which  can  be  de- 
rived by  suitable  processing  of  the  experimental  data.  Although 
a certain  amount  of  data  processing  can  be  done  in  real  time  the 
decision  has  been  made  to  record  the  IF  waveform  for  later  proc- 
essing and  for  use  in  the  simulation  equipment.  This  means  that 
most  of  the  processing  procedures  can  be  decided  upon  at  a later 
date.  However,  some  initial  consideration  of  what  has  to  be 
measured  is  in  order.  A simple  model  of  wideband  man-made  inter 
ference  is  the  sum  of  an  impulsive  component  and  a Gaussian  proc 
ess.  Both  components  can  be  expected  to  be  nonstationary.  It 
seems  reasonable  to  measure  at  least  the  following  quantities. 


1.  Average  power  over  each  measurement  interval 

2.  Amplitude  probability  distribution  function  (pdf)  of 
the  peaks 

3.  Duration  pdf  of  the  impulses  at  given  amplitude  levels 

4.  Pulse  spacing  pdf 

5.  Crossing  rate  pdf 


In  addition  to  the  above  quantities  some  simpler  parameters 
such  as  average  and  rms  values  after  detection  can  be  measured 
in  real  time.  It  is  worth  pointing  out  that  instrumentation 
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Figure  1.  Diagram  of  the  Major  Tasks  and  Problem  Areas  and  the 
Logical  Flow  of  the  Data  in  the  Measurement  Program 


similar  to  that  being  studied  for  the  noise  measurements  can  be 
used  in  the  future  for  the  investigation  of  channel  occupancy. 

By  this  we  mean  the  evaluation  of  the  probability  that  a given 
frequency  region  be  occupied  at  a certain  level  by  other  trans- 
mitters. Because  of  the  limited  spectral  region  available  for 
land  mobile  communication  it  is  also  important  to  determine  the 
minimum  allowed  spatial  separation  of  stations  transmitting  in 
the  same  frequency  region  in  order  to  avoid  appreciable  cochannel 
interference.  This  type  of  measurement  can  also  be  performed  by 
the  system  under  consideration. 

A simple  block  diagram  of  the  noise  measuring  system  is 
shown  in  Figure  2. 

CHANNEL  CHARACTERIZATION  EXPERIMENT 

We  now  consider  the  experimental  work  planned  for  study  of 
the  second  component  of  interference  caused  by  the  fading  dis- 
persive characteristics  of  the  urban  channel.  This  component 
arises  from  vehicle  motion  through  a continually  changing  multi- 
path  environment.  It  is  of  equal  importance  to  measure  these 
deleterious  effects  and  to  analytically  model  them.  Measurement 
is  important  because  it  provides  a realistic,  repeatable  replica 
of  the  communication  channel  which  can  be  used  in  optimizing  or 
refining  a system  design,  or  in  evaluating  contending  candidate 
systems  at  different  times  on  an  equal  basis.  Analytic  work  in 
this  area  is  important  in  defining  the  form  of  a channel  simula- 
tor even  when  the  problem  is  too  complex  for  a detailed  mathemat- 
ical analysis,  and  in  providing  the  basic  design  for  suitable 
transmitter  waveforms  for  use  in  the  fading  environment. 

Based  on  extensive  research  carried  out  over  the  last  decade 
in  the  area  of  characterization  of  fading  dispersive  channels,  it 
is  now  practical  to  build  a channel  simulator  capable  of  present- 
ing to  a candidate  communication  system  the  simulated,  stored 
replica  of  the  distortion  characteristics  of  any  desired  channel 
displaying  multipath  delay,  Doppler  spread  and  additive  noise. 
This  channel  playback  approach  effectively  brings  the  field  con- 
ditions into  the  laboratory  and  makes  possible  the  comparative 
evaluation  of  competing  techniques  under  identical  conditions. 
This  stored  channel  simulation  can  represent  an  estimate  of  the 
channel  characteristics  based  upon  mathematical  modeling  or  can 
in  fact  be  a faithful  reproduction  of  the  channel  under  operating 
conditions.  The  availability  of  both  approaches  is  desirable  for 
a comprehensive  simulator  facility. 

The  proposed  technique  involves  two  procedures  with  corre- 
sponding equipment.  The  first  procedure  involves  field  tests 
employing  specially  designed  transmitter  probing  signals  which 
cover  the  frequency  bands  of  interest  in  the  planned  program 
areas.  The  probing  signals  are  processed  in  the  moving  vehicle 
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Figure  2.  Basic  Receiver  Block  Diagram 


to  extract  the  time-variant  system  function  of  the  propagation 
environment  or  medium  and  record  it  on  a multi-channel  tape  re- 
corder. In  addition,  the  additive  noise  would  be  recorded.  A 
tape  library  could  be  built  up  which  would  adequately  cover  all 
environments  for  which  systems  are  being  developed.  The  replay 
procedure  is  carried  out  using  these  tape  records.  The  time- 
varying  channel,  with  the  additive  noise,  would  be  played  back 
repeatedly  to  enable  tests  of  various  systems  to  be  made  under 
the  same  field  conditions.  Basic  block  diagrams  of  the  signal 
processing  operations  carried  out  in  the  channel  prober  and  play- 
back systems  are  shown  in  Figures  3 and  4. 

Figure  3 defines  the  basic  signal  processing  operations  of 
the  probing  system.  A periodic  probing  signal  with  complex  en- 
velope z (t)  is  filtered  at  the  transmitter  by  IF  and  RF  filters 
prior  to  transmission.  All  the  filtering  operations  in  the 
transmitter  are  lumped  together  as  one  filter,  called  the  trans- 
mitter filter,  which  has  impulse  response  h<p(t).  The  propagation 
medium  is  represented  by  the  complex  time-variant  impulse  re- 
sponse h(t,x)  and  the  additive  noise  by  the  complex  process  n(t). 
All  the  linear  operations  in  the  receiver  prior  to  the  correla- 
tion operations  are  lumped  together  into  one  receiver  filter  with 
impulse  response  h^(t).  The  output  of  the  filter  is  fed  to  a 
number  of  parallel  correlators,  only  one  of  which  is  shown  in  the 
figure.  Complex  notation  is  used  to  describe  the  correlation 
operation  as  a multiplication  of  the  receiver  filter  output  by 
the  complex  conjugate  of  a shifted  probing  signal  followed  by  a 
complex  low  pass  filtering  operation. 

Channel  playback  is  achieved  with  a tapped  delay  line  as 
shown  in  Figure  4.  The  complex  time-varying  multipliers  are  just 
equal  to  the  low  pass  filter  outputs  of  the  correlators  which 
yield  estimates  of  the  impulse  response  for  uniformly  spaced  val- 
ues of  delay  separated  by  x seconds.  Filtering  operations  neces- 
sary to  avoid  aliasing  in  the  digital  realization  of  the  playback 
channel  are  shown  in  dashed  lines. 

OUTLINE  OF  THE  FOLLOWING  SECTIONS 

The  specifications  for  all  the  major  equipments  are  given 
in  the  next  section.  The  planning  of  the  experimental  procedure, 
such  as  the  number  of  measurements,  frequency,  and  the  length  of 
each  measurement  will  be  discussed.  Finally,  a milestone  for 
FY72  is  given. 

EQUIPMENT  SPECIFICATIONS 

RECEIVER 

The  function  of  the  receiver  is  as  follows:  To  make  both 

wideband  or  narrowband  measurements  of  man-made  noise,  spectrum 
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Figure  3.  Definition  of  Basic  Signal  Processing  Operations  in  Prober  Systems 
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Figure  4.  Channel  Playback  Schematic  Diagram 


occupancy  and  channel  characterization  of  the  time  varying  propa 
gation  media.  These  functions  are  carried  out  separately  in  time 
The  main  receiver  characteristics  are  the  following: 


1.  Three  selectable  rf  frequencies  and  their  respective 
tunable  frequency  range: 


Frequency  (MHz) 


Range  (MHz) 


149.95 

416.6 

902.2 


135-165 

400-480 

850-950 


2.  Selectable  i-f  predetection  BPF  (3  dB) : 

2,1  and  0.2  MHz 

3.  Available  outputs: 

a.  I-F  waveform 

b.  Video  amplitude  for  direct  magnetic  tape  recording 

c.  Two  orthogonal  components  of  the  phase  information, 
for  direct  magnetic  tape  recording 

d.  RMS  monitoring  and  recording  for  either  magnetic 
tape  or  strip  chart  recording  . 

The  receiver  should  have  a wide  dynamic  range  and  sensitiv- 
ity to  detect  any  conceivable  man-made  noise  (intentional  and 
unintentional)  in  the  urban  and  suburban  environment. 


Detailed  Specifications: 

Input  Impedance: 

50  ohms,  VWSR  < 1.5 

RF  Frequency: 

1.  149.95  MHz 

Preselector  rf  passband  (3  dB) : 4.3  MHz 

Tunable  range:  135-165  MHz 

Image  rejection:  80  dB  (min) 

Input  power  range:  -104  to  -19  dBm  over 

1 MHz  bandwidth 

2.  416.6  MHz 

Preselector  rf  passband  (3  dB) : 5 MHz 

Tunable  range:  400-480  MHz 

Image  rejection:  80  dB  (min) 

Input  power  range:  -104  to  -20  dBm  over 

1 MHz  bandwidth 
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3. 


902.2  MHz 

Preselector  rf  passband  (3  dB) : 18  MHz 

Tunable  range:  850-950  MHz 

Image  rejection:  80  dB  (min) 

Input  power  range:  -108  to  -23  dBm  over 

1 MHz  bandwidth 


Noise  Figure  (max) 


Frequency 


9 dB 
9 dB 
6 dB 


149.95  MHz 
416.6  MHz 
902.2  MHz 


Local  Oscillators:  In  addition  to  three  oscillators 

for  the  three  center  frequencies , an  external  ref- 
erence oscillator  port  should  be  provided. 


Stability:  ±5  kHz 


I-F  Frequency:  30  MHz 


I-F  Passband:  Selectable  and  switchable  according  to 

the  following  predetection  bandwidth: 

1.  2 MHz  ± 10% 

2.  1 MHz  ± 10% 

3.  0.2  MHz  ± 10% 


Rise  Time: 


1 

2X  bandwidth 


Overshoot:  <10  percent 


Outputs:  All  outputs,  with  the  exception  of  i-f  output, 

should  have  a maximum  compressed  dynamic  range  of 
30  dB  in  order  to  match  the  dynamic  range  capabil- 
ity of  an  analog  tape  recorder. 


1.  I-F  Output: 

Accuracy:  ±1  dB 

Maximum  level:  10  dBm/MHz 

Minimum  level:  -80  dBm/MHz 

I-F  Buffer  amplifier  should  be  provided  with  50  ohms 
output  impedance. 

2.  Baseband  Output: 

Accuracy:  ±1  dB 

Maximum  level:  5 V 

Minimum  level:  0.1  V 

Compressed  dynamic  range:  30  dB 

Output  impedance:  75  ohms 
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3.  Phase  Information  Output: 

Two  orthogonal  components  of  the  phase  information,  or 
alternatively,  the  in-phase  and  quadrature  components 
of  the  compressed  i-f  waveform. 

Phase  accuracy:  ±10° 

Level:  1 Vrms  nominal 

Phase  channel  center  frequency:  0.75  MHz 

Bandwidth  (3  dB) : 1 MHz 

The  output  bandwidth  should  be  compatible  with  the 
bandwidth  of  the  tape  recorder  (400  Hz  to  1.5  MHz). 

4.  RMS  output: 

5.  Time  constants:  Variable  in  steps  of  0.01,  0.1,  1 and 

10  seconds. 

Figure  5 shows  the  block  diagram  of  the  receiver. 

TRANSMITTER 

The  transmitter  has  the  function  of  transmitting  coded  sig- 
nals (binary  PM)  for  measuring  the  multipath,  frequency  broaden- 
ing, pulse  broadening  and  other  parameters  while  the  receiving 
station  is  moving  through  urban  and  suburban  environments.  It 
should  be  capable  of  transmitting  three  selected  frequencies, 
separately  in  time.  The  main  transmitter  characteristics  are  the 
following : 

1.  Three  selectable  rf  frequencies: 

149.95  MHz 
416.6  MHz 
902.2  MHz 

2.  Output  power:  20  W (min) 

Detailed  Specifications: 

Frequency:  149.95  MHz 

416.6  MHz 
902.2  MHz 

Inputs : 

1.  Source  Input 

Frequency:  79.95,  346.6  and  381.1  MHz 

Power  level:  0 dBm 

Input  impedance:  50  ohms,  VSWR  <1.5 
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Figure  5.  Block  Diagram  of  the  Receiver 


2.  L.O.  Input 

Frequency:  70  MHz 

Power  level:  -7  dBm 

Input  impedance:  50  ohms 

VSWR  <1.5 


Wideband  Mixers 

Operating  frequency  range:  140-460  MHz 

Output  level:  16  dBm  (min) 

Outputs : 

1.  149.95  MHz 

Bandwidth  (3  dB) : 5 MHz 

Power:  20  W (CW) 

2.  416.6  MHz 

Bandwidth  (3  dB) : 5 MHz 

Power:  20  W (CW) 

3.  902.2  MHz 

Bandwidth  (3  dB) : 10  MHz 

Power:  20  W (CW) 

Linearity:  +1.5  dB  for  all  amplifiers 

The  block  diagram  is  shown  in  Figure  6. 

MAGNETIC  TAPE  RECORDER  1 

This  magnetic  tape  recorder  has  the  purpose  of  recording 
all  the  information  on  man-made  noise.  The  information  to  be 
recorded  is  the  phase,  the  compressed  baseband  amplitude  and  the 
rms  level  for  further  processing.  It  should  have  the  tape 
speeds:  120  ips  (max)  to  1-7/8  ips  (min).  The  main  specifica- 

tions for  the  tape  recorder  are  as  follows: 

Number  of  channels:  7 

Flutter  (peak-to-peak) : 0.15%  at  120  ips 

Dynamic  interchannel  time-displacement  error:  0.2  ys  at 

120  ips 

Capstan  speed  accuracy:  ±0.01  percent 

Instantaneous  time  base  error:  less  than  +500  ns  at  120  ips 
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Figure  6.  Block  Diagram  of  the  Receiver 


Direct  record/reproduce  channel 

Input  sensitivity:  0.2  to  10  Vrms 

Input  impedance:  75  ohms 

Output  impedance:  75  ohms 

S/N  ratio:  30  dB  at  120  ips 

Bandwidth:  400  Hz  to  1.6  MHz  at  120  ips 

Frequency  modulation  record/reproduce  channel 
Input  sensitivity:  0.5  to  12.5  V peak 

Input  impedance:  20,000  ohms 

Output  impedance:  75  ohms 

S/N  ratio:  50  dB  at  120  ips 

Bandwidth:  dc  to  80  kHz  at  120  ips 

MAGNETIC  TAPE  RECORDER  II 

The  function  of  this  magnetic  tape  recorder  is  to  record  the 
outputs  of  a 10  tap-delay  line.  Both  in-phase  and  quadrature 
components  of  the  tap  outputs  will  be  recorded.  The  tape  record- 
er should  have  the  tape  speeds:  60  ips  (max)  to  1-7/8  ips  (min). 

The  main  specifications  for  the  tape  recorder  are  the  fol- 
lowing : 


No.  of  channels:  14 

Capstan  speed  accuracy:  ±0.05%  maximum 

Flutter  (peak-to-peak) : 0.7%  at  lh  ips 

Bandwidth:  DC  to  1250  Hz  at  ips 

S/N  ratio:  46  dB  at  lh  ips 

Input  level:  1 Vrms 

Input  impedance:  75  ohms 

Output  impedance:  1000  ohms 

Dynamic  range:  35  dB 

DATA  PROCESSOR 

The  function  of  the  data  processor  is  to  analyze  test  data 
resulting  from  the  man-made  noise  measurements,  both  real  time 
while  the  van  is  in  mobile  operation,  and  in  more  detail  follow- 
ing the  taking  of  test  data.  The  data  processor  consists  of  the 
following  major  components: 

1.  Interface  between  the  magnetic  tape  recorder  or 
receiver  video  output  and  a mini-digital  computer. 

2.  A mini-digital  computer  which  has  the  following  main 
subsystems : 

a.  Central  processor 

b.  Four  4K  16-bit  read/write  memories 
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c.  High  speed  paper  tape  reader  (300  cps)  and  punch 
(50  cps) 

d.  Analog  to  digital  convertor 

3.  Graphic  display  system 

Character  plotting  rate:  5000  characters/sec 

Point  plotting  rates:  100K  pt/sec 

Appendix  A gives  a detailed  description  of  the  data  process 
ing  procedure. 

In  addition  to  the  above  major  measuring  instruments,  a mo- 
bile van  was  furnished  and  modified.  The  van  will  be  used  as  a 
mobile  laboratory  which  will  house  the  measuring  equipment  such 
as  receiver,  tape  recorders,  data  processor,  etc.  The  measure- 
ments will  be  made  with  two  operators  in  the  van.  The  locations 
and  states  (stationary  and  moving)  of  the  van  will  be  discussed 
in  the  following  section. 

DESIGNS  OF  EXPERIMENTS 

NOISE  MEASUREMENT 

The  characterization  of  the  man-made  noise  present  in  the 
urban  environment,  and  therefore  the  design  of  an  appropriate 
experiment  to  derive  such  characterization,  appears  to  be,  at 
first,  a rather  hopeless  task  because  of  the  multivarious  vari- 
ability of  that  noise  with  a number  of  variables  such  as  time, 
space,  hour  of  the  day,  traffic  conditions,  etc.  In  addition, 
the  difficulty  of  the  task  is  compounded  by  the  fact  that  some 
of  those  variations,  for  example  spatial  variations,  are  seldom 
smooth  owing  to  occasional  sharp  discontinuities  in  the  distribu 
tion  of  noise  sources  and  to  the  presence  of  standing  wave  noise 
patterns.  Fortunately,  however,  a very  large  number  of  similar 
experiments  have  already  been  performed  by  other  investigators 
and  we  can  benefit  from  their  findings  and  try  to  improve  our 
knowledge  of  the  urban  noise  in  the  iterative  fashion  that  is 
typical  of  statistical  experimentation. 

The  first  important  subdivision  of  the  characterization 
problem  comes  from  recognizing  that  there  are  basically  three 
types  of  terminals  in  mobile  communication: 

1.  stationary  vehicles 

2.  moving  vehicles 

3.  base  stations. 

We  will  thus  break  the  design  of  our  experiment  into  three  cate- 
gories corresponding  to  the  three  types  of  terminals. 


20 


For  stationary  vehicles  we  can  start  with  a concept  put 
forward  in  the  FCC  study,  that  of  "activity  zone".  The  FCC  Re- 
port divides  the  total  urban  and  suburban  environment  into  areas 
characterized  essentially  by  their  level  of  man-made  noise.  Gen- 
eral categories  defined  in  the  study  were:  clear  areas,  resid- 

ences, businesses,  industrial  sites,  transportation  areas,  in- 
stitutional sites,  automotive  activities.  For  the  purpose  of 
locating  most  of  these  zones  in  a large  city  like  the  Greater 
Boston  area,  use  can  be  made  of  an  Urban  Atlas  recently  published 
This  atlas  contains  detailed  information  on  the  characteristics 
of  20  large  U.S.  cities  including  Boston.  Population  density  is, 
for  example,  given  with  a resolution  grid  of  250  m X 250  m and 
the  different  grids  of  each  city  are  divided  in  five  categories 
(50-200,  200-500,  500-1200,  1200-3600,  3600  and  more  people/grid) 
In  the  same  document  the  following  subdivisions  are  also  identi- 
fied for  a large  city:  Airport,  Institutional,  Large  Institu- 

tional, Park,  Industrial,  Commercial. 

Examination  of  the  FCC  measurements  and  some  additional 
thought  suggests  the  following  locations  where  measurements 
should  be  made  with  a stationary  vehicle. 

1.  Clear  areas 

2.  Residential  - separation  between  homes  500  ft. 

3.  Residential  - separation  less  than  30  ft.  but  detached 

4.  Residential  - row  houses 

5.  Apartment  buildings  - 2 stories 

6.  Apartment  buildings  - 6 stories  and  higher 

7.  Office  buildings  - 8 stories  and  more 

8.  Industrial  plant 

9.  High  power  lines  (220  kV) 

10.  Low  power  lines  (33  kV) 

11.  Power  plant 

12.  Hospitals 

13.  Main  thoroughfare 

14.  Busy  intersection 

15.  Large  institution 

16.  Airport 

It  is  expected  that  these  "activity  zones"  will  be  sufficiently 
different  in  the  amount  and  type  of  man-made  noise  present,  to 
justify  this  subdivision.  The  experiment  will,  of  course,  con- 
firm or  modify  this  assumption.  However,  even  if  the  noise 
characteristics  were  found  to  be  very  similar  for  two  or  more 
categories,  the  equivalence  of  different  types  of  sites  (say  a 
power  line  and  a congested  traffic  situation)  will  certainly 
constitute  important  information.  It  is  also  expected  that  for 
each  site  there  will  be  at  least  one  "control  variable"  to  take 
into  consideration:  for  example,  the  hour  of  the  day  in  a resi- 

dential area,  the  traffic  at  a busy  intersection  (smooth  traffic 
versus  a traffic  jam),  weather  conditions  for  a power  line,  etc. 


21 


Economy  of  experimentation  suggests  the  consideration  of  only 
two  levels  for  this  variable.  It  seems  also  advisable  to  repli- 
cate the  same  experiment  at  least  twice  for  each  type  of  site  in 
order  to  gain  some  insight  into  the  validity  of  the  assumption 
of  similarity.  It  must  also  be  mentioned  that,  owing  to  the  fact 
that  automotive  noise,  whenever  present,  is  an  important  source 
of  noise,  it  will  be  probably  impossible,  or  at  least  very  dif- 
ficult, to  measure  the  noise  characteristics  of  some  of  the  sites 
mentioned  without  including  automotive  noise.  This,  however, 
must  be  considered  a fact  of  life  (office  buildings  are  generally 
located  in  areas  of  high  traffic)  and  it  is  therefore  appropriate 
and  relevant  to  measure  the  total  noise  present  at  those  sites. 

The  total  number  of  measurements  required  will  in  general  be 
N = C X K X M 


where 


C = number  of  activity  areas  identified 

K = number  of  levels  of  the  control  variable 

M = number  of  replicates  for  the  same  activity  area 

From  the  above  considerations  it  seems  that  a minimum  of  64  mea- 
surements is  necessary  for  stationary  vehicles.  The  number  of 
sites  is  practically  coincident  with  the  number  of  measurements 
since,  in  general,  changing  the  level  of  the  control  variable 
will  prevent  using  a site  for  more  than  one  measurement. 

The  time  to  be  spent  on  each  measurement  depends  on  several 
factors.  In  general,  it  will  be  the  sum  of  two  terms:  time  nec- 

essary for  preparatory  and  concluding  procedures  (calibration, 
validation  of  measurements,  etc.)  and  the  measurement  time  itself. 
The  former  will  be  dictated  by  practical  considerations,  the  lat- 
ter is  related  to  the  assumed  stationarity  of  the  environment 
and  to  the  accuracy  with  which  we  estimate  probability  density 
functions  and  cumulative  density  functions.  Here  we  must  dis- 
tinguish between  on-line  processing  and  off-line  processing.  It 
seems  that  as  far  as  on-line  processing  is  concerned,  what  one 
can  determine  is  the  cumulative  distribution  on  a point-by-point 
basis,  assuming  stationarity  of  the  underlying  distribution.  One 
gets,  in  essence,  an  estimate  of  p^,  i.e.,  the  probability  that 
a certain  variable,  say  amplitude,  exceeds  the  level  £.  The 
estimate  is  clearly  given  by  x^/n  where  x^  is  the  number  of 
samples  which  exceed  the  level  and  n is  the  total  number  of 
samples.  The  random  variable  p^,  for  each  level,  can  be  thought 
as  the  population  parameter  in  a binomial  distribution.  The  bi- 
nomial distribution  can  be  approximated  by  a Gaussian  distribu- 
tion when  n becomes  large.  Since  the  mean  of  the  population 
x^/n  is  p^  and  its  variance  is  ppqf/n  where  q^  = 1-p^  the 
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standard  form  of  the  variable  x^/n  is 

z = (Vn)  ~ pt 
1 V p*Vn 


A confidence  interval  for  the  estimate  of  can  easily  be  con- 
structed. The  1-a  confidence  interval  for  2p  is 


< 


<Vn>  - p>. 

Vp*vn 


Z /o 
a/2 


The  1-a  confidence  interval  for  p can  be  obtained  by  solving 
the  previous  inequality  for  p^.  Since , however,  we  do  not  know 
the  standard  deviation  which  is  a function  of  p^  we  must  also 
estimate  it  from  the  data.  The  most  logical  estimate  is 


Est 


J¥-  / 


(x£/n) [1  - (x£/n) ] 


n 


and  therefore  the  confidence  interval  for  p is 
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The  number  n,  from  the  sampling  theorem,  is  given  by  the  expres- 
sion 


n = 2WT 


where  W is  the  bandwidth  of  the  receiver  and  T the  duration  of 
the  measurement  for  that  particular  level.  Thus,  if  one  assumes 
stationarity , a clear  tradeoff  exists  between  the  accuracy  re- 
quired (note  that  the  accuracy  is  a function  of  p^)  and  the  num- 
ber of  levels  that  one  is  willing  to  consider.  If  one  does  not 
want  to  assume  stationarity,  but  an  estimate  is  required  of  the 
validity  of  this  assumption,  the  total  time  must  be  also  sub- 
divided in  intervals  corresponding  to  possible  time  constants 
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of  stationarity.  To  get.  a feeling  for  the  order  of  magnitude 
involved  assume  a total  measurement  time  (per  frequency)  of  100 
sec.  Assume  that  we  want  to  test  stationarity  in  10  intervals 
of  10  sec.  each  and  that  we  quantify  the  pdf  in  10  levels.  The 
time  available  for  each  level  is  therefore  1 sec.  Assume  a band 
width  of  1 MHz.  Then  n = 2 X 106.  Assume  also  that  we  are  test 
ing  a level  of  low  probability , say  x/n  = 10” 5 . The  rms  error 
in  estimating  p^  is 


n 


10-6 


and  the  relative  error  becomes 


3.16  X 


10~6 

10“5 


3.16  X 10"1  . 


Note  that  if  the  hypothesis  of  stationarity  is  sustained,  one 
can  obviously  "pool"  the  estimates  of  variance  and  reduce  it 
further . 

Similar  results  hold  for  off-line  processing  where  one  de- 
termines a pdf  simultaneously  for  a chosen  number,  say  k,  of 
quantization  intervals.  One  has  then,  in  essence,  a multinomial 
distribution  with  k classes,  where  f]_,  f2,  . ...  f^  are  the  ex- 
perimental frequencies  and  pj_,  P2,  . ..,  Pk  the  true  values  of 

the  probabilities.  A simple  application  of  a theorem  due  to 
Fischer  shows  that  for  large  n (i.e.,  n >>  1,  n >>  fj_) 

Var  (fj_  n)  ~ 1 

= — for  each  i . 

(fi  n)2  f f 

The  characterization  of  the  urban  noise  for  a terminal  on  a 
moving  vehicle  should,  in  principle,  require  the  duplication  of 
the  measurements  taken  at  a stationary  site.  It  seems,  however, 
that  the  main  difference  between  a moving  vehicle  and  a station- 
ary one  arises  when  the  moving  vehicle  is  immersed  in  traffic. 
This  point  of  view  was  implicitly  taken  in  the  FCC  measurements 
where  the  only  measurements  from  moving  vehicles  were  made  while 
riding  in  downtown  traffic  and  on  a busy  suburban  thoroughfare. 
These  measurements  should  be  duplicated,  perhaps  adding  a few 
additional  ones  in  other  "activity  zones"  to  compare  the  results 
with  those  obtained  from  a stationary  vehicle. 
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Measurements  of  urban  noise  from  base  stations  differ  from 
the  ones  from  stationary  vehicles,  essentially  because  of  the 
height  of  the  receiving  antenna.  Furthermore,  most  base  stations 
are  presently  located  either  downtown  or  in  busy  suburban  dis- 
tricts. It  is  suggested  that  the  top  of  the  TSC  building  be 
used  for  a set  of  such  measurements.  Two  or  three  additional 
measurements  with  antennas  located  at  various  heights  will  prob- 
ably suffice  for  the  characterization  of  the  base  station  noise. 
Remarks  similar  to  the  ones  made  for  the  accuracy  achievable 
apply  here  also. 

SPECTRUM  OCCUPANCY 

The  functions  of  the  spectrum  occupancy  measurement  are  the 
following : 

1.  To  determine  the  average  rate  of  occupancy  in  the  land 
mobile  broadcasting  bands  150-160  MHz  and  450-470  MHz. 

2.  To  develop  techniques  for  determining  the  interference 
from  a large  distance,  say  a city  50  miles  away. 

3.  To  determine  the  center  frequency  and  i-f  bandwidth  to 
be  used  for  noise  measurement,  hence  avoiding  the  con- 
fusion of  two  distinctive  types  of  interference,  i.e., 
narrowband  - cw  and  wideband  impulsive. 

Site  Selection 


For  the  implementation  of  the  first  two  measurement  func- 
tions, the  initial  choice  of  the  roof  of  TSC  management  building 
seems  to  be  an  appropriate  one.  It  is  one  of  the  highest  build- 
ings in  a radius  of  two  miles. 

The  area  covered  by  transmission  from  the  management  build- 
ing may  include  the  Greater  Boston  Area,  i.e.  within  Route  128. 
This  seems  to  be  a good  start  initially  and  any  further  expan- 
sion of  the  terminals  will  be  decided  by  the  results  of  initial 
measurements  taken  from  the  roof  of  TSC  management  building. 

For  the  implementation  of  the  last  measurement  function  the 
sites  selected  are  similar  to  the  sites  selected  for  noise  mea- 
surements. The  measurements  will  be  taken  from  the  mobile  van. 
Because  of  the  proximity  of  the  antenna  to  the  ground  and  to  the 
urban  environment,  streets,  buildings,  expressway  etc.,  it  is 
not  expected  that  the  characteristics  of  the  spectrum  occupancy 
at  the  two  frequency  bands,  150-165  and  450-470  MHz,  will  be  the 
same  at  these  sites.  Accordingly,  the  total  number  of  the  sites 
is  the  same  as  the  number  of  activity  areas  defined  in  the  noise 
measurements . 
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Monitoring 


In  establishing  the  measurement  procedure  a HP  spectrum  an- 
alyzer (HP-855B/851B)  and  a recording  system  described  in  Appen- 
dix A will  be  used  to  monitor  the  spectrum  at  the  site. 

Three  frequency  bands  will  be  monitored:  135-165,  400-480 

and  850-950  MHz.  The  total  spectrum  is  to  be  monitored  at  the 
site  in  21  bands  of  10  MHz  each.  Each  band  will  be  recorded 
only  for  a 5-minute  run,  with  a complete  cycle  through  all 
twenty-one  bands  requiring  slightly  less  than  two  hours.  This 
pattern  will  be  repeated  four  times  each  day.  In  addition  to 
daytime  monitoring,  nighttime  monitoring  may  be  needed  if  the 
results  of  initial  daytime  monitoring  so  indicate.  The  duration 
of  nighttime  monitoring  will  be  decided  upon  after  analysis  of 
the  initial  data. 

CHANNEL  MEASUREMENT 

This  section  describes  the  planning  of  the  urban  channel 
measurements  in  the  Greater  Boston  Area.  The  site  selection  of 
the  transmitting  terminals,  the  characteristic  regions  of  urban 
and  suburban  environments,  and  the  average  time  for  individual 
measurements  are  discussed. 

Transmitting  Terminals 


Presently  three  different  types  of  terminals  are  planned. 
These  three  terminals  are:  a)  the  top  of  a highrise  building  in 

urban  area;  b)  the  top  of  a low  rise  building  in  urban  area;  and 
c)  a high  tower  in  a suburban  area.  For  the  above  three  ter- 
minals we  have  selected  the  following: 

1.  The  roof  of  the  TSC  management  building  because  it  is 
one  of  the  highest  buildings  in  the  area. 

2.  The  roof  top  of  TSC  Technology  building,  because  it  is 
a typical  low  rise  building  in  the  area. 

3.  The  transmitting  tower  of  New  England  Telephone  Co.  in 
Prospect  Hill,  Waltham,  a suburb  of  Boston. 

These  terminals  with  the  exception  of  2,  we  believe,  would 
represent  typical  terminals  for  all  land  mobile  base  stations. 
With  the  use  of  a low  terminal  (type  2),  we  hope  to  determine 
the  effect  of  transmit  antenna  height  in  urban  communications. 

Zoning  Areas 

In  the  urban  and  suburban  environments  one  can  characterize 
zones  according  to  structures  of  the  following  types: 
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1.  High  rise  buildings:  In  this  area  all  the  high  rise 

buildings  are  congested  together.  An  example  of  that 
is  the  downtown  Boston  area  near  Government  Center. 

2.  High  rise  buildings  scattered  among  low  rise  buildings: 
In  this  area  a few  high  rise  buildings  are  situated 
among  the  low  rise  buildings.  Typical  examples  are  the 
Back  Bay  area  of  Boston  and  the  area  near  TSC  in  Cam- 
bridge . 

3.  Low  rise  buildings:  In  this  area  there  is  no  high  rise 

building  over  ten-stories.  Brookline,  the  western  part 
of  Cambridge,  and  Somerville  are  a few  typical  areas  in 
the  Greater  Boston  area. 

4.  Individual  houses  located  on  a plain:  This  is  a typi- 

cal suburban  area  where  all  the  houses  are  located  on 
individual  lots.  Most  of  the  suburbs  of  Boston  fit  in 
this  category. 

5.  Houses  located  among  the  hills:  This  is  another  type 

of  the  suburbs.  There  are  a number  of  places  within 
the  Greater  Boston  area  which  fit  this  description,  for 
example,  Belmont  Hill. 

Measurement 


Initially,  five  to  ten  runs  in  each  zone  for  each  terminal 
are  planned  at  all  three  frequencies,  149.95,  416.6  and  902.2 
MHz.  Each  run  will  take  about  10  to  15  minutes.  The  basis  for 
measuring  time  is  the  requirement  that  each  measurement  should 
cover  a distance  of  three  to  five  miles,  on  the  assumption  that 
the  van  is  moving  between  10  and  30  miles  per  hour. 

MILESTONES  FOR  FY’72 

During  FY'71  we  were  involved  in  getting  to  know  the  prob- 
lems specific  to  urban  communication  systems,  command  and  control 
of  urban  vehicles,  and  voice  and  digital  communication  in  these 
channels.  We  found  that  there  are  three  major  interference 
sources  which  degrade  a communication  channel: 

a.  Man-made  noise 

b.  Co-channel  interference 

c.  Multipath  due  to  reflections  off  buildings  and 
terrain. 

With  a view  to  characterizing  the  above  three  interference 
sources,  a measurement  program  was  initiated  in  FY'71.  During 
FY'71  we  have  reviewed  the  past  work  on  measurement  and  charac- 
terization of  the  electromagnetic  interference,  have  planned  and 
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designed  a versatile  measurement  scheme  for  both  man-made  noise 
and  multipath  interference,  and  have  initiated  procurement  of  all 
the  necessary  equipment  to  establish  a mobile  laboratory.  This 
laboratory  includes  a data  processor  and  channel  analyzer  to  be 
used  for  characterizing  the  inter ference/  and  for  evaluating  pres- 
ent and  future  communication  systems. 

For  FY'72  the  emphasis  will  be  placed  on  actual  measure-" 
ments,  analysis  and  simulation.  The  following  table  (Table  1) 
shows  the  time  schedule  of  various  tasks. 

INSTALLATION  AND  CHECKOUT  OF  THE  MOBILE  VAN 

The  mobile  van  was  delivered  at  the  end  of  FY!71.  During 
the  first  two  and  a half  months  of  FY'72,  the  following  tasks 
are  planned: 

a.  Installation  of  antennas 

b.  Layout  of  connecting  cables  between  the  equipment 

c.  Arrangement  equipment  and  console  such  that  one 
operator  can  conduct  the  noise  measurement. 

Checkout  of  the  mobile  van  includes  the  following: 

a.  Performance  of  the  van 

b.  Performance  of  the  cooling  and  heating  systems 

c.  Performance  of  the  5 kW  power  generator 

d.  EM  noise  generated  by  the  van  and  its  accessories 
and  find  the  proper  way  to  suppress  them  if  neces- 
sary . 

CHANNEL  OCCUPANCY  MEASUREMENTS 

While  the  van  is  being  checked  out,  the  channel  measure- 
ments can  be  conducted  concurrently.  The  test  instrument  re- 
quired is  a spectrum  analyzer  which  is  presently  available  at 
the  center.  The  channel  occupancy  measurements  serve  three  func- 
tions as  mentioned  in  the  section  Spectrum  Occupancy.  These 
measurements  should  be  made  before  conducting  the  measurements 
on  man-made  noise.  We  plan  to  make  these  measurements  from  mid- 
July  to  the  end  of  September.  The  total  period  is  about  two  and 
a half  months. 

CALIBRATION  AND  TEST  INSTRUMENTS 

The  major  experiment  components  are  expected  during  the 
first  few  months  of  FY'72.  Test  and  calibration  of  these  compo- 
nents are  required  before  field  measurements  are  taken.  Further- 
more, TSC  personnel  may  need  to  be  trained  to  operate  the  instru- 
ments. A proper  procedure  of  data  taking  will  be  planned. 
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Table  1.  Time  Schedule  for  Various  Tasks 
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NOISE  MEASUREMENTS 


Initial  noise  measurements  are  expected  to  start  sometime 
in  October  1971.  The  time  required  to  collect  necessary  and  suf- 
ficient data  for  reasonable  statistical  modeling  and  simulating 
is  estimated  at  approximately  three  to  four  months.  Therefore, 
we  expect  the  measurements  will  be  completed  by  the  end  of  March 
1972.  Note  that  the  channel  measurement  will  be  conducted  con- 
currently, though  the  channel  measurements  and  noise  measurements 
are  not  being  made  simultaneously.  The  estimated  total  time  has 
taken  into  consideration  winter  weather  conditions  in  the  Boston 
area . 

CHANNEL  MEASUREMENTS 

The  time  schedule  for  channel  characterization  measurements 
is  the  same  as  that  for  the  noise  measurement. 

DATA  ANALYSIS  AND  MODELING 

Data  analysis  will  begin  as  soon  as  the  first  measurements 
are  made.  The  data  analysis  will  interact  with  actual  data  col- 
lection as  to  the  location,  time,  format  of  the  data,  etc.  The 
initial  data  analysis  will  emphasize  first  order  statistics  such 
as  pulse  amplitude  distribution,  pulse  duration  distribution, 
pulse  separation  distribution  and  rms  values.  Further  analysis 
will  involve  second  order  statistics  such  as  correlation  with 
traffic  flow,  zoning  variations,  time  of  the  day,  type  of  the 
station  terminals  and  others.  The  period  of  the  data  analysis 
will  carry  to  the  end  of  FY'72. 

The  major  effort  of  modeling  will  begin  after  the  prelimin- 
ary analysis  has  been  completed.  Because  we  have  reviewed  most 
of  the  modeling  work  which  has  been  done  in  the  past,  it  is  our 
belief  that  any  further  work. on  modeling  will  be  meaningless  un- 
til a certain  amount  of  data  has  been  collected  and  analyzed 
using  the  new  measuring  equipment  which  is  substantially  differ- 
ent than  any  used  before  in  studying  the  urban  channel.  The 
major  differences  include  wideband  techniques,  large  dynamic 
range,  and  different  frequency  ranges. 

SIMULATOR 

The  contractor,  Signatron,  Inc.,  has  the  chief  responsibil- 
ity for  designing  and  building  a channel  analyzer  for  the  fading 
dispersive  channel  with  provision  for  including  the  additive 
noise.  The  delivery  date  of  the  simulator  is  approximately  the 
end  of  1971  (calendar  year)  . During  the  interim,  TSC  personnel 
will  keep  close  contact  with  Signatron,  Inc.  personnel  to  insure 
that  contract  requirements  are  met  according  to  schedule.  During 
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the  ensuing  six  months,  TSC  personnel  will  be  trained  to  test  and 
operate  the  simulator.  The  tape  recorded  data  for  both  channel 
and  noise  will  be  used  for  evaluation  of  existing  equipment  on 
an  equal  basis  and  for  the  design  and  evaluation  of  new  signal- 
ling methods  and  new  command  and  control  equipment. 

APPENDIX  A 

INTRODUCTION 

This  paper  describes  the  noise  measurement  I/O  Controller 
for  the  PDP-11. 

Data  in  the  sections  Timing  and  Measurement  Accuracy  should 
be  sufficient  to  determine  the  validity  of  the  measurements  and 
to  plan  algorithms  for  data  evaluation  and  analysis.  Programming 
provides  information  for  writing  programs  which  use  the  control- 
ler. The  section  Controller  Operation  describes  the  controller 
hardware  and  operation. 

In  the  future  this  document  will  be  expanded  to  include 
logic  diagrams,  maintain  information,  and  sample  programs. 

This  is  a preliminary  description  and  instruction  manual 
for  the  noise  measurement  I/O  controller.  All  specifications 
and  parameters  are  preliminary  and  subject  to  change  during  de- 
vice construction. 

DESCRIPTION  OF  I/O  CONTROLLER  AND  MEASUREMENT  PROCESS 

The  I/O  controller  is  a direct  memory  access  interface  to 
the  PDP-11.  It  is  connected  to  the  PDP-11  Unibus.  The  device 
has  been  designed  to  measure  the  significant  parameters  of  an 
impulse  noise  random  process  and  transfer  them  into  the  PDP-11 
memory . 

Three  noise  process  parameters  are  measured  for  each  noise 
impulse;  elapsed  time  since  last  pulse,  duration  of  pulse,  and 
amplitude  of  pulse.  These  parameters  are  measured,  and  then 
transferred  directly  to  the  PDP-11  memory  without  any  program 
interrupt.  The  program  is  interrupted  only  when  a buffer  is 
filled  with  three  word  records,  one  record  from  each  noise  im- 
pulse. Certain  error  conditions  can  also  cause  an  interrupt; 
these  will  be  described  later. 

TIMING 


The  ability  of  the  device  to  generate  a representation  of  a 
noise  process  is  strongly  dependent  on  its  timing  characteris- 
tics. These  are  in  turn  dependent  on  the  computer  core  memory 
organization  and  Unibus  timing.  Complete  understanding  of  the 
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remainder  of  this  section  requires  an  understanding  of  the  PDP-11 
Unibus  and  possible  memory  organizations;  for  this  information 
see  the  PDP-11  manuals.  A flow  chart  of  the  I/O  controller  tim- 
ing and  control  sequence  is  shown  in  Figure  7.  It  should  be 
noted  that  timing  is  relatively  complex  and  will  typically  be 
dependent  on  previous  measured  data.  The  A to  D conversion  time 
tAD  will  be  11  ys  (approximately)  in  initial  operation  of  the 
controller  with  the  DEC  AD01  analog  to  digital  converter.  This 
will  be  the  limiting  factor  for  controller  operation.  The  A to 
D converter  can  be  replaced  by  an  8 bit  1 ys  converter  at  a rea- 
sonable cost  in  which  case  memory  write  time  (see  Figure  8)  will 
become  the  principal  speed  limitation.  In  the  initial  controller 
configuration  a time  of  16  ys  to  process  each  impulse  would  be  a 
good  estimate.  This  time  can  be  reduced  to  4 ys. 

MEASUREMENT  ACCURACY 

The  peak  amplitude  measurement  error  comes  from  two  sources: 
A to  D converter  error  and  peak  value  sampler  errors.  The  peak 
value  is  sampled  by  waiting  a preset  time  interval  after  the 
crossing  of  a threshold.  This  is  based  on  the  assumption  that 
the  noise  pulse  has  a much  broader  bandwidth  than  the  measure- 
ment system  and  that  the  rise  time  of  the  sampler  input  is  pro- 
portional only  to  measurement  system  bandwidth.  If  measurement 
data  fails  to  confirm  this  approximation  a more  complex  peak  de- 
tector can  easily  be  installed.  A to  D converter  accuracy  is  10 
bits  which  is  far  in  excess  of  that  required  for  this  experiment. 
The  sample-and-hold  module  supplied  with  the  AH01  does  not  have 
sufficient  bandwidth  to  allow  the  .01%  accuracy  of  the  A to  D 
converter  to  be  achieved.  This  error  source  will  have  to  be  ex- 
amined during  device  calibration. 

The  time  measurements  are  made  by  counting  the  pulses  gen- 
erated by  a clock  in  the  I/O  controller.  The  counter  cannot 
measure  durations  of  less  than  200  ns  but  this  will  not  be  a de- 
tectable limitation  on  system  performance.  The  clock  pulse  rate 
can  be  up  to  25  MHz  but  it  will  initially  be  set  to  4 MHz  (250  ns 
resolution) . The  counter  will  round  off  measured  times  to  the 
next  lowest  integer  time  count.  All  time  intervals  of  less  than 
200  ns  will  be  indicated  as  200  ns;  counter  error  will  be  plus 
or  minus  one  half  count  plus  clock  accuracy  for  times  up  to  4096 
counts.  For  times  above  4096  counts  (initially  1024  ys)  accur- 
acy is  plus  or  minus  .05%  plus  clock  error. 

The  time  measurements  require  a large  dynamic  range.  To 
suppress  overflow  problems  an  unsigned  floating  point  number  rep- 
resentation was  selected.  This  results  in  a word  format  as  shown 
in  Figure  9.  The  largest  count  before  overflow  is  215*  ( 2 1 2 -1 ) = 
134  184  960.  This  will  initially  correspond  to  33.5  seconds. 

The  floating  point  representation  will  require  no  arithmetic 
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Figure  7.  Measurement  Action  Sequence  and  Timing  Data 
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ABBREVIATIONS  ARE: 

AC  = ADDRESS  COUNTER 

WC  = WORD  COUNTER 

t = PROCESSING  TIME  FOR  BLOCK 

t = MEMORY  WRITE  TIME 
in 

t,  = UNIBUS  DELAY  TIME 
b 


Figure  8.  Data  Transfer  Operations  for  Timing  Estimation 
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BIT 


BIT 


15  14  13  12 

11  10  9876543210 

DATA  WORD  FROM  COUNTER 

EXPONENT 

UNSIGNED 


0 < EXP  < 15 


FRACTION 


0 IF  < 2048 

1 IF  > 2048 
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UNNORMALIZED  IF  < 2048 
NORMALIZED  IF  > 2048 


0 < FRACTION  < 4095 


NUMBER  = FRACTION  * 2 
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15 

14  13  12  11  10  9 8 7 6 5 

4 3 2 1 0 

DATA  WORD  FROM  A TO  D CONVERTER 

0 

MOST  LEAST 

SIGNIFICANT  SIGNIFICANT 
DIGITS  DIGITS 

INSIGNI- 

FICANT 

DIGITS 

0 1 1 1 1 

Figure  9.  Data  Word  Formats 
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operations  for  conversion  to  integer  format  when  the  count  is 
less  than  4096*  This  will  be  the  situation  most  of  the  time. 

GROWTH  POTENTIAL 

Provision  has  been  made  to  transfer  an  extra  data  word  in 
the  I/O  controller  if  an  additional  parameter  must  be  measured. 
This  will  require  significant  controller  modifications.  The 
time  counters  cannot  easily  be  modified,  but  by  changing  the 
clock  frequency  a large  range  of  resolutions  and  accuracies  can 
be  provided.  The  A to  D converter  can  be  easily  changed  or  re- 
placed by  another  measurement  device.  The  pulse  detector  and 
peak  detector  circuits  can  be  easily  modified. 

Two  dc  logic  levels  which  can  be  set  by  the  program  are 
available  for  setting  modes  or  gains. 

A pulse  is  generated  each  time  input  noise  impulses  occur 
too  closely  spaced  in  time  for  the  controller  to  measure  one  of 
them.  These  pulses  can  be  counted  by  an  external  counter  or  the 
controller  can  be  modified  to  interrupt  the  CPU  when  this  condi- 
tion occurs. 

PHYSICAL  DESCRIPTION 

The  I/O  controller  is  mounted  on  one  DEC  5h  inch  high  cir- 
cuit card  rack.  Its  power  supply  occupies  3 \ inches  of  rack 
space.  The  controller  should  be  mounted  next  to  the  AD-l-D 
Analog  to  Digital  Converter  in  the  computer  rack.  Total  power 
consumption  is  less  than  70  watts.  No  power  is  drawn  from  the 
computer . 

PROGRAMMING 

The  controller  contains  three  registers  which  are  address- 
able by  the  program.  They  are:  Address  Count  (AC) , Word  Count 

(WC) , and  status  (S) . These  registers  are  described  below. 

Word  Count  - Location  767774,  write  only 

This  register  contains  the  buffer  size  in  words.  It  need 
only  be  loaded  when  the  controller  is  initialized  or  buffer  size 
is  changed. 

Address  Count  - Location  767772,  read/write 

This  register  contains  the  address  to  which  the  next  data 
word  will  be  transferred.  At  the  end  of  a buffer  the  address 
will  be  one  beyond  the  last  word  of  the  buffer.  If  an  error  in- 
terrupt has  occurred  this  value  tells  how  many  good  words  of  data 
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have  been  transferred.  Loading  this  register  always  causes  the 
device  status  to  be  set  to  "run  enabled,"  and  resets  Word  Count 
to  its  initial  value  (the  buffer  size) . 

Status  - Location  767770 

The  bits  of  the  status  register  are: 

Int  Enb; 

Bit  6, 

Read/Write , 

Initial  state  0, 

Function:  enables  end  of  buffer  and  error  interrupts 

Run; 


Bit  7 , 

Read  only, 

Set  by  loading  AC, 

Reset  on  any  error  or  end  of  buffer, 

Function:  enables  data  transfer 

Initial  state  0 

Time  Out; 

Bit  15, 

Read  only. 

Set  by  hardware  or  address  range  error, 

Reset  by  loading  AC, 

Initial  state  0. 

In  addition  to  these  three  bits  other  bits  are  described  in 
the  maintenance  section.  Status  is  set  to  its  initial  condition 
by  the  PDP-11  power-on  sequence  or  by  the  RESET  instruction  (re- 
sets all  devices) . 

In  addition  to  these  registers  the  Analog  to  Digital  Con- 
verter status  register  must  be  set  to  000030.  The  status  regis- 
ter is  location  776070.  This  location  is  reset  to  000000  by 
power  on  or  RESET.  The  word  000030  sets  the  converter  gain  to 
2.5  v full  scale. 

Interrupt  Vectors 

At  the  end  of  a buffer  the  interrupt  vector  is  location 
000170.  If  any  error  occurs  the  interrupt  vector  will  be  000174. 
If  both  error  and  end  of  buffer  occur  simultaneously  the  vector 
is  000174.  If  an  error  occurs  the  value  of  ac  should  be  saved 
and  then  ac  set  to  a new  buffer.  The  old  AC  value  can  then  be 
used  immediately  to  recover  all  good  data. 
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CONTROLLER  OPERATION 

Overview 

The  Unibus  is  described  in  the  PDP-11  Handbook.  A function- 
al block  diagram  of  the  controller  is  shown  in  Figure  10.  The 
measurement  and  data  transfer  parts  of  the  controller  are  asyn- 
chronous and  essentially  independent  of  each  other.  The  indi- 
vidual blocks  in  Figure  10  are  described  below. 

Measurement  Control  Logic 


The  measurement  control  logic  consists  of  a precision 
Schmidtt  trigger  circuit  which  acts  as  the  pulse  detector  fol- 
lowed by  timing  and  sequencing  circuits  to  control  the  measure- 
ment. A block  diagram  of  this  circuit  is  shown  in  Figure  11. 

Measurement  Devices 


Performance  of  the  floating  point  counter  was  described  in 
the  chapter  Technical  Discussion.  The  counter  can  count  at  a 
25  MHz  maximum  rate  but  it  takes  approximately  200  ns  for  the 
count  to  propagate  to  the  end  of  the  count  chain.  When  the  count 
period  ends,  a one  shot  times  out  the  200  ns  propagation  delay. 
The  counter  contents  are  then  transferred  into  one  of  two  buffer 
registers  (pulse  duration  or  pulse  spacing)  and  the  counter  is 
reset  to  the  value  that  corresponds  to  the  elapsed  transfer  time 
(200  ns  approximately) . The  counter  is  always  counting  except 
when  a data  transfer  occurs.  A counter  block  diagram  is  shown 
in  Figure  12. 

The  Analog  to  Digital  converter  is  a standard  DEC  AD01-D;it 
is  triggered  through  its  external  input  (EXT  IN  A) . 

Address  Counter 


The  address  counter  and  logic  contains  three  registers;  one 
for  current  address,  one  for  current  word  count,  and  one  for  buf- 
fer size.  When  the  controller  is  to  be  operated,  the  buffer  size 
register  (WC)  must  be  loaded  first.  When  the  current  address 
counter  is  loaded  from  the  Unibus,  the  current  word  - counter  is 
loaded  from  the  buffer  size  register.  Each  single  word  data 
transfer  results  in  the  address  count  being  incremented  by  two 
and  the  word  count  being  decremented  by  two.  End  of  buffer  is 
signaled  when  the  word  count  reaches  zero. 

This  structure  was  selected  to  allow  buffers  to  be  swapped 
with  one  CPU  instruction.  Any  buffer  size  or  buffer  location  can 
be  used;  additional  instructions  are  necessary  only  if  errors  oc- 
cur or  if  buffer  size  must  be  changed.  The  Address  Counter  is 
described  in  Figure  13. 
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Figure  11.  Measurement  Control  Logic  Block  Diagram 
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Figure  12.  Floating  Point  Counter  Block  Diagr 
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Data  Transfer  Control  Logic 


The  data  transfer  control  logic  is  relatively  standard.  It 
is  built  around  the  DEC  M782  Interrupt  Control  Module  and  the 
DEC  M104  Address  Decoder  Module.  The  circuit  is  described  in 
detail  in  the  maintenance  section. 
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REPRINT  OF  IN-HOUSE  PUBLICATIONS 


URBAN  EM  NOISE  MEASUREMENT,  Jan  25,  1971,  P.  Yoh  and  R.  Esposito 
INTRODUCTION 

During  the  past  several  decades  many  measurements  of  man- 
made electromagnetic  noise  have  been  made  at  HF,  VHF  and  UHF  for 
different  environments  such  as  urban,  suburban  and  rural  areas. 
However,  because  of  the  instrumentation  used  and  of  the  limited 
scope  of  most  of  the  measurements,  the  results  are  often  of 
limited  use  and  sometimes  misleading  for  the  communication  en- 
gineer. Basically,  a communication  system  planner  would  like  to 
know  the  detailed  characteristics  of  the  background  noise  in 
order  to  evaluate  its  effects  on  system  performance.  This  in- 
formation includes  also  the  time  and  space  variations  of  the 
noise  as  well  as,  for  prediction  purposes,  the  correlation  with 
factors  such  as  population  density,  vehicular  traffic,  power 
consumption  and  topographical  characteristics.  The  requirements 
are  more  stringent  for  digital  communication  systems  which,  for 
example,  demand  a more  thorough  specification  of  the  statistical 
characteristics  of  the  noise.  By  and  large  this  information  is 
not  currently  available  especially  for  wideband  systems,  i.e., 
for  systems  whose  bandwidth  exceeds  10  KHz.  Therefore,  the 
purpose  of  the  experiment  under  consideration  is  to  measure 
broadband,  wide  dynamic  range,  impulsive-type  noise  in  selected 
regions  of  the  VHF  and  UHF  band.  The  bandwidth  of  the  measuring 
system  has  been  tentatively  selected  to  be  1 MHz.  The  predomi- 
nant noise  is  expected  to  be  generated  by  man-made  machines  such 
as  automobile  engines,  industrial  equipment,  etc. 

The  reasons  for  this  experimental  program  are  the 

following: 

1.  To  study  experimental  techniques  for  the  measurement  of 
wideband  man-made  noise  in  the  VHF-UHF  region  and 
establish  a capability  for  the  routine  performance  of 
such  measurements. 

2.  To  identify  the  major  sources  of  radiation  that  will 
interfere  with  the  performance  of  any  transportation 
system  being  planned  in  the  frequency  region  mentioned 
above. 

3.  To  determine  the  information  needed  for  the  design  and 
evaluation  of  any  communication  or  command  and  control 
system  of  mass  transit  or  land,  air,  water  and  highway 
communication.  The  combination  with  the  measurements 
of  the  impulse  response  of  the  medium  this  information 
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will  provide  a complete  characterization  of  the  urban 
and  suburban  channel. 

4.  To  determine  the  radiation  characteristics  of  the 

existing  and  proposed  transportation  systems  such  as 
electric  cars,  trolley  cars,  guided  way,  electric 
trains,  etc.  Figure  1 indicates  the  major  tasks  and 
problem  areas  and  the  logical  flow  of  data  of  the 
measurement  program. 

NOISE  MEASUREMENTS 

One  of  the  first  tasks  is,  of  course,  the  determination  of 
the  statistical  quantities  describing  the  noise  which  can  be 
derived  by  suitable  processing  of  the  experimental  data.  Al- 
though a certain  amount  of  data  processing  can  be  done  in  real 
time  the  decision  has  been  made  to  record  the  IF  waveform  for 
later  processing  and  for  use  in  the  simulation  equipment.  This 
means  that  most  of  the  processing  procedures  can  be  decided  upon 
at  a later  date.  Some  initial  considerations  on  what  has  to  be 
measured  is,  however,  in  order.  A simple  model  of  the  wideband 
man-made  interference  is  the  sum  of  an  impulsive  component  and 
of  a Gaussian  process.  Both  components  can  be  expected  to  be 
nonstationary.  It  seems  necessary  to  measure  at  least  the 
following  quantities : 

1.  Average  power  over  each  measurement  interval 

2.  Percentage  of  power  in  the  impulsive  component  and 
Gaussian  component 

3.  Amplitude  pdf  of  the  peaks 

4.  Duration  pdf  of  the  impulses  at  given  amplitude  levels 

5.  Pulse  spacing  pdf 

6.  Crossing  rate  pdf. 

In  addition  to  the  above  quantities  some  simpler  parameters 
such  as  average  and  rms  values  after  detection  can  be  measured 
in  real  time.  It  is  worth  pointing  out  that  instrumentation 
similar  to  the  one  being  studied  for  the  noise  measurements  can 
be  used  in  the  future  for  the  investigation  of  channel  occupancy. 
By  this  we  mean  the  evaluation  of  the  probability  that  a given 
frequency  region  be  occupied  at  a certain  level,  by  other  trans- 
mitters. Because  of  the  limited  spectral  region  available  for 
land  mobile  communication  it  is  also  important  to  determine  the 
minimum  allowed  spatial  separation  of  stations  transmitting  in 
the  same  frequency  region  in  order  to  avoid  appreciable  co-channel 
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Figure 


interference.  This  type  of  measurements  can  also  be  performed 
by  the  system  under  consideration. 

A simple  block  diagram  of  the  noise  measuring  system  is 
shown  in  Figure  2 . 

The  following  is  a preliminary  discussion  of  the  specifica- 
tions for  the  various  blocks  of  the  system  and  of  the  possible 
tradeoffs . 

Frequencies  of  Operation 


The  measurements  will  be  taken  in  the  neighborhood  of  150, 
450  and  900  MHz.  Inquiries  have  been  made  to  the  FAA  to  deter- 
mine whether  a 1 MHz  slot  free  from  intelligible  interference  is 
available  in  the  greater  Boston  area  in  the  frequency  regions 
mentioned. 

Antenna 


In  deciding  on  the  type  of  antenna  to  be  used  a number  of 
questions  should  be  considered  such  as  polarization,  antenna 
pattern,  gain,  frequency  response  characteristics  and  mutual 
coupling  if  there  are  other  antennas  close  by.  All  of  the  above 
factors  depend  on  the  type  of  measurements.  It  seems  that  for 
the  spatial  distribution  expected  from  urban  noise  under  differ- 
ent conditions  the  omnidirectional  short  dipole  will  be  a good 
candidate.  A directional  antenna  would,  instead,  be  necessary 
if  one  wants  to  establish  the  location  of  the  noise  sources. 

Also  antennas  with  different  polarization  are  needed  if  polari- 
zation effects  must  be  studied.  Note  that  some  of  the  measure- 
ments available  indicate  little  difference  between  horizontal 
and  vertical  polarization.  A number  of  closely  spaced  antennas 
would  be  needed  if  one  wants  to  map  the  spatial  distribution  of 
the  noise  sources.  Thus,  the  initial  specifications  for  the 
antennas  are: 

Frequency:  150  , 450  and  900  MHz  within  + 10% 

Polarization:  vertical  over  ground  plane 

Antenna  Pattern:  Omnidirectional  in  azimuth  and  cos  0 in 

elevation,  where  0 is  evaluation  angle 


Type : 


Short  dipole 


Gain : 


3 dB 
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Figure 


Considerations  for  future  extensions  of  the  experiment  are: 

1.  Simultaneous  measurements  of  different  polarizations 

2.  Simultaneous  measurements  of  all  three  frequency  bands 

3.  Array  measurements  at  the  same  frequency,  including  the 
study  of  mutual  coupling  effects 

4.  Change  of  antenna  pattern 
RECEIVER  ~ RF,  Mixer,  L.O.  and  IF  State 

The  two  essential  quantities  in  the  overall  design  of  the 
receiver  are  sensitivity  and  dynamic  range.  These  quantities 
are  related  since  a high  sensitivity  requires  a large  dynamic 
range  if  the  external  noise  has  large  variations.  The  system 
sensitivity  is  defined  as  the  level  at  which  external  noise 
power  equals  the  receiver  noise  power.  The  external  noise  has 
both  natural  and  man-made  origin.  The  natural  noise  is  galactic, 
atmospheric  and  solar.  However,  in  the  frequency  range  of 
interest  this  noise  is  almost  negligible  with  the  exception  of 
periods  of  high  solar  activity.  The  total  contribution  of  these 
sources  can  be  3 dB  above  the  receiver  thermal  noise  level. 

Thus,  if  low  level  man-made  noise  is  at  150  MHz  and  negligible 
above  this  frequency  of  interest,  the  sensitivity  of  the  re- 
ceiver should  be  as  good  as  possible.  Since  the  thermal  noise 
level  corresponding  to  thermal  noise  for  a reference  temperature 
of  288  K and  1 MHz  bandwidth  is  144  dBW  (or  -114  dBM)  the  re- 
ceiver must  be  able  to  detect  any  signal  above  a level  of  about 
-140  dBW.  For  an  input  resistance  of  50  ohms  this  implies  a 
detectable  noise  level  of  0.7  yV.  The  minimum  detectable  noise 
field  can  be  easily  determined  by  using  the  capture  area  of  the 
antenna. 
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where  A is  the  wavelength. 

The  previous  considerations  suggest  that  the  noise  figure 
of  the  RF  stage  should  be  rather  low,  of  the  order  of  3 or  4 dB 
if  possible.  It  seems  also  desirable  to  make  the  bandwidth  of 
the  RF  stage  the  same  as  the  IF  bandwidth  in  order  to  prevent 
overloading  of  the  mixer.  The  RF  stage  should  also  have  some 
gain  (10  dB  seems  a good  ballpark  number)  in  order  to  prevent 
overall  sensitivity  degradation  caused  by  mixer  and  local 
oscillator  noise.  The  overall  sensitivity  of  the  receiver  should 
be  1 yV/m  or  better. 
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The  next  question  is  the  determination  of  the  dynamic  range 
required.  This  is  perhaps  the  most  difficult  parameter  to 
determine  since  it  involves  an  estimate  of  the  unknown  amplitude 
distribution  of  the  man-made  noise.  Note  that  the  dynamic  range 
can  always  be  expressed  in  percentiles  of  the  amplitude  distri- 
bution. For  example,  ITS  measurements  of  amplitude  probability 
density  of  man-made  noise  made  in  San  Antonio,  Texas  at  50  MHz 
with  a bandwidth  of  10  KHz  indicate  a dynamic  range  of  86  dB 
between  the  10~6  and  0.99  probability  levels.  It  is  often 
stated  that  the  dynamic  range  varies  roughly  as  10  logqo  (band- 
width) and  therefore,  the  dynamic  range  for  1 MHz  should  be  106 
dB  (with  the  further  assumption  that  the  dynamic  range  is  in- 
dependent of  the  frequency  of  operation) . However,  a closer 
look  at  the  argument  which  yields  the  10  logqo  b relation  shows 
that  it  is  not  valid  when  the  bandwidth  is  very  large  if  the 
number  of  incoming  impulses  is  not  exceedingly  large.  In  fact, 
one  can  argue  that  a wideband  receiver  can  resolve  the  individual 
pulses  and  only  very  rarely  will  the  impulses  have  enough  over- 
lap to  cause  a significant  increase  in  the  required  dynamic 
range.  Since  each  car  radiates  a few  hundred  impulses/sec  and 
only  very  close  by  cars  contribute  to  the  upper  part  of  the 
dynamic  range  the  argument  seems  reasonable.  Then,  what  is 
important  for  the  determination  of  the  dynamic  range  is  the 
maximum  peak  value  of  impulses  received  from  a car  located  very 
close  to  the  receiving  equipment.  Some  measurements  made  in 
England  with  a 2.5  MHz  receiver  indicate  a maximum  received 
field  of  56  mV/m  at  150  MHz  when  the  radiating  car  is  30  feet 
away  from  the  receiver.  The  maximum  received  field  is  about 
5 dB  higher  at  450  MHz.  This  value,  converted  to  a 1 MHz  band- 
width, yields  a value  of  about  80  dB  of  dynamic  range  based  on 
an  overall  system  noise  figure  of  7 dB  at  150  MHz  and  a slightly 
lower  value  at  450  MHz.  Since  we  are  interested  in  events  of 
as  low  a probability  as  10 ”6,  the  effective  dynamic  range  must 
be  increased  to  take  into  account  the  possibility  of  overlapping 
of  impulses  coming  from  different  cars.  This  increase  can  be 
estimated  to  be  of  the  order  of  10  dB.  Therefore,  the  overall 
dynamic  range  seems  to  be  of  the  order  of  90  dB. 

It  is  our  intention  to  perform  some  simple  measurements 
with  instrumentation  available  in  the  laboratory  to  ascertain 
whether  our  estimate  of  the  dynamic  range  required  is  correct. 

Recording  and  Data  Processing 


As  mentioned  before  part  of  the  processing  will  be  done  in 
real  time  and  part  by  computer.  This  tradeoff  is  very  important 
because,  owing  to  the  large  dynamic  range  and  bandwidth  used  in 
the  experiment,  some  severe  requirements  are  imposed  on  the 
recording  and  data  processing  system.  To  show  this,  assume  a 
dynamic  range  of  90  dB  and  a 1 MHz  bandwidth.  In  order  to  use 
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digital  processing  the  data  must  be  sampled  at  a rate  of  2 MHz, 
while  10  bits  are  required  for  10%  accuracy  (1  bit  is  the  sign) . 
Thus  the  measuring  system  generates  data  at  a rate  of  20  Mbit/sec. 
Assume  now  that  each  processing  operation  takes  30  computer 
cycles  and  that  the  computer  cycle  is  1 ysec.  Assume  also  that 
the  tape  is  played  back  into  the  computer  at  a speed  4 times 
lower  than  the  recording  speed.  For  100  sec  of  data  the  required 
processing  time  for  each  operation  is  thus 

20.10^*30.10  ^ x 100  x 4=2.4.10^  sec  = 66  hrs , 

a large  number  indeed.  One  would  still  require  over  16  hours  if 
the  tape  were  to  be  played  back  at  the  recording  speed.  Since 
most  of  the  measurements  we  want  to  make,  at  least  initially, 
probability  density  functions,  the  use  of  high  speed  threshold 
counters,  operating  over  selected  time  intervals  and  energy 
levels  must  be  considered  for  real  time  processing. 

A simple  example  of  the  possible  real  time  processor  is 
shown  below  in  Figure  3. 

IF  AMP 


Figure  3 . 


Recorder 
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Calibration 


This  system  shall  have  the  necessary  equipment  to  accu- 
rately and  conveniently  calibrate  all  channels  for  both  noise 
and  channel  measurement,  so  hence  the  external  noise  or  signal 
relative  to  the  terminal  of  an  equivalent  lossless  antenna 
system  is  measured  in  absolute  units  to  accuracies  of  + 3 dB. 
The  receiver,  detectors,  and  A-D  converter  output  shall  be 
calibrated  to  measure  to  accuracies  of  + 1 dB. 

The  absolute  calibration  on  the  channel  measurement  is 
needed  so  that  the  simulation  study  can  be  conducted  success- 
fully . 

Operation 

The  equipment  must  be  grouped  together  in  a console  so 
that  it  may  be  monitored  by  a single  seated  operator  who  can 
conveniently  listen  to  the  noise  channel  and  make  note  on  dis- 
crete signals  when  they  occur.  The  operator  must  be  able  to 
observe  the  system  outputs  to  ascertain  the  system  is  function- 
ing properly.  In  case  of  malfunction  or  others  the  operator 
should  try  to  adjust,  eliminate  or  note. 

Procedure  of  Measurement 

It  is  tentatively  planned  that,  initially,  a single  fre- 
quency measurement  will  be  made  in  the  greater  Boston  area. 
Because  of  the  large  variety  of  the  zones:  commercial,  in- 

dustrial, residential,  transportation  terminals,  high  rise 
buildings,  etc.,  it  presents  a great  opportunity  to  measure  the 
noise  characteristics  of  a city.  However,  a single  measurement 
can  be  made  at  any  given  time.  Hence,  it  is  strongly  suggested 
to  have  a detailed  plan  on  when  and  where  the  measurements 
shall  be  made.  The  following  criteria  may  be  useful  for  the 
planning : 

1.  Vehicle  volume 

2.  Population  density 

3.  Geometrical  distribution  of  high  rise  building 

4.  Geographical  and  topological  features  of  the  city 

5.  Zoning  characteristics 
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Abstract . - A program  to  measure  VHF  and  UHF  urban  radio 
channels  for  playback  and  analysis  in  the  laboratory  is  des- 
cribed. The  measured  data  will  be  used  to  develop  analytic 
models  for  channel  behavior.  A 50  KHz  bandwidth  tapped  delay 
line  real-time  channel  simulator  is  being  built.  Both  the 
simulator  and  the  models  will  be  used  for  the  design  and  evalua- 
tion of  systems  for  communication  and  vehicular  monitoring  in 
the  urban  environment. 

Presented  at  1971  International  Telemetry  Conference, 
Washington,  D.C.  Proc.  International  Telemetering  Conf.  Vol. 

VII  pp.  166-171,  International  Foundation  for  Telemetering, 
Woodland  Hill,  California,  Reprinted  by  permission  of  the  Instru- 
ment Society  of  America,  publisher  of  the  proceedings. 
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The  authors  are  at  the  U.S.  Dept,  of  Transportation,  Transporta- 
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SUMMARY 

This  paper  outlines  a program  to  measure  VHF  and  UHF  urban 
radio  channels  for  playback  and  analysis  in  the  laboratory. 

These  measurements  are  required  for  the  design  of  high  data  rate 
communication  equipment  which  makes  efficient  use  of  the  radio 
spectrum.  We  include  in  this  category  data  communication  and 
surveillance  equipment  necessary  to  meet  the  demands  of  the 
variety  of  command  and  control  functions  envisioned  for  trans- 
portation systems  by  federal,  state,  and  local  governments. 
Included  in  the  paper  is  a qualitative  description  of  the  channel 
followed  by  a discussion  of  the  measurement  techniques  and  hard- 
ware that  will  be  used.  The  goal  of  the  measurement  program  is 
to  develop  both  analytic  and  experimental  tools  for  the  design 
and  evaluation  of  advanced  communication  techniques  specialized 
for  signalling  over  the  urban  radio  channel.  These  techniques 
will  include  modulation  waveform  design,  coding,  and  the  design 
of  receiver  structures  including  those  which  adapt  to  measured 
channel  behavior.  Future  technology  will  permit  complex  hardware 
designs  to  be  used  in  the  implementation  of  these  techniques 
while  maintaining  low  system  cost.  The  same  techniques  and 
approach  used  here  for  the  urban  channel  can  be  applied  to  other 
radio  channels  that  occur  in  transportation  systems.  In  partic- 
ular, they  are  applicable  to  data  links  for  air-ground  and  high 
speed  ground  transportation  system  uses. 

Introduction 

Requirements  for  communication  in  an  urban  environment  have 
become  much  more  stringent  in  recent  years.  Where  traditionally 
voice  channels  have  been  sufficient  to  satisfy  user  needs,  new 
requirements  have  emerged  for  data  communication  and  surveillance 
to  support  proposals  for  command  and  control  for  transportation 
systems  made  by  federal,  state  and  local  governments.  For  ex- 
ample, automatic  vehicle  monitoring  and  digital  data  links  have 
been  proposed  for  both  public  and  private  vehicles  such  as 
buses,  taxis,  private  automobiles,  and  special  services  including 
police  and  fire  (this  latter  category  may  also  require  some  form 
of  secure  communication) . Implicit  in  these  system  requirements 
are  demands  for  increased  communication  capacity  — high  data 
rates  and  efficient  utilization  of  the  rapidly  disappearing 
available  radio  spectrum.  In  order  to  provide  both  the  reliable 
and  efficient  communication  needed,  one  must  understand  in  detail 
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the  urban  communication  environment.  The  measurement  of  the 
properties  of  this  radio  channel  is,  therefore,  not  academic. 

It  is  integral  to  the  design  and  evaluation  of  communication 
equipment  matched  to  the  job  in  both  performance  and  cost.  The 
urban  radio  channel  environment  will  be  qualitatively  described 
and  a program  will  be  outlined  which  will  systematically  measure 
and  playback  the  channel  for  communication  system  design. 

PROGRAM  DESCRIPTION 

There  are  two  basic  technological  problem  areas  in  urban 
radio  communication.  The  available  spectral  space  is  severely 
limited  and  will  not  grow  appreciably  in  the  future.  At  the 
same  time  both  the  number  of  users  and  the  type  and  variety  of 
services  required  are  growing  with  time.  The  radio  propagation 
medium  in  urban  areas  is  characterized  by  severe  multipath  dis- 
tortion and  impulsive  noise  which  makes  the  design  of  high  per- 
formance non-voice  communication  systems  difficult. 

Short-term  solutions  for  increasing  radio  traffic  handling 
capacity  are  to  open  new  frequency  bands,  use  narrower  band- 
width channels,  and  to  promote  good  operating  procedures  for 
minimizing  radio  channel  usage  by  individual  users.  These  tech- 
niques will  not  provide  the  several  orders  of  magnitude  increase 
in  traffic  handling  capacity  that  will  be  needed  over  longer 
time  periods.  The  only  efficient  solution  to  this  problem 
appears  to  be  the  use  of  semi  or  fully  automatic  digital  data 
transmission  with  prespecified  message  formats  and  relatively 
complex  polling  or  scheduling  procedures  for  control  of  channel 
usage.  The  efficient  transmission  of  digital  data  on  urban 
radio  channels  is  difficult  because  such  channels  exhibit  com- 
plex multipath  propagation,  typically  without  any  direct  path 
between  transmitter  and  receiver.  This  means  that  for  reliable 
digital  data  transmission,  signaling  waveforms  and  modem  struc- 
tures must  be  carefully  designed  to  match  specific  channel 
behavior.  Some  sort  of  coding  will  almost  certainly  be  required. 
A detailed  understanding  of  urban  multipath  propagation  is  also 
required  for  design  or  evaluation  of  surveillance  systems  which 
locate  vehicles  by  measurement  of  radio  propagation  time.  In 
addition  to  the  multipath  problem,  radio  noise  in  cities  is 
especially  complex  and  highly  structured.  Receivers  designed 
for  white  gaussian  noise  will  not  perform  efficiently.  By  taking 
advantage  of  the  urban  noise  structure,  performance  will  be 
improved . 

There  are  a variety  of  techniques  for  designing  communica- 
tion systems  which  adapt  to  the  time  varying  behavior  of  the 
urban  radio  channel.  These  range  from  adaptive  equalizing 
channel  measuring  communication  systems  to  systems  which  allo- 
cate channel  usage  time  to  individual  users  according  to  their 
measured  performance. 
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In  an  adaptive  equalizing  system,  the  transfer  function  of 
the  channel  over  the  frequency  range  of  interest  must  be  measured 
and  an  inverse  filter  generated  at  either  end  of  the  link.  Adap- 
tive equalizing  systems  have  the  ability  to  allow  high  data  rate 
transmission  on  badly  dispersive  channels.  They  are  not  a total 
cure  for  time  dispersive  channel  problems  in  that  they  clearly 
can  not  increase  the  signal  energy  transferred  over  a fading 
channel.  This  in  effect  limits  their  usefulness  to  relatively 
broad  bandwidth  systems  in  which  the  signalling  waveforms  have 
sufficient  range  (or  equivalently  time)  resolution  to  allow  some 
of  the  multipath  scattering  structure  of  the  channel  to  become 
measurable.  Adaptive  equalizing  systems  can  only  be  designed  if 
a broad  bandwidth  characterization  of  the  multipath  channel  is 
available  for  the  required  analysis.  The  construction  of  such 
a system  is  made  much  easier  by  the  ability  to  simulate  the  chan- 
nel in  the  laboratory.  There  are  a variety  of  adaptive  equaliz- 
ing concepts  using  different  schemes  for  channel  multipath 
measurement,  and  different  schemes  for  apportioning  hardware  be- 
tween transmitter  and  receiver.  Several  different  techniques 
will  be  evaluated. 

Other  procedures  for  allowing  a vehicular  communication 
network  to  adapt  to  propagation  conditions  involve  dynamic  chan- 
nel scheduling  algorithms  which  attempt  to  take  account  of  the 
spatially  non-homogeneous  nature  of  the  urban  radio  channel. 

Also,  error-detecting  or  correcting  codes  with  two-way  data  links 
can  allow  a system  to  adapt  to  the  error  performance  of  individ- 
ual network  participants.  There  are  a large  variety  of  such 
systems  which  should  be  evaluated  to  determine  their  potential 
for  providing  reliable  and  economical  communication.  The  per- 
formance of  this  type  of  system  depends  on  the  time  variation  of 
channel  behavior.  Such  systems  can  be  expected  to  provide  good 
performance  if  reasonably  large  numbers  of  users  are  present.  A 
broad  bandwidth  channel  characterization  is  required  for  evalua- 
tion or  design  of  any  candidate  system  of  this  type. 

Error-detecting  or  correcting  codes  can  be  considered  for 
several  purposes  on  the  radio  channels  of  interest  here.  A few 
of  the  most  promising  uses  are:  control  of  errors  due  to  noise 

(a  classical  application) , control  of  errors  due  to  multiple 
users  interfering  with  each  other,  providing  redundant  informa- 
tion for  channel  measurement  while  transmitting  data,  and  facil- 
itating the  simple  generation  of  waveforms  with  good  properties 
for  transmission  over  the  multipath  channel  (e.g.,  spectral 
shaping) . The  channel  characterization  experiments  will  provide 
the  required  tools  for  evaluation  of  coding  techniques  for  these 
purposes . 

Security  from  both  interception  and  jamming  is  also  desired 
for  a variety  of  civil  communication  systems  (e.g.,  police). 
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It  is  much  easier  to  provide  secure  coding  for  digital  communi- 
cations than  it  is  for  voice  communications. 

In  order  to  design  and  evaluate  digital  communication 
systems  a detailed  understanding  of  the  urban  noise  and  propaga- 
tion environment  is  required.  At  present  a two-fold  program  for 
channel  characterization  and  measurement  is  planned.  An  attempt 
will  be  made  to  generate  analytic  models  for  channel  characteri- 
zation which  will  be  useful  for  system  design.  In  addition,  a 
real-time  hardware  simulator  is  being  developed  which  will  be 
used  to  play  back  the  channel  using  the  recorded  measurement 
data  permitting  the  design  and  evaluation  of  proposed  modems  in 
the  laboratory  (see  Figure  1) . The  next  phase  in  the  program 
will  be  to  evaluate  a variety  of  modulation  waveforms,  coding 
techniques,  and  receiver  structures  including  systems  which 
adapt  to  measured  channel  behavior.  The  channel  data  to  be 
measured  will  be  selected  to  provide  a realistic  test  of  modem 
performance. 

URBAN  RADIO  CHANNEL  DESCRIPTION 

The  urban  environment  is  composed  of  a variety  of  noise 
sources  (e.g.,  ignition  systems  of  motor  vehicles,  electrical 
machinery,  high  voltage  power  transformers,  etc.).  In  addition, 
the  physical  character  of  urban  areas  with  their  dense  construc- 
tion present  a variety  of  paths  for  transmitted  radio  signals  — ■ 
often  with  the  absence  of  a direct  path  from  transmitter  to 
receiver.  These  two  situations  and  their  interaction  qualita- 
tively describe  the  urban  radio  channel.  The  description  of 
the  channel  can  therefore  be  divided  into  two  parts:  1)  addi- 

tive noise  composed  of  high-level  impulsive  man-made  noise  and 
lower-level  background  noise;  2)  multipath,  a multiplicative 
noise  due  to  the  presence  of  multiple  signalling  paths  and 
causing  time-and- frequency  dispersion  of  transmitted  signals. 

MULTIPATH  CHANNEL  MEASUREMENT 

A characterization  of  the  urban  radio  multipath  channel  is 
being  developed  which  is  experimentally  exact  for  use  in  the 
design  of  systems  with  bandwidths  up  to  50  KHz  and  which,  with 
appropriate  care,  will  be  useable  out  to  150  KHz.  This  should 
be  adequate  for  the  design  and  evaluation  of  systems  developed 
for  the  present  nominal  bandwidth  allocations  and  also  allow 
sufficiently  broad-band  channel  characterization  to  study  some 
channel  multiplexing  strategies.  The  final  choice  of  channel 
characterization  bandwidth  was  based,  of  course,  on  the  need  to 
provide  reasonable  system  design  flexibility  at  a moderate  cost. 
The  equipment  for  measurement  and  hardware  simulation  of  the 
urban  multipath  channel  is  being  designed  and  built  for  DOT/TSC 
by  Signatron,  Inc.. 
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Figure  1.  Communication  System  Simulator 


The  multipath  measurement  technique  being  used  is  a "pulse- 
probing" technique  which  determines  the  band-limited  impulse 
response  of  the  multipath  channel.  Previous  measurements  have 
indicated  that  the  time  duration  of  the  channel  impulse  response 
determined  by  the  multipath  delay  spread  between  the  shortest 
and  the  longest  propagation  paths  in  the  urban  environment  is  on 
the  order  of  40  ysec.  A probing  signal  with  an  impulse-like 
time  auto-correlation  function  is  transmitted  over  the  urban 
radio  channel  and  the  channel  output  is  measured  by  correlation 
techniques,  sampled,  and  recorded.  A conceptual  picture  of  this 
process  is  shown  in  Figure  2a. 

To  simulate  the  multipath  channel,  a filter  with  the  re- 
corded channel  impulse  response  is  generated  by  using  a tapped 
delay  line.  The  effects  of  the  multipath  channel  on  an  arbi- 
trary signalling  waveform  can  be  duplicated  in  the  laboratory 
by  applying  this  same  signal  to  the  input  of  this  tapped  delay 
line  channel  model  and  varying  its  tap  gains  according  to  the 
prerecorded  channel  data.  This  process  is  shown  in  Figure  2b. 

The  mathematical  theory  of  dispersive  channel  characterization 
has  been  extensively  studied  and  is  described  in  References  1, 

2,  and  3.  Also,  a survey  of  the  variety  of  measurement  techni- 
ques for  these  channels  can  be  found  in  Reference  4. 

The  actual  channel  measurement  transmitter  hardware  uses  a 
511  bit  pseudo-random  sequence  of  phase-modulated  5 ysec  long 
raised  cosine  pulses  as  its  sharp  auto-correlation  function 
probing  signal.  The  channel  measurement  receiver  uses  a bank  of 
10  pulse  compression  filters  to  continuously  sample  the  esti- 
mated channel  impulse  response.  These  samples  are  then  recorded 
for  use  as  the  10  channel  delay  line  model  tap  gains.  Each  tap 
gain  is  a low  pass  complex  random  process  with  its  bandwidth 
proportional  to  the  rate  of  change  of  the  multipath  environment, 
or,  equivalently,  to  the  doppler  spread  of  the  urban  radio 
channel.  Measurement  data  taken  by  others  has  indicated  that 
doppler  spreads  of  under  a hundred  cycles  will  be  experienced  in 
moving  ground  vehicles. 

NOISE  MEASUREMENTS 

Urban  radio  channels  typically  exhibit  high-level  impulsive 
man-made  noise  as  well  as  low-level  background  noise.  The  noise 
environment  will  be  recorded  over  a wide  bandwidth  (up  to  1 MHz) 
and  with  extremely  large  dynamic  range  (85  dB) . While  only  a 
relatively  narrow  band,  lower  dynamic  range  characterization  of 
the  additive  noise  (approximately  100  KHz)  is  needed  for  the 
channel  simulator,  the  more  detailed  noise  structure  is  neces- 
sary for  analysis  and  modelling  efforts.  To  support  these  latter 
efforts,  a PDP-11  mini-computer  along  with  a hardware  interface 
which  is  presently  being  built  will  be  used  for  real-time  control 
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(Simplified) . 
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Figure  2b.  Channel  Measurement  and  Playback 
Technique  (Simplified) . 


of  the  actual  noise  measurement.  First-order  statistics  as  well 
as  tests  of  homogeneity  will  be  computed  for  real-time  evalua- 
tion of  the  validity,  quality,  and  type  of  measurement  data 
being  taken;  histograms  will  be  computed  as  well  as  conventional 
moments  such  as  rms  noise  power  which  alone  are  inadequate. 
Various  bookkeeping  functions  will  also  be  provided  by  the  com- 
puter such  as  logging  of  the  measurements  and  any  other  parts  of 
the  laborious  measurement  task  which  can  be  easily  automated. 

An  attempt  will  also  be  made  in  the  near  future  to  automate 
parts  of  the  multipath  measurement. 

Further  analysis  of  the  measured  noise  will  be  performed 
back  in  the  laboratory.  These  include  a second-order  analysis 
and  further  tests  of  homogeneity  of  the  measured  data.  The 
details  of  this  second  level  of  analysis  remain  yet  to  be  clear- 
ly defined.  However,  preliminary  studies  indicate  that  exten- 
sive use  can  be  made  of  the  homogeneity  test  to  determine 
stationarity , classification  and  separation  of  noise  phenomena, 
intervals  of  quasi-s tationarity , and  to  aid  in  the  estimation  of 
second-order  statistics  of  the  noise  process.  Also,  some  pre- 
liminary work  is  being  done  to  provide  foundations  for  modelling 
both  multipath  and  additive  noise  for  the  urban  channel. 

Several  candidate  techniques  for  taking  advantage  of  the 
noise  structure  have  been  proposed.  These  include  several 
multiplexing  strategies,  burst  error-correcting  codes,  noise 
canceller,  and  two-way  data  link  error  information  feedback. 

The  noise  measurements  will  provide  the  tools  for  the  evaluation 
of  these  concepts. 
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SURVEY  OF  MULTIPATH  AND  NOISE  MEASUREMENTS  IN  CITIES 


ABSTRACT 

This  paper  summarizes  the  available  information  on  the 
statistical  characteristics  of  the  urban  communication  channel. 
This  includes  the  results  of  measurements  of  the  propagation 
properties  of  urban  channels  as  well  as  the  results  of  measure- 
ments of  man-made  and  natural  electromagnetic  noise.  Statisti- 
cal models  of  the  urban  multipath  channel  based  on  these  data 
are  discussed.  Some  models  of  the  various  noise  processes  are 
described.  Finally,  the  experimental  program  and  supporting 
studies  of  the  Transportation  Systems  Center  of  the  U.S. 
Department  of  Transportation  are  presented. 
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INTRODUCTION 


An  important  aspect  of  any  controlled  transportation  system 
is  the  communication  link.  Over  this  link,  telemetrical , com- 
mand and  control  signals  are  sent  to  the  controlled  vehicles, 
and  information  on  vehicular  position,  passenger  load,  reports 
of  difficulties,  etc.,  are  sent  to  a central  command  and  control 
station.  Because  of  the  requirements  of  automatic  vehicular 
monitoring  and  control  and  the  need  for  making  more  effective 
use  of  the  limited  number  of  available  radio  channels  the  trend 
in  urban  communication  is  from  voice  towards  digital  communica- 
tion. Reliable  digital  data  transmission  requires  more  complete 
knowledge  of  channel  characteristics  than  does  the  highly  re- 
dundant voice  transmission.  For  this  reason,  the  measurement 
and  characterization  of  the  urban  communication  channel  have 
received  an  increasing  amount  of  attention  in  recent  years. 

In  most  urban  communication  systems,  information  is  trans- 
mitted to  various  terminals  over  links  in  the  VHF  band  (30-60 
MHz  and  150  MHz)  and  UHF  band  (450-470  MHz) . The  transmitted 
signals  are  subjected  to  interference  consisting  of  two  com- 
ponents: an  additive  noise  component  and  a time-variant  dis- 

persive fading  component.  The  latter  is  caused  by  vehicle 
motion  through  a continually  varying  multipath  environment.  The 
additive  noise  in  the  urban  communication  channel  consists  of 
electromagnetic  disturbances  generated  by  other  communication 
equipments,  power  lines,  electrical  equipments,  automotive 
ignition  systems,  etc.,  as  well  as  Gaussian  noise  contributed 
by  the  receiver  itself.  It  is  known  that  the  noise  generated 
unintentionally  by  the  various  man-made  sources  is  of  an  im- 
pulsive type  with  characteristics  similar  to  that  of  the  atmos- 
pheric noise  which  affects  the  performance  of  communication 
systems  at  lower  frequencies. 

This  report  summarizes  the  available  information  on  the 
statistical  characteristics  of  the  urban  communication  channel. 

The  techniques  used  and  results  obtained  in  measurements  of  pro- 
pagation properties  of  urban  channels  are  described,  together 
with  available  data  on  multipath  delay  and  Doppler  spread  of 
the  channel.  Several  statistical  models  of  the  urban  multipath 
channel  based  on  these  experimental  data  are  discussed.  The 
statistical  characteristics  of  impulse  noise  in  the  urban 
channel  are  discussed.  Models  of  the  noise  process  are  described. 

The  information  summarized  here  concentrates  on  character- 
istics of  the  urban  channel  which  are  relevant  to  the  design  of 
reliable  communication  systems  in  the  urban  environment.  It 
will  be  evident  from  this  summary  that  data  obtained  in  measure- 
ments made  to  date  on  the  multipath  channel  and  the  additive 
noise  process  in  the  urban  environment  are  inadequate  for  the 
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complete  characterization  of  the  urban  channel  and  are  often 
misleading  for  communication  system  design.  We  shall  point  out 
in  this  survey  areas  in  which  useful  work  could  be  done  to  im- 
prove our  methods  of  reliable  information  transmission  in  the 
urban  environment.  Measurements  of  the  multipath  urban  channel 
and  the  additive  man-made  noise  characteristics  which  will  be 
made  by  the  Transportation  Systems  Center  of  DOT  will  be  de- 
scribed. On  the  basis  of  these  measurements,  working  models  of 
the  noise  process  will  be  constructed  and  a channel  simulator 
will  be  built  capable  of  presenting  to  a candidate  communication 
system  the  distortion  characteristics  of  the  urban  channel. 

PROPAGATION  PROPERTIES  OF  THE  URBAN  CHANNEL 

The  problem  of  analyzing  the  propagation  properties  of  the 
urban  environment  and  the  urban  channel  characterization  has  been 
approached  by  various  investigators  on  both  theoretical  and  ex- 
perimental bases.  Early  work  on  urban  radio  propagation  has,  for 
the  most  part,  gathered  data  on  fading  distributions  and  Doppler 
spreads  at  particular  frequencies  (1-7)  . These  empirical  data 
obtained  by  transmitting  spectrally  pure  or  narrowband  CW  sig- 
nals, are  useful  in  predicting  the  performance  of  narrowband 
systems  used  in  voice  communication  links.  To  predict  the  be- 
haviour of  digital  communication  systems,  measurements  of  the 
multipath  properties  of  urban  channels  have  been  conducted  (8- 
12) . Experimental  data  on  multipath  propagation  of  pulse  sig- 
nals have  led  to  several  models  of  the  multipath  channel.  In 
this  section,  we  shall  discuss  the  results  of  the  measurements 
made  using  CW  and  pulse  techniques  and  the  associated  theoreti- 
cal models. 

MEASUREMENTS  OF  FADING  STATISTICS  USING  CW  TECHNIQUES 

In  a typical  mobile-radio  system,  one  station  is  at  a fixed 
location,  while  the  other  is  mobile.  The  received  signal  is  the 
time-varying  super-position  of  the  direct  wave  and  reflected  and 
scattered  waves.  When  the  transmitted  signal  is  narrowband, 
the  multipath  effects  of  the  propagation  medium  manifest  them- 
selves as  fluctuations  in  signal  strength  which  depend  on  the 
relative  strengths  and  phases  of  the  component  waves.  The  amp- 
litude distribution  of  the  received  signal,  therefore,  depends 
on  the  characteristics  of  the  propagation  path,  while  the  rate 
of  fading  is  a function  also  of  the  speed  of  the  mobile  station. 

Signal  Amplitude  Distribution:  It  has  been  generally 

observed  that  the  amplitude  distribution  of  the  received 
signal  is  independent  of  the  carrier  frequency  of  the  trans- 
mitted signal  over  a broad  frequency  range  (6) . By  trans- 
mitting narrowband  signals  (of  approximately  10  KHz  bandwidth) 
at  150,  450  and  900  MHz,  Young  (13)  has  found  that 
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the  amplitudes  of  the  received  signals  over  small  areas  (300  ft 
to  1000  ft  in  extent)  in  New  York  City  are  Rayleigh  distributed 
when  measured  with  a mobile  receiver.  The  mean  values  of  the 
received  signal  amplitudes  decrease  with  the  distance  between 
the  sample  area  and  the  transmitter.  Signal  amplitudes  were  also 
found  to  be  Rayleigh  distributed  at  several  other  frequencies 
(14,  15,  7)  on  suburban  streets  and  in  rural  areas  when  the 
receiver  was  far  from  the  transmitter  (at  a distance  greater  than 
4 km) . The  measured  signal  amplitude  being  Rayleigh  distributed 
suggests  that  the  received  signals  in  these  experiments  are 
composed  of  scattered  waves  of  low  strengths.  For  short  trans- 
mission distance  in  cities  and  for  transmission  over  open  areas, 
the  received  signal  amplitude  was  observed  to  be  Rice  distributed 
(and  lognormally  distributed  (6)).  This  behaviour  of  the  re- 
ceived signal  amplitude  distribution  suggests  the  presence  of 
relatively  strong  direct  signal  component. 


Doppler  Power  Spectra:  By  transmitting  an  838  MHz  CW 

signal  from  a fixed  based  station  and  recording  the  fading 
waveforms  of  the  received  signal  on  a vehicle  moving  at  15  mph 
on  city  streets  (New  Providence,  N.J.),  Ossanna  (3)  calculated 
the  Doppler  power  spectrum  of  the  received  signal.  The  power 
spectra  obtained  on  different  streets  differ  somewhat  in  details. 
In  general,  a sharp  cut-off  in  the  baseband  spectrum  at  twice 
the  maximum  Doppler  shift  frequency  was  observed  in  all 
measured  spectra. 


A representative  spectrum  in  shown  in  Fig.  1.  The  Doppler 
power  spectrum  was  also  measured  at  910  MHz  in  a suburban  en- 
vironment (1) . By  plotting  the  received  signal  amplitudes  as  a 
function  of  Doppler  shift,  it  is  possible  to  detect  distinct 
propagation  paths  with  different  Doppler  shifts. 

Frequency  Correlation  of  Signals:  To  determine  the  maximum 

usable  bandwidth  in  mobile  radio,  it  is  necessary  to  know  the 
correlation  between  signals  at  different  frequencies.  Since  the 
correlation  bandwidth  of  a channel  is  inversely  related  to  its 
multipath  spread,  by  measuring  the  covariance  of  two  signals  as 
a function  of  their  frequency  separation,  one  can  also  obtain 
the  distribution  of  time  delays  experienced  by  the  component 
waves  of  the  received  signal.  The  only  known  empirical  data  on 
the  covariance  of  signals  were  obtained  by  Ossanna  in  a measure- 
ment at  860  MHz  in  a suburban  environment  for  carrier  frequency 
separations  of  0.1,  0.5,  1.0  and  2.0  MHz  (13).  The  experimental 
data  points  fit  roughly  a theoretical  curve  corresponding  to  an 
exponentially  distributed  multipath  delay  with  a mean  of  1/4  y 
sec . 
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Power  spectrum  in  dB 


0 


Frequency  in  cycles/second 


-■■  ■ Theoretical  spectrum 

Measured  spectrum 

Speed  of  the  vehicle:  15.8  miles  per  hour. 

Angle  between  the  direction  to  the  fixed  station 
and  the  direction  of  vehicle  travel:  83.3°. 


Figure  1.  Typical  Doppler  Power  Spectrum 

(Relative  Amplitudes  are  Arbitrary) , 
Figure  10  in  (37) . 
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Spatial  Correlation;  Using  a space  diversity  system  to 
receive  frequency-modulated  signals , Rustako  (7)  measured  the 
probability  distribution  of  the  predetection  combined  signals 
of  one  to  four  channels  with  1/4 , 3/4,  and  5/4  wavelength 
antenna  spacing.  For  a single  channel,  the  signal  fading  was 
nearly  Rayleigh  distributed.  As  more  channels  were  combined, 
the  distribution  became  different  from  Rayleigh  with  noticeable 
reduction  in  fading  depth.  Moreover,  the  fading  was  observed  to 
decrease  with  antenna  spacing.  The  increase  in  diversity  im- 
provement with  antenna  spacing  indicates  correlation  between 
received  signal  and  coupling  between  antennas  at  these  spacings . 

Theoretical  Models  for  CW  Case:  A number  of  phenomenologi- 

cal models  of  the  urban  propagation  channels  have  been  proposed 
to  explain  the  experimentally  observed  characteristics  of  the 
received  mobile-radio  signals.  These  models  assume  that  the 
received  field  is  composed  of  a number  of  interfering  waves. 
Fluctuations  in  received  signal  amplitude  occur  because  of  the 
motion  of  the  mobile  station  through  the  spatial  interference 
pattern.  Depending  on  the  properties  of  the  interfering  waves, 
the  models  can  be  classified  into  two  types:  reflection  models 

and  scattering  models. 

The  reflection  model  was  constructed  by  Ossanna  (3)  to 
predict  theoretically  the  Doppler  power  spectra  of  the  fading 
waveforms  observed  in  measurements  in  a suburban  environment. 

The  model  assumes  that  reflections  occur  at  the  flat  sides  of 
houses,  modeled  as  randomly  placed  vertical  plane  reflectors. 
The  incident  and  reflected  waves  form  a standing  wave  pattern. 
Using  the  radio-carrier  frequency,  vehicle  speed,  and  the  angle 
between  the  direction  of  vehicle  motion  and  the  direction  to 
the  transmitter  as  parameters,  theoretical  Doppler  power  spectra 
were  computed.  In  general,  the  theoretical  spectra  are  in  good 
agreement  with  those  obtained  in  the  measurement  discussed  above 
except  at  low  frequencies  and  at  frequencies  around  the  maximum 
Doppler  shift  frequency.  The  discrepancies  are  due  to  the  fact 
that  the  theoretical  model  failed  to  consider  the  effects  of 
shadowing  by  buildings,  ground  reflection,  and  non-randomness  of 
the  reflectors. 

Somewhat  better  match  of  the  gross  behaviours  of  the  theoretical 
and  experimental  spectra  were  obtained  by  Clarke  (4)  using  a 
scattering  model.  The  model  assumes  that  the  received  field  is 
made  up  of  a number  of  horizontally  traveling  plane  waves,  in- 
cluding scattered  waves,  whose  azimuthal  angles  of  arrival  and 
phases  are  statistically  independent.  The  phases  of  the  com- 
ponent scattered  waves  are  evenly  distributed.  The  probability 
that  a component  wave  will  occur  in  the  azimuthal  angular  sector 
(a,  at  da)  and  the  probability  density  functions  of  the  time- 
delays  of  the  component  waves  can  be  chosen  for  different  en- 
vironments. The  received  field  may  also  contain  a direct  wave 
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component  whose  phase  and  strength  depend  on  the  distance  and 
attenuation  of  the  direct  propagation  path  between  the  receiver 
and  transmitter. 

With  this  model,  it  is  possible  to  show  that  when  the 
strength  of  the  direct  wave  is  weak  compared  to  the  scattered 
waves,  the  envelope  of  the  total  received  field  is  Rayleigh 
distributed  and  its  phase  is  evenly  distributed,  as  observed  in 
experiments.  When  there  is  a strong  direct  wave  present,  the 
resulting  envelope  and  phase  will  no  longer  be,  respectively, 
Rayleigh  and  uniformly  distributed.  The  phase  and  amplitude 
distributions  are  those  derived  by  Rice  for  a sine  wave  plus 
random  noise.  Clarke  was  also  able  to  derive  the  spatial  auto- 
correlation and  covariance  functions  of  the  electric  and  magnetic 
fields  for  the  particular  case  in  which  the  arrival  angles  of 
the  component  waves  are  uniformly  distributed.  For  antennas 
spaced  at  1/4,  3/4,  and  5/4  wavelength  apart,  the  spatial  cor- 
relation coefficients  of  the  signal  amplitudes  were  found  to  be 
0.25,  0.06,  and  0.03,  respectively.  The  correlation  of  signals 
at  different  frequency  separations  may  be  derived  from  the  pro- 
bability density  function  of  the  time-delays  between  the  com- 
ponent waves . 

Gilbert  (16)  examined  several  scattering  models  in  which 
the  angles  of  arrival  of  the  component  waves  are  uniformly  dis- 
tributed. It  was  shown  that  when  the  number  of  component  waves 
is  sufficiently  large,  all  models  in  which  the  phases  of  the 
component  waves  are  evenly  distributed  are  equivalent,  as 
expected. 

Pulse  Measurements  of  Studies  of  Young  and  Lacy:  To 

determine  the  maximum  bandwidth  or  the  shortest  time  interval 
that  may  be  used  satisfactorily  in  some  wideband  multi-channel 
voice  communication  systems,  a measurement  of  the  multipath 
characteristics  of  the  urban  channel  was  conducted  at  450  MHz 
in  New  York  City  by  Young  and  Lacy  (19).  In  the  experiment, 
pulses  of  1/2  ysec  duration  were  transmitted  from  the  base 
station.  The  channel  response  was  detected  by  the  receiver  on 
the  mobile  unit.  No  information  on  the  absolute  time-of-arrival 
was  obtained.  Multipaths  whose  path  losses  are  between  6 and 
12  dB  more  than  that  of  the  main  path  (the  path  corresponding 
to  the  strongest  component  in  the  received  signal)  may  occur  in 
a delay  spread  period  of  15  to  20  ysec.  It  was  observed  that 
multipaths  appear  and  disappear  and  their  strengths  change 
rapidly  as  the  mobile  unit  moves  at  speeds  10  to  15  mph  through 
the  city.  Unfortunately,  no  quantitative  data  were  obtained  in 
the  experiment. 

Schmid  Propagation  Model:  Based  on  the  results  of  Young 

and  Lacy,  Schmid  (8)  constructed  a statistical  propagation  model 
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to  estimate  the  probability  of  occurrence  of  multipath.  The 
model  predicts  that  for  a constant  distance  between  the  trans- 
mitter and  the  receiver,  the  multipath  strength  decreases  with 
increasing  multipath  delay.  And  for  a given  delay,  the  probabil- 
ity of  occurrence  of  a multipath  increases  with  the  distance 
between  the  transmitter  and  the  receiver.  That  is,  the  width  of 
the  probability  of  multipath  occurrence  as  a function  of  excess 
path  delay  increases  with  the  distance.  (This  has  not  been 
supported  by  the  experiment  results  obtained  in  the  Teknekron 
experiment  described  next.)  Although  the  Schmid  model  offers 
some  rough  indication  of  the  multipath  properties  of  the  channel, 
it  is  of  little  value  for  communication  system  design. 

Measurement  and  Theoretical  Models  of  Urban  Channel  by 
Teknekron:  Recently , for  the  analysis  and  design  of  a wideband 

radio-location  system,  Teknekron  (9,  17)  conducted  more  exten- 
sive measurements  of  the  multipath  porperties  of  the  urban  chan- 
nel. The  measurements  were  designed  to  determine  the  parameters 
of  a statistical  model  which  considers  the  urban  channel  as  a 
time-invariant  linear  filter.  The  impulse  response  of  the  filter 
was  assumed  to  be  of  the  form 


y (t)  = ^ ak  e1^  6 (t-tk) 
k=o 


where  tQ  is  the  delay  of  the  direct  signal.  - tQ,  k = 1,  2, 

represent  the  excess  delays  of  the  multipath  signals.  The  sets 
of  path  variable  {a^}  and  (tk)  are  statistically  independent  of 
the  carrier  phases  4>k  of  the  paths  were  assumed  to  be  statisti- 
independent,  evenly  distributed  random  variables.  Studies  of 
the  statistics  of  the  path  strengths  a^  in  measurements  using 
narrowband  signals  4,  6,  8,  2 indicate  that  over  a small,  sample 
reception  area  the  path  strengths  a^  have  Rayleigh  or  Rice 
distributions,  while  for  larger  areas  the  strengths  become  log- 
normally  distributed.  The  purpose  of  the  experiment  was  to 
determine  the  statistics  of  the  path  variables  a^  and  t^  and 
their  statistical  dependence. 

In  the  experiment,  100  nanosecond  pulses  at  each  of  the 
three  frequencies  488  MHz,  1280  MHz,  and  2920  MHz  were  trans- 
mitted simultaneously  from  fixed  transmitters.  The  responses  of 
the  channel  at  all  three  frequencies  were  received  by  the  mobile 
unit  (not  in  motion)  with  the  video-detected  response  displayed 
on  an  oscilloscope.  The  scope  was  triggered  t seconds  after  the 
transmitters  were  pulsed,  where  x was  set  to  be  slightly  less 
than  the  smallest  direct  path  delay  expected  in  the  reception 
area.  The  scope  traces  were  photographed  and  then  reduced  on 
an  optical  scanning  - table.  Each  of  the  peaks  and  inflection 
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points  in  the  video  response  trace  was  taken  to  be  a distinct 
multipath.  Its  occurrence  time,  measured  with  respect  to  the 
occurrence  time  of  the  first  peak  or  inflection  point  in  the 
trace,  is  the  excess  delay  of  the  path.  Its  amplitude  was  taken 
to  be  the  path  strength. 

Data  obtained  in  San  Francisco  in  areas  of  various  multipath 
propagation  characteristics  indicated  that,  with  minor  excep- 
tions, the  multipath  properties  of  the  urban  environment  are 
similar  for  all  three  frequencies.  Figure  2 shows  the  probabil- 
ity of  occurrence  of  one  or  more  multipaths  within  +50  feet  of 
the  excess  range  for  various  delay  ranges  at  1280  MHz.  This 
result  shows  that  the  multipath  spread  depends  almost  totally  on 
the  local  environment  of  the  receiver  contrary  to  that  predicted 
by  Schmid's  theoretical  model. 

The  empirical  data  on  excess  delays  t^  - tQ ; k = 1 , 2 , 3 
of  the  multipath  signals  indicates  that  the  sequence  t^;  k = 1, 

2,  roughly  form  a Poisson  process.  Multipaths  with  small  excess 
delays  tend  to  be  more  clustered  than  a Poisson  sequence  would 
dictate.  Better  fit  between  the  empirical  distribution  and 
Poisson  distribution  was  observed  for  multipaths  of  larger  ex- 
cess delays. 

The  multipath  strengths  a^  were  found  to  be  lognormally 
distributed  with  means  and  standard  deviations  dependent  on  the 
corresponding  multipath  delay  t^.  Figure  3 shows  the  dependence 
of  the  mean  of  a^  as  a function  of  the  excess  range  for  different 
types  of  environments.  The  strengths  of  different  paths  were 
found  to  be  correlated.  The  correlation  coefficients  of  the 
path  strengths  are  yet  to  be  determined. 
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Figure  2.  Probability  of  Given  Multipath  Occurence 
for  Various  Range.  (Fig.  4 in  [9]). 
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Figure  3.  Dependence  of  Average  Multipath  Strength 
on  Excess  Range  (Fig.  7 in  [9]). 


NOISE  IN  THE  URBAN  COMMUNICATION  CHANNEL 


The  noise  external  to  the  receiving  antenna  is  composed  of 
emissions  from  natural  sources  and  electromagnetic  pollution  due 
to  man-made  sources.  Based  on  the  average  noise  power  levels, 
the  nautral  noise  processes  contirbute  the  dominant  noise  at  the 
extremes  of  the  radio  spectrum  in  urban  environment.  In  addi- 
tion to  the  thermal  radiation  from  the  local  environment,  natural 
noise  consists  of  atmospheric  radio  noise,  galactic  radio  noise, 
and  sky  noise.  Most  communication  systems  in  the  urban  environ- 
ment operate  in  the  VHF  and  UHF  bands.  In  these  systems,  the 
effects  of  man-made  electromagnetic  pollution  are  automotive 
ignition  systems,  power  lines  and  electric  equipments.  Clearly, 
the  contribution  of  each  class  of  sources  depends  on  the  radia- 
tion spectrum  and  the  number  of  sources  of  that  class  present  in 
the  environment.  In  urban  and  suburban  environments,  at  fre- 
quencies below  20  MHx,  the  dominant  man-made  sources  seem  to  be 
power  lines,  and  electric  equipments  (such  as  r-f  stabilized  arc 
welders,  brush-type  electric  motors,  large  electric  contacts, 
etc.).  At  frequencies  above  20  MHz,  the  most  important  noise 
source  is  automotive  ignition  systems. 

The  characteristics  of  noise  generated  by  various  types  of 
noise  sources  have  been  studied  extensively  (18-26) . The  rela- 
tive field  strengths,  statistical  characteristics,  frequency  and 
distance  dependences  of  these  noise  fields  have  been  reported 
comprehensively  in  several  review  papers  (26-29),  there  is  no 
need  to  describe  them  in  detail  here.  Generally  speaking,  in- 
cidental man-made  noise  is  impulsive  in  waveform  and  aperiodic 
in  occurrence.  When  recorded  with  wideband  receivers,  the  wave- 
forms of  these  noises  were  observed  to  be  bursts  of  narrow 
pulses.  (For  the  power  line  noise,  the  durations  of  the  pulses 
were  observed  to  be  in  the  tens  of  nanoseconds  with  the  bursts 
lasting  several  milliseconds.  Automotive  ignition  noise  occurs 
as  nanosecond  pulses  in  bursts  lasting  for  tens  of  microseconds 
to  several  milliseconds.)  Since  the  pulses  have  very  narrow 
waveforms,  the  amplitude  of  the  very  high  frequencies  of  the 
spectra  remains  large  even  though  the  rates  of  occurrence  of  the 
noise  pulses  are  low.  Figure  4 shows  the  median  operating  noise 
power,  Fam , in  dB  above  kTQB  (k  is  Boltzman's  constant,  TQ  is 
288.37  degrees  Kelvin  and  B is  the  effective  noise  bandwidth  of 
the  receiving  system  in  Hz)  from  man-made  noise  sources  together 
with  that  from  natural  noise  sources  and  these  values  are  the 
levels  expected  from  an  omnidirectiona , short,  lossless,  vertical 
antenna  near  the  surface  (29) . The  composite  frequency  spectra 
of  major  man-made  sources  are  shown  in  Figure  5.  Figures  6 and 
7 show  the  typical  values  of  the  average  occurrence  rates  of 
pulses  for  power  line  noise  and  automotive  ignition  noise, 
respectively  (30) . Figure  8 shows  the  typical  behaviour  of  the 
amplitude  probability  distribution  of  man-made  noise  when  ob- 
served over  a narrow  band  (31) . 
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Figure  4 . Median  Values  of  Radio  Noise  Power  (Omnidirectional 
Antenna  Near  Surface) . 
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Figure  5.  Composite  Frequency  Spectra  of  Noises  from  Major 
Man-Made  Sources.  (Fig.  2 in  [26]). 
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Figure  6.  Power  Line  Noise  in  Clear  Weather.  (Graph 
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Figure  7.  Automotive  Ignition  System  Noise.  (Graph 
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Figure  8.  Typical  Amplitude  Probability 

Distribution  of  Man-Made  Noise  [3]. 
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The  amount  of  published  information  concerning  measurements 
of  man-made  noise  in  the  urban  environment  is  truly  remarkable. 
Due  to  the  difference  in  the  scope  of  the  measurements,  in  the 
techniques  used,  and  in  the  detail  with  which  they  have  been 
reported,  the  comparison  and  characterization  of  these  results  is 
often  an  impossible  task.  The  report  by  the  Joint  Technical 
Advisory  committee  of  the  IEEE,  Subcommittee  63.1.3  on  unintended 
radiation  (29) , represents  a good  attempt  to  summarize  these 
results.  It  contains  most  of  the  relevant  information  on  the 
measurements  made  on  man-made  noise  prior  to  1968.  The  report 
also  contains  some  conclusions  on  the  completeness  of  the  avail- 
able information  on  man-made  noise  and  some  recommendations  for 
future  work.  Not  included  in  the  report  are  some  measurements 
made  in  Japan  at  VHF  and  UHF  with  very  wideband  receivers  (32, 
33).  Few  additional  measurements  have  been  reported  after  1968. 
Of  special  interest  are  the  ones  made  by  the  Radiotelevisione 
Italiana  at  850  MHz  and  12  GHz  (34)  and  the  Institute  of  Tele- 
communication Science  (31) . 

There  is  little  doubt  that  man-made  urban  noise  is  a pro- 
blem worthy  of  serious  consideration  both  for  analog  and  digital 
communication.  Especially  in  digital  communications,  even  a 
modest  level  of  impulse  noise  can  seriously  degrade  the  perform- 
ance of  system  not  designed  to  withstand  this  type  of  inter- 
ference. (For  example,  it  was  reported  (35)  that  an  essentially 
error-free,  low  data  rate  (2400  bps)  digital  system  whose  per- 
formance was  degraded  to  an  error  rate  of  10~2  during  the  inter- 
val of  time  corresponding  to  the  passing  by  of  a single  automo- 
bile. The  theory  presented  in  (36)  makes  it  clear  that  the  pro- 
blem is  especially  acute  for  low  data  rate  systems) . Unfortun- 
ately, there  are  very  few  data  and  validated  models  available 
to  adequately  describe  the  wideband  characteristics  of  man-made 
noise  for  the  purpose  of  reliable  communication  systems  design. 
Measurements  of  man-made  noise  have  been  made  almost  solely  with 
narrow  band  receivers . Data  obtained  in  the  measurements  are 
not  sufficient  for  determining  the  effects  of  the  noise  on  the 
performance  of  wideband  systems.  The  statistical  parameters 
measured  have  been  very  simple  ones  such  as  first  and  second 
moments,  first  order  probability  density  functions  of  the  de- 
tected noise  envelope.  Most  of  the  results  are  presented  either 
as  average  noise  power  levels  with  averaging  time  on  the  order 
of  seconds  or  minutes,  or  as  percentage  of  the  total  number  of 
pulses  which  exceeds  a given  level.  Examples  of  these  results 
are  given  in  Figures  4-8. 

Very  little  work  has  been  done  to  date  on  the  problem  of 
modeling  the  man-made  noise  process.  The  only  man-made  noise 
models  are  several  simple  one-  and  two-dimensional  ones  proposed 
by  Ottesen  (37).  In  these  models,  the  field  amplitude  at  any 
receiver  site  is  the  random  phasor  sum  of  the  field  generated 
by  sources  scattered  randomly  in  space  and  overlapping  in  time. 
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By  adjusting  the  parameters  of  the  theoretical  model,  it  is 
possible  to  match  the  theoretical  amplitude  probability  distri- 
butions of  the  noise  envelope  to  those  obtained  experimentally. 
However,  these  models  are  inadequate  for  characterization  of  the 
noise  process  in  the  problem  of  reliable  communication  system 
design. 

Scientists  at  ITS  have  proposed  to  use  some  of  the  models 
for  atmospheric  noise  to  model  man-made  noise  (31) . The  latter, 
being  of  impulsive  nature,  behaves  very  similarly  to  atmospheric 
noise  at  least  when  observed  over  a narrow  band.  The  atmospheric 
noise  was  modeled  by  Furutsu  and  Ishida  (38)  as  a Poisson  - 
Poisson  process  in  the  sense  that  bursts  of  pulses  occur  at  ran- 
dom and  the  times  of  arrival  of  the  pulses  in  the  bursts  are  also 
Poisson  distributed.  An  explicit  method  is  given  in  (38)  to 
evaluate  theoretical  amplitude  probability  density  functions. 
However,  to  apply  the  model  to  urban  man-made  noise,  more  ana- 
lysis is  needed  to  investigate  whether  the  conditions  prevailing 
in  the  urban  environment  allow  analytical  simplifications  of  the 
probability  density  function.  The  Hall  model  (39)  is  the  pro- 
duct of  two  independent  zero  mean  random  processes.  One  of 
these  processes  is  Gaussian.  A refinement  of  this  model  for 
modeling  the  nam-made  noise,  which  will  be  discussed  in  this 
conference  by  Middleton  (40) , is  part  of  the  program  to  be 
carried  out  at  the  Transportation  Systems  Center. 
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PROGRAM  OF  THE  TRANSPORTATION  SYSTEMS  CENTER 


The  U.S.  Department  of  Transportation  at  its  Transportation 
Systems  Center  in  Cambridge,  Massachusetts  has  initiated  a pro- 
gram of  characterization  and  measurement  of  the  urban  communica- 
tion channel  with  the  long-range  goal  of  improving  the  perform- 
ance of  communication  systems  in  the  urban  environment  by 
establishing  a capability  for  comparative  evaluation  of  competing 
components,  subsystems  and  systems,  on  an  equal  basis,  taking 
into  account  the  effects  of  the  dispersive  channel  characteris- 
tics as  well  as  the  additive  noise.  In  addition,  supporting 
theoretical  studies  are  being  performed  particularly  in  the  area 
of  physical  modeling  of  noise  sources.  The  two  main  areas  of 
experimental  work  are  the  programs  of  noise  measurement  and 
channel  characterization  and  recording.  Included  in  the  latter 
is  the  fabrication  of  a channel  playback  simulator  which  will 
allow  the  comparative  evaluations,  previously  reffered  to,  to  be 
carried  out.  We  will  discuss  these  two  areas  of  experimental 
work  with  emphasis  on  the  specification  of  the  equipment  and  the 
design  of  the  experiment,  as  far  as  present  program  status  allows. 
Some  comments  about  impulsive  noise  models  will  also  be  made. 

Noise  Measurement  Program:  A communication  system  planner  would 

like  to  know  the  detailed  characteristics  of  the  background 
noise  in  order  to  evaluate  its  effects  on  system  performance. 

This  information  includes  the  time  and  space  variation  of  the 
noise  as  well  as,  for  prediction  purposes,  its  correlation  with 
factors  such  as  population  density,  vehicular  traffic,  power 
consumption  and  topographical  characteristics.  As  has  been 
mentioned  the  requirements  are  more  stringent  for  digital  com- 
munication systems  which  require  a more  thorough  specification 
of  the  statistical  characteristics  of  the  noise.  By  and  large 
this  information  is  not  currently  available  and  this  lack  of 
design  data  hampers  designers  of  wideband  systems,  i.e.,  systems 
whose  bandwidth  exceeds  10  kHz.  Therefore,  the  purpose  of  the 
experiment  under  consideration  is  to  measure  broadband,  wide 
dynamic  range,  impulsive-type  noise  in  selected  regions  of  the 
VHF  and  UHF  band.  The  bandwidth  of  the  measuring  system  has 
been  tentatively  selected  to  be  1 MHz.  The  predominant  noise  is 
expected  to  be  generated  by  man-made  machines  such  as  automobile 
engines,  industrial  equipment,  etc. 

The  reasons  for  this  experimental  program  are  the  following: 

1.  To  study  experimental  techniques  for  the  measurement  of 
wideband  man-made  noise  in  the  VHF-UHF  region  and 
establish  a capability  for  the  routine  performance  of 
such  measurements. 
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2.  To  identify  the  major  sources  of  radiation  that  will 
interfere  with  the  performance  of  any  transportation 
system  being  planned  in  the  frequency  region  mentioned 
previously . 

3.  To  determine  the  information  needed  for  the  design 
and  evaluation  of  any  communication  or  command  and 
control  system  of  mass  transit  or  land,  air,  water  or 
highway  communication.  In  combination  with  channel 
probing  measurements  this  information  will  provide  a 
complete  characterization  of  the  urban  and  suburban 
channel . 

4.  To  determine  the  radiation  characteristics  of  existing 
and  proposed  transportation  systems  such  as  electric 
cars,  trolley  cars,  electric  trains,  etc.  Figure  9 
gives  an  overview  of  the  entire  program  showing  the 
relationships  between  measurements,  data  processing 
and  the  expected  outputs  of  the  program. 

Measurement  Parameters:  One  of  the  first  tasks  is,  of  course, 

the  determination  of  the  statistical  quantities  describing  the 
noise  which  can  be  derived  by  suitable  processing  of  the  ex- 
perimental data.  Although  a certain  amount  of  data  processing 
can  be  done  in  real  time  the  decision  has  been  made  to  record 
the  IF  waveform  for  later  processing  and  for  use  in  the  simu- 
lation equipment.  This  means  that  most  of  the  processing  pro- 
cedures can  be  decided  upon  at  a later  date.  Some  initial  con- 
sideration on  what  has  to  be  measured  is,  however,  in  order. 

A simple  model  of  the  wideband  man-made  interference  is  the  sum 
of  an  impulsive  component  and  of  a Gaussian  process.  Both  com- 
ponents can  be  expected  to  be  nonstationary.  It  seems  neces- 
sary to  measure  at  least  the  following  quantities. 

1.  Average  power  over  each  measurement  interval 

2.  Amplitude  pdf  of  the  peaks 

3.  Duration  pdf  of  the  impulses  at  given  amplitude  levels 

4.  Pulse  spacing  pdf 

5.  Crossing  rate  pdf 

In  addition  to  the  above  quantities  some  simpler  parameters 
such  as  average  and  rms  values  after  detection  can  be  measured 
in  real  time.  It  is  worth  pointing  out  that  instrumentation 
similar  to  that  being  studied  for  the  noise  measurements  can  be 
used  in  the  future  for  the  investigation  of  channel  occupancy. 

By  this  we  mean  the  evaluation  of  the  probability  that  a given 
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Figure  9.  Flow  Chart  of  Major  Tasks,  Problem  Areas  and 

Logical  Flow  of  Data  for  the  Measurement  Program. 


frequency  region  be  occupied  at  a certain  level,  by  other  trans- 
mitters. Because  of  the  limited  spectral  region  available  for 
land  mobile  communication  it  is  also  important  to  determine  the 
minimum  allowed  spatial  separation  of  stations  transmitting  in 
the  same  frequency  region  in  order  to  avoid  appreciable  cochannel 
interference.  This  type  of  measurement  can  also  be  performed  by 
the  system  under  consideration. 

A simple  block  diagram  of  the  noise  measuring  system  is 
shown  in  Figure  10 . 

Design  of  Noise  Measurement  Experiment:  The  characterization  of 

the  man-made  noise  present  in  the  urban  environment,  and  there- 
fore the  design  of  an  appropriate  experiment  to  derive  such 
characterization,  appears  to  be,  at  first  a rather  hopeless  task 
because  of  the  multifarious  variability  of  that  noise  with  a 
number  of  variables  such  as  time,  space,  hour  of  the  day,  traffic 
conditions  etc.  In  addition,  the  difficulty  of  the  task  is  com- 
pounded by  the  fact  that  some  of  those  variations,  for  example 
spatial  variations,  are  seldom  smooth  owing  to  occasional  sharp 
discontinuities  in  the  distribution  of  noise  sources  and  to  the 
presence  of  standing  wave  noise  patterns.  Fortunately,  however, 
a very  large  number  of  similar  experiments  have  already  been 
performed  by  other  investigators  and  we  can  benefit  from  their 
findings  and  try  to  improve  our  knowledge  of  the  urban  noise  in 
the  iterative  fashion  that  is  typical  of  statistical  experimenta- 
tion . 


The  first  important  subdivision  of  the  characterization 
problem  comes  from  recognizing  that  there  are  basically  three 
types  of  terminals  in  mobile  communication: 

1.  stationary  vehicles 

2.  moving  vehicles 

3.  base  stations 

We  will  thus  break  the  design  of  our  experiment  in  three  cate- 
gories corresponding  to  the  three  types  of  terminals. 

For  stationary  vehicles  we  can  start  with  a concept  put 
forward  in  the  FCC  study  (30),  that  of  "activity  zone".  The  FCC 
Report  divides  the  total  urban  and  suburban  environment  into 
areas  characterized  essentially  by  their  level  of  man-made 
noise.  General  categories  defined  in  the  study  were:  a)  clear 

areas,  b)  residences,  c)  businesses,  d)  industrial  sites,  e) 
transportation  areas,  f)  institutional  sites,  g)  automotive 
activities.  For  the  purpose  of  locating  most  of  these  zones  in 
a large  city  like  the  Greater  Boston  area  use  can  be  made  of  an 
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Figure  10.  Basic  Receiver  Block  Diagram 


Urban  Atlas  recently  published  (41) . This  atlas  contains  de- 
tailed information  on  the  characteristics  of  20  large  U.S. 
cities  including  Boston.  Population  density  is , for  example, 
given  with  a resolution  grid  of  250  m x 250  m and  the  different 
grids  of  each  city  are  divided  in  five  categories  (50-200,  200- 
500,  500-1200,  1200-3600  and  more  people/grid).  In  the  same 
document  the  following  subdivisions  are  also  identified  for  a 
large  city:  Airport,  Institutional,  Large  Institutional,  Park, 

Industrial,  Commercial. 

Examination  of  the  FCC  measurements  and  some  additional 
thought  suggests  the  following  locations  where  measurements 
should  be  made  with  a stationary  vehicle. 

1.  Clear  areas. 

2.  Residential  - separation  between  homes  500  ft. 

3.  Residential  - separation  less  than  30  ft.  but  detached. 

4.  Residential  - row  houses. 

5.  Apartment  buildings  - 2 stories. 

6.  Apartment  buildings  - 6 stories  and  higher. 

7.  Office  buildings  - 8 stories  and  more. 

8.  Industrial  plant. 

9.  High  power  lines  (220  KV) . 

10.  Low  power  lines  (33  KV) . 

11.  Power  plant. 

12.  Hospitals. 

13.  Main  thoroughfare. 

14.  Busy  intersection. 

15.  Large  institution. 

16.  Airport. 

It  is  expected  that  these  "activity  zones"  will  be  sufficiently 
different  in  the  amount  and  type  of  man-made  noise  present,  to 
justify  this  subdivision.  The  experiment  will,  of  course,  con- 
firm or  modify  this  assumption.  However,  even  if  the  noise 
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characteristics  were  to  be  found  very  similar  for  two  or  more 
categories , the  equivalence  of  different  types  of  sites  (say  a 
power  line  and  a congested  traffic  situation)  will  certainly 
constitute  important  information.  It  is  also  expected  that  for 
each  site  there  will  be  at  least  one  "control  variable"  to  take 
into  consideration:  for  example,  the  hour  of  the  day  in  a 

residential  area,  the  traffic  at  a busy  intersection  (smooth 
versus  a traffic  jam),  weather  conditions  for  a power  line,  etc. 
Economy  of  experimentation  suggests  the  consideration  of  only 
two  levels  for  this  variable.  It  seems  also  advisable  to  repli- 
cate the  same  experiment  at  least  twice  for  each  type  of  site  in 
order  to  gain  some  insight  into  the  validity  of  the  assumption 
of  similarity.  It  must  also  be  mentioned  that,  owing  to  the 
fact  that  automotive  noise,  whenever  present,  is  an  important 
source  of  noise,  it  will  be  probably  impossible,  or  at  least 
very  difficult,  to  measure  the  noise  characteristics  of  some  of 
the  sites  mentioned  without  including  automotive  noise.  This, 
however,  must  be  considered  a fact  of  life  (office  buildings  are 
generally  located  in  areas  of  high  traffic)  and  it  is  therefore 
appropriate  and  relevant  to  measure  to  total  noise  present  at 
those  sites. 

The  total  number  of  measurements  required  will  in  general  be 

N = C X K X M 

C = number  of  activity  areas  identified 

K = number  of  levels  of  the  control  variable 

M = number  of  replicates  for  the  same  activity  area 

From  the  above  considerations  it  seems  that  a minimum  of  64 
measurements  are  necessary  for  stationary  vehicles.  The  number 
of  sites  is  practically  coincident  with  the  number  of  measure- 
ments since,  in  general,  changing  the  level  of  the  control 
variable  will  prevent  using  a site  for  more  than  one  measurement. 

The  characterization  of  the  urban  noise  for  a terminal  on 
a moving  vehicle  should,  in  principle,  require  the  duplication 
of  the  measurements  taken  at  a stationary  site.  It  seems,  how- 
ever, that  the  main  difference  between  a moving  vehicle  and  a 
stationary  one  arises  when  the  moving  vehicle  is  immersed  in 
traffic.  This  point  of  view  was  implicitly  taken  in  the  FCC 
measurements  where  the  only  measurements  from  moving  vehicles 
were  made  while  riding  in  downtown  traffic  and  on  a busy  sub- 
urban thoroughfare.  These  measurements  should  be  duplicated, 
perhaps  adding  a few  additional  ones  in  other  "activity  zones" 
to  compare  the  results  with  those  obtained  from  a stationary 
vehicle. 
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Measurements  of  urban  noise  from  base  stations  differ  from 
the  ones  from  stationary  vehicles,  essentially  because  of  the 
height  of  the  receiving  antenna.  Furthermore,  most  base  stations 
are  presently  located  either  downtown  or  in  busy  suburban  dis- 
tricts. It  is  suggested  that  the  top  of  the  TSC  building  be 
used  for  a set  of  such  measurements.  Two  or  three  additional 
measurements  with  antennas  located  at  various  heights  will 
probably  suffice  for  the  characterization  of  the  base  station 
noise.  Remarks  similar  to  the  ones  made  for  the  accuracy 
achievable  apply  here  also. 

Channel  Characterization  Experiment:  We  now  consider  the  experi- 

mental work  planned  for  study  of  the  second  component  of  inter- 
ference caused  by  the  fading  dispersive  characteristics  of  the 
urban  channel.  This  component  arises  from  vehicle  motion  through 
a continually  changing  multipath  environment.  It  is  of  equal 
importance  to  measure  these  deleterious  effects  and  to  analyti- 
cally model  them:  measurement  is  important  because  it  provides 

a realistic,  repeatable  replica  of  the  communication  channel 
which  can  be  used  in  optimizing  or  refining  a system  design,  or 
in  evaluating  contending  candidate  systems  at  different  times  on 
an  equal  basis;  analytic  work  in  this  area  is  important  in  de- 
fining the  form  of  a channel  simulator  even  when  the  problem  is 
too  complex  for  a detailed  mathematical  analysis,  and  in  pro- 
viding the  basic  design  for  suitable  transmitter  waveforms  for 
use  in  the  fading  environment. 

Based  on  extensive  research  carried  out  over  the  last 
decade  in  the  area  of  characterization  of  fading  dispersive 
channels  it  is  now  practical  to  build  a channel  simulator  capable 
of  presenting  to  a candidate  communication  system  the  simulated, 
stored  replica  of  the  distortion  characteristics  of  any  desired 
channel  displaying  multipath  delay,  Doppler  spread  and  additive 
noise.  This  channel  playback  approach  effectively  brings  the 
field  conditions  into  the  laboratory  and  makes  possible  the  com- 
parative evaluation  of  competing  techniques  under  identical 
conditions.  This  stored  channel  simulation  can  represent  an 
estimate  of  the  channel  characteristics  based  upon  mathematical 
modeling  or  can  in  fact  be  a faithful  reproduction  of  the  channel 
under  operating  conditions.  The  availability  of  both  approaches 
is  desirable  for  a comprehensive  simulator  facility. 

The  proposed  technique  involves  two  procedures  with  cor- 
responding equipment.  The  first  procedure  involves  field  tests 
employing  specially  designed  transmitter  probing  signals  which 
cover  the  frequency  bands  of  interest  in  the  planned  program 
areas.  The  probing  signals  are  processed  in  the  moving  vehicle 
to  extract  the  time-variant  system  function  of  the  propagation 
environment  or  medium  and  record  it  on  a multi-channel  tape 
recorder.  In  addition  the  additive  noise  would  be  recorded.  A 
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tape  library  could  be  built  up  which  would  adequately  cover  all 
environments  for  which  systems  are  being  developed.  The  replay 
procedure  is  carried  out  using  these  tape  records.  The  time- 
varying  channel , with  the  additive  noise,  would  be  played  back 
repeatedly  to  enable  tests  of  various  systems  to  be  made  under 
the  same  field  conditions.  Basic  block  diagrams  of  the  signal 
processing  operations  carried  out  in  the  channel  prober  and 
playback  systems  are  shown  in  Figures  11  and  12. 

Figure  11  defines  the  basic  signal  processing  operations  of 
the  probing  system.  A periodic  probing  signal  with  complex 
envelope  z (t)  is  filtered  at  the  transmitter  by  IF  and  RF  filters 
prior  to  transmission.  All  the  filtering  operations  in  the 
transmitter  are  lumped  together  as  one  filter,  called  the  trans- 
mitter filter,  which  has  impulse  response  hT  (t)  . The  propaga- 
tion medium  is  represented  by  the  complex  time-variant  impulse 
response  h(t,x).  and  the  additive  noise  by  the  complex  precess 
n ( t ) . All  the  linear  operations  in  the  receiver  prior  to  the 
correlation  operations  are  lumped  together  into  one  receiver 
filter  with  impulse  response  h^(t).  The  output  of  the  filter 
is  fed  to  a number  of  parallel  correlators,  only  one  of  which 
is  shown  in  the  figure.  Complex  notation  is  used  to  describe 
the  correlation  operation  as  a multiplication  of  the  receiver 
filter  output  by  the  complex  conjugate  of  a shifted  probing 
signal  followed  by  a complex  low  pass  filtering  operation. 

Channel  play  back  is  achieved  with  a tapped  delay  line  as 
shown  in  Figure  12.  The  complex  time-varying  multipliers  are 
just  equal  to  the  low  pass  filter  outputs  of  the  correlators 
which  yield  estimates  of  the  impulse  response  for  uniformly 
spaced  values  of  delay  separated  by  seconds.  Filtering  opera- 
tions necessary  to  avoid  aliasing  in  the  digital  realization  of 
the  playback  channel  are  shown  in  dashed  lines. 
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Figure  11.  Definition  of  Basic  Signal  Processing  Operations 
in  Prober  System. 
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Figure  12.  Channel  Playback  Schematic  Diagram. 
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p.12.  4 lines  after  Eq.  (3.1):  Rj  = Rj  - R' . 
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W1 (a) . 
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Eq.  (4.52a):  delete  and  write  Cy  = (y-lja^ 
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p . 44 . 
(cont ' d) 


pg. 45 . 


p . 46  : 


Eq.  (4.53):  delete  right  hand  side,  y > k+2,  and 
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(y 


-i)  r(l+1)r(x^i) 

°o  ;•  y+r 
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2 Lines  below  Eq.  (4.53) : 
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amplitudes " . 
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ENVIRONMENTS 
I.  FOUNDATIONS *  * 


by 


David  Middleton** 
[Jan.  31,  1971] 


ABSTRACT 


Theoretical  foundations  for  an  analytical  model  of  urban  radio-noise  en- 
vironments are  presented,  at  a level  of  generality  broad  enough  to  include  the 
pertinent  physical  and  statistical  elements  which  critically  influence  the  tem- 
poral and  statistical  character  of  such  interference  in  radio  receivers.  The 
control  roles  of  geometry,  directionality  (beam  patterns),  waveforms  (from 
the  interfering  sources),  motion  (doppler),  source  density  and  distribution, 
and  secular  variations  of  source  parameters,  as  well  as  the  radiation  mechan- 
ism, are  specifically  developed.  The  basic  statistical  model  involves  indepen- 
dent sources,  in  space  (and  in  time),  leading  to  Poisson  radiation  fields  and  a 
Poisson  process  X(t)  in  a typical  receiver.  This  received  process  X(t)  can  be 
considered  the  sum  of  a Gauss  process  (by  the  Central  Limit  Theorem),  repre- 
senting the  cumulative  effects  of  a large  number  of  sources,  none  of  which  is 
very  large  vis-a-vis  the  rest,  and  a Poisson  process,  produced  by  those  few 
strong  transients,  which  when  present,  dominate  the  background.  The  process 
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Ill 


X(t)  is  essentially  stationary  for  periods  comparable  to  the  secular  period  of 
changes  in  traffic  intensity  and  flow,  which  permits  the  construction  of  usefully 
large  experimental  ensembles  from  which  to  estimate  the  process  statistics. 

An  equivalent,  semi-empirical,  but  more  analytically  tractable  model, 
similar  to  that  introduced  by  Hall  for  impulsive  atmospheric  noise,  but  used 
here  for  independent  sources,  is  constructed.  This  model  is  represented  by 
X(t)=a(t)  Z(t),  where  a,  Z are  independent  processes,  both  zero-mean,  and  Z 
is  regarded  as  gaussian.  The  first-order  p.d.  of  X(t)  is  determined  empiri- 
cally for  large  |x|  (i.e.,  |a|-*<»),  which  in  turn  specifies  the  parameters  of 
the  first-order  p.d.  of  a.  The  equivalent  model  is  postulated  to  have  the  same 
moments  (of  first-  and  higher-order)  as  the  basic  Poisson  model,  which  now 
permits  us  to  include  explicitfy  the  details  of  the  governing  physical  mechan- 
isms in  the  formulation.  The  present  study  serves  also  as  a guide  for  the  ex- 
perimentation needed  to  establish  the  model. 
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STATISTICAL-PHYSICAL  MODELS  OF  URBAN  RADIO-NOISE 

ENVIRONMENTS 

I.  FOUNDATIONS  * 


by 


David  Middleton** 
[Jan.  31,  1971] 


1.  INTRODUCTION 

The  present  paper  is  the  first  of  a series  of  technical  studies  whose  ulti- 
mate aim  is  to  provide  a joint  analytical  and  experimental  treatment  of  radio- 
noise  interference  occurring  in  the  urban  environment.  The  purpose  of  this 
treatment  is  to  determine  quantitatively  the  physical  character  of  such  electro- 
magnetic noise,  its  primary  sources,  its  effects  on  the  reception  of  desired 
signals,  and  methods  of  improving  the  design  and  performance  of  receiving 

systems,  in  order  to  mitigate,  minimize,  or  remove,  this  class  of  electro- 

1,  2 

magnetic  pollution.  This  is,  of  course,  not  a new  problem.  ’ However,  in 
order  to  achieve  these  desired  results,  and  to  take  advantage  of  the  many  ad- 
vances in  communication  theory  and  technology  of  recent  years,  much  more 
fully  developed  physical  and  analytical  models  of  these  noise  processes  appear 
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necessary  on  the  one  hand,  and  more  sophisticated  and  sensitive  criteria  and 
experiments  are  needed,  on  the  other,  with  theory  and  experiment  reinforcing 
each  other  at  each  appropriate  stage. 

The  aims  of  this  proposed  series  of  papers  may  now  be  stated  in  more 
detail.  We  wish  to  construct  an  analytical  model  of  typical  urban  radio-noise 
environments,  which  will  include  the  following  critical  features: 

(i)  the  geographical  (i.e.,  geometric)  distribution  of  the  various  inter- 
fering sources  vis-a-vis  a typical  receiver; 

(ii)  the  physical  character  of  the  noise  fields  generated  by  these  sources; 

(iii)  the  conversion  of  these  noise  fields  into  waveforms  upon  reception  by 
a typical  receiver; 

(iv)  a statistical  model  which  is  phenomenologically  broad  enough  to  ac- 
count for  the  principal  observed  interference  effects,  and  which  has 
sufficient  physical  structure  to  enable  us  to  identify  and  incorporate 
the  pertinent,  individual  interference  mechanisms; 

(v)  anafytical  evaluation  of  the  "statistical  parameters"  of  (iv),  which  in- 
corporate the  physical  structure  of  the  model.  Efy  "statistical  param- 
eters" here  we  mean  not  only  those  structural  parameters  of  the  fields 
and  waveforms  which  may  in  our  ensemble  model  have  statistical  dis- 
tributions, e.g. , amplitudes,  frequency,  phases,  epochs,  signal  dura- 
tion, etc.,  but,  principally,  the  "statistics"  of  the  noise  model.  These 
are  the  intensify,  the  covariance  functions  of  the  process  (and  their  as- 

3 

sociated  intensify  spectra),  moments  of  higher-order  and  higher  degree  , 
first-order  probability  distributions,  and  where  possible,  higher-order 
distributions. 

(vi)  various  special  features,  such  as  doppler  (sources  and  receiver  in  rel- 
ative motion;  directionality  effects  (beam  patterns);  mixed  classes  of 
interference,  i.e.,  a variety  of  different  interference  mechanisms;  the 
effects  of  buildings,  traffic  patterns,  diurnal  and  seasonal  changes, 
weather,  etc. 

The  noise  mechanisms  we  are  considering  here  are,  of  course,  man-made, 
and  as  such  must  be  identified  and  classified  as  to  their  relative  importance  as 
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interfering  agents.  Broadly  speaking,  we  may  say  that  these  man-made  proc- 
esses consist  of  a mixture  of  impulsive  and  background  noise,  with  the  former 
a markedly  impulsive  component  consisting  of  a few,  discrete  signals.  For 
the  present,  we  shall  focus  our  attention  on  ignition  noise,  which  appears  to  be 
the  most  common  culprit  where  communication  involving  land-mobile  trans- 
portation systems,  at  VHF  and  UHF,  are  concerned. 

Two  principal  approaches  are  indicated: 

(I)  measurement  of,  and  model-building  for,  the  wide-band  noise,  as  it 
actually  occurs,  and 

(II)  measurement  and  modelling  of  this  interference,  after  it  has  been  cor- 
rected to  comparatively  narrow-band  noise  in  a typical  receiver. 

4 

In  the  former  case  (I),  essentially  nothing  adequate  appears  to  have  been  done, 
while  for  the  latter  (II),  the  modelling  and  data  need  to  be  extended  to  include 
and  establish  the  validity  and  specific  role  of  many  of  the  physical  elements  cited 
in  (i)-  (vi)  above.  In  both  instances  more  sophisticated  models  and  criteria  are 
needed. 

Collateral  aims  of  this  series  of  studies  are: 

(1)  To  provide  structure  and  guidance  for  the  experimental  investigations 
that  must  be  undertaken  to  give  us  the  desired  empirical  knowledge  of 
urban  radio  noise.  This  is  needed,  in  order  to  relate  our  models  to 
reality  and  to  make  the  kind  of  measurements  which  will  enable  us  to 
predict  the  effects  of  this  electromagnetic  pollution  on  receiver  per- 
formance and  design. 

(2)  To  provide  a framework  for  simulation,  wherever  and  whenever  the 
analytic  results  are  intractable  and/or  the  reservoir  of  data  is  incom- 
plete or  too  expensive  or  impracticable  to  obtain. 

(3)  Once  adequate  noise  models  have  been  constructed  and  tested,  to  apply 
the  results  to  the  design  of  optimal  (or  near-)optimal  receivers  for 
maximum  immunity  against  this  class  of  interference.  This  includes 
the  calculation  and  measurement  of  receiver  performance  in: 

(a)  urban  radio-noise,  without  the  desired  signals; 

(b)  urban  radio-noise,  with  desired  signals. 
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The  new  results  of  this  study  (Part  I)  are: 

A.  The  construction  of  a primarily  physical  model  of  the  propagation  fields 
arising  from  the  spatially  and  temporally  distributed  noise  sources.  This 
provides  us  with  the  explicit  waveforms  at  the  receiver,  and  allows  us  to 
include  the  necessary  factors  of  geometry,  source  characteristics,  and  so 
on.  Earlier  approaches  have  not  usually  made  this  anatomization  of  their 
models.  We  emphasize  that  it  is  essential  to  include  the  waveforms  of  the 
interference,  since  it  is  ultimately  waveform,  i.e.,  a time-phenomenon, 
which  registers  upon  the  observer  as  interference. 

B.  The  development  of  the  so-called  "Basic  Statistical  Model"  (abbreviated 
"BSM"),  which  provides  the  framework  for  incorporating  quantitatively 
the  principal  mechanisms  and  features  noted  above  in  (i)~  (vi).  This  model 
consists  of  the  sum  of  two  independent  noise  fields.  The  first  is  a high 
density  Poisson  field,  with  (asymptotically)  normal  statistical  properties, 
representing  the  general  interference  background.  The  second  is  a low 
density  Poisson  field,  which  represents  the  few  but  distinct  and  effectively 
high  level  sources  of  impulsive  noise  interference.  From  the  physical 
description  of  this  composite  ambient  field  we  can  derive  the  detailed  struc- 
ture of  the  mean  intensity  and  covariance  functions  ( and  higher-order  mo- 

3 

ments,  and  in  principle,  the  probability  densities  as  well ),  which  consti- 
tute the  "statistical  parameters"  of  the  statistical  model  [cf.  (iv),  (v) 
above] . 

5 

C.  Theuseof  amodified  version  of  Hall’s  atmospheric  noise  model,  which  we 
shall  call  here  the  "Equivalent  Statistical  Model"  (ESM).  This  is  a plausible, 
quasi-empirical  model,  which  is  used  to  replace  the  exact,  but  generally  in- 
tractable Basic  Model  in  the  analytic  determination  of  probability  densities. 

In  essence,  it  replaces  the  sum  of  a Gauss  and  Poisson  process  by  a proc- 

g 

ess  equivalent  to  it,  at  least  in  the  first-order  distribution.  This  equiva- 
lent process  is  postulated  (and  later  justified)  to  be  a Gaussian  process  in 
product  with  a relatively  slowly- varying  process,  whose  first-order  dis- 
tribution is  asymptotically  matched  at  large  amplitudes  to  the  observed, 
empirical  distribution  (at  these  large  levels)  of  the  urban  noise  in  question. 
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The  key  innovation  here  is  that  we  now  determine  the  moments  of  higher 
3 

degree  of  this  quasi-empirical,  equivalent  process  directly  from  the 
corresponding  moments  of  the  Basic  Model  (BSM).  Thus,  for  example, 
the  covariance  function  of  the  former  is  specified  by  that  of  the  latter, 
which  is  itself  explicitly  structured  according  to  the  underlying  physics  of 
the  actual,  ambient  interference  field.  And  this  field,  in  turn,  we  have 
ab  initio  determined,  according  to  (A)  above.  [Certain  important  differ- 
ences between  this  Equivalent  Statistical  Model  (ESM  ) and  Hall's  Atmos- 
5 

pheric  Model  will  be  discussed  more  fully  in  the  following  sections.] 

The  approach  described  in  (A)~  (C)  above  is  in  part  based  on  recent  work 
7 8 

of  the  author  ’ , with  appropriate  modifications  and  extensions  ( see  Sections 
3 & 4 below).  The  organization  of  this  paper  (Part  I)  is  accordingly:  Section  2 
outlines  the  principal  steps  in  the  general  formulation  of  our  models;  Section  3 
is  devoted  to  the  physical  aspects,  and  Section  4,  to  the  statistical  features,  of 
both  the  BSM  and  ESM.  Section  5 concludes  with  a brief  discussion  of  the  prin- 
cipal results  and  indicates  some  of  the  next  steps  to  be  taken.  In  Part  I here, 
our  presentation  is  what  might  be  called  a "detailed  summary,"  where  we  re- 
serve to  later  papers  (Parts  n,  in,  etc.)  any  desired  developments  in  depth. 
Finally,  although  we  consider  only  ignition  noise  specifically  here,  it  will  be 
clear  from  Secs.  2-4  that  this  in  no  fundamental  way  restricts  the  generality 
of  our  treatment. 


2.  STEPS  IN  THE  FORMULATION: 

We  next  outline  the  principal  stages  in  the  construction  of  our  statistical 
model.  A somewhat  more  detailed  account  is  given  in  Sections  3 and  4,  follow- 
ing. Our  intent  here  is  to  provide  a general  overview  before  developing  the 
analytical  structure. 

We  begin  with  the  determination  of  the  interference  field  N_(  t,  R)  and  the 

waveform  X(t)  of  the  received  interference,  after  it  has  been  processed  by  the 

(N) 

receiving  aperture  (antenna),  whose  beam  pattern  is  £7  ( v,  f).  Letting  T , 

g •h«iVI 

T be  the  general  operators  representing  respectively  the  effects  of  the 
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medium  — here  the  set  of  interfering  noise  ( = signal ) sources 
the  receiving  aperture  — we  may  write  generally 


(2.1) 


for  the  received,  "ambient"  noise  entering  the  subsequent  linear  (and  nonlinear) 


In  this  way  we  can  link  the  various  stages  in  the  generation  and  propagation  of 
the  interference  field  to  the  final  observed  effect.  For  the  present,  however, 
our  attention  is  directed  toward  the  waveform  X(t). 

th 

We  proceed  as  follows:  we  consider  a typical  j — source  of  ambient  radio 
noise  Sj(t,R)^  and  calculate  its  field  at  some  point  P(R).  We  assume,  for  the 
time-being  here  to  simplify  the  analysis,  a scalar  source,  i.e.,  one  that  effec- 
tively represents  that  particular  component  of  the  actual,  vector  field  that  the 

receiving  aperture  principally  responds  to,  so  that  the  resulting  field  is  a 
12 

scalar  quantity.  Next,  we  additively  combine  a number  K (j=l,  ...,K)  of  the 
fields  of  these  sources  to  produce  the  composite  field  that  impinges  upon  our  re- 
ceiver. These  sources  are  independently  operative  and  distributed  in  space  — 
for  our  purposes  of  mobile-land  transportation  (e.g. , ignition  noise)  — on  a 
surface,  namely  the  streets  of  the  urban  region  in  question.  Knowing  T , we 
then  have  the  desired  resultant  waveform,  X(t),  cf.  (2.1).  By  suitable  attention 
to  the  form  and  structure  of  S (t,R).,  and  to  the  relative  locations  of  these 


sources,  in  the  particular  region  and  under  the  conditions  appropriate  to  the 
traffic  state  at  the  time  (of  day,  week,  year,  etc.),  we  can  incorporate  the  gov- 
erning physical  and  geometrical  characteristics  of  the  interference  situation. 

This  completes  the  first  major  step  in  the  construction  of  the  model.  The 

second  major  step  is  to  introduce  a statistical  mechanism.  This  is  done  by  pos- 

8 

tulating  that  the  sources  have  a Poisson  distribution  in  space.  Moreover,  each 
source  represents  an  "event"  in  time:  the  generation  of  a noise  field  whose 


(N) 

stages  T' 
~K0 


of  the  receiver.  Thus,  the  observed  interference  is 


(N)  (N) 

= Rq  =AR£M 


(2.2) 
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physical  character  is  basically  transient  and  impulsive,  the  epoch  of  initia- 
tion of  the  transient  being  a random  parameter  vis-a-vis  the  receiver.  Other 
random  parameters  may  be  introduced:  amplitude,  frequency,  duration,  for 
example,  which  allow  us  to  treat  a variety  of  ambient  sources,  which  may  be 
operating  concurrently,  but  which  may  arise  from  different  mechanisms.  In 
any  case,  the  result  is  that  the  interference  field  N^(t,  R)  is  a Poisson  stochas- 
tic field,  and  that  X(t)  is  correspondingly  a Poisson  random  process,  since  the 

aperture  ( ~T  \ operates  linearly  on  the  incident  field. 

^ A II  J 

A further  structure  can  be  imposed.  From  our  experience  so  far,  it  is 
reasonable  to  classify  the  noise  sources  into  two  groups:  those  that  are  not  re- 
solvable in  the  receiver  and  which  together  form  a "background”  noise  field,  and 
those  that  are  sufficiently  strong  to  appear  individually  and  discretely.  The  for- 
mer are  very  numerous,  of  all  signal  strengths  and  locations,  but  none  such  that 

anyone  dominates  the  others  vis-a-vis  the  receiver.  The  resultant  is  then  essen- 

14 

tially  a "high-density"  Poisson  process,  whose  statistics  are  asymptotically 

gaussian.  The  latter  few,  however,  retain  their  individual  character  and  remain 

14 

a Poisson  process,  of  low-density,  and  may  exhibit  the  very  large  amplitude 
excursions,  when  "on,"  of  typical  impulse  noise.  Thus,  we  have  a sum  of  gauss 
and  Poisson  fields  and  processes: 


N][(t,R)  = G(t,R)  + P(t,R)  , 


(X(t)}  = G(t)  + P(t)  , 


(2.3a) 

(2.3b) 


respectively. 

The  key  assumption  here  is  the  independence  of  the  various  sources  with 

respect  to  one  another.  [This  is  entirely  reasonable  for  the  urban  noise  model, 

where  one  ignition  acts  independently  of  another,  but  is  generally  not  the  case 

for  atmospheric  noise,  where  the  interference  tends  to  come  in  clusters  of  im- 
15 

pulses,  so  that  the  purely  Poisson  assumption  is  no  longer  valid.  ] If  we  let 


U ( j = 1, ... , K ) be  the  j — waveform,  attributable  to  the  j — source,  as  it  appears 
in  the  receiver  following  the  aperture,  we  can  accordingly  divide  X(t)  into  these 
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two  components: 

K S S 

x«-£vl>  Bi+EuDj’  (2-4) 

3 3 3 

where  the  first  is  the  sum  of  the  large  number  ^K^-*  of  relatively  weak  and 
comparable  sources  producing  the  gaussian  background,  and  the  second  repre- 
sent the  few  0,  1, 2,  3,  say  j , strong  sources  which  discretely  impinge  on 

the  receiver,  and  which  belong  to  the  low-density  Poisson  process.  This,  then, 
is  the  essence  of  our  nBasic  Statistical  Model"  (BSM),  which  we  shall  regard  as 
the  fundamental  model,  by  which  all  other  equivalent  and  approximate  models 
are  to  be  calibrated,  as  we  shall  indicate  presently.  Statistically,  the  BSM  [for 
X(t)]  is  described  as  the  sum  of  independent  Gauss  and  Poisson  processes. 

Although  the  BSM  appears  to  be  the  correct  form  of  noise  model  here,  it 
presents  a fundamental  technical  disadvantage:  it  is  analytically  complex,  when 
we  attempt  to  develop  explicit,  closed-form  expressions  for  the  probability  den- 
sities describing  the  process,  densities  which  are  needed  not  only  in  estimating 
and  predicting  its  properties,  but  also,  ultimately,  in  determining  the  perfor- 
mance of  receiving  systems  subject  to  such  interference. 

One  approach,  which  we  shall  consider  in  the  course  of  the  investigation, 
is  directly  to  employ  simulation,  based  on  the  general  expressions  developed 
below  (Sec.  4)  for  the  BSM.  A second  approach,  which  is  very  attractive  ana- 
lytically, is  to  use  an  "equivalent,"  mixed-process  model  (the  ESM),of  the  type 

5 

suggested  by  Hall's  Atmospheric  Noise  model,  but  differing  from  it  in  a number 
of  basic  assumptions.  Specifically,  we  replace  (2.3a,b)  by  the  mixed  (or  "mod- 
ulated") process 


{X (t) } = fa(t)  Z(t)} 


(2.5) 


where  a(t)  is  a slowly- varying,  comparatively  broad-band  random  process,  which 
effective^  sets  the  unit,  or  scale,  of  X(t),  and  whose  statistical  properties  are 
to  be  determined  empirically.  To  begin  with,  this  is  done  by  choosing  W^(a),  the 
first-order  p.d.  of  a,  so  that  it  gives  the  same  first-order  p.d.  of  X for  large 
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amplitudes  (|x|  -» <»),  as  is  observed  empirically.  Higher-order  densities  are 
to  be  similarly  constructed,  from  whatever  higher-order  statistical  infor- 
mation can  be  obtained,  in  conjunction  with  the  physically  derived  moments  of 
the  Poisson  process. 

Here  Z(t)  is  a (zero-mean)  gaussian  random  process,  independent  of  a(t), 

whose  statistical  character  is  to  be  found  from  the  basic  Gauss-plus-Poisson 

model  (BSM)  above,  and  whatever  empirical  data  are  required  to  relate  that 

model  to  specific  noise  situations.  As  we  shall  show  presently  (Sec.  4),  at 

least  the  expected  first-order  p.  d.  of  X is  gaussian  for  small  amplitudes,  which 

1 6 

supports  our  postulate  of  the  gaussian  character  of  Z(t).  For  the  important 
first-  and  second-order  moments  we  have 

(X>  = (a(t)>  • (Z(t)>  = 0 ; <Z(t)>=0,  <a(t)>  = 0 (2.6) 


for  the  mean  values  (this  last  from  physical  considerations,  Sec.  4),  and  for 
the  covariance  functions, 


Kx(*i-,a)  “ <x(ti)x(t2))  = *JS>h)  kz(V2)  = Kx(trt2)  • <2-7> 


BSM 

(<x>=0) 


where  p is  the  covariance  function  of  the  BSM,  whose  detailed  structure 
is  determined  by  the  Poisson  form  of  the  process  and  which  explicitly  incorpo- 

17 

rates  the  basic  waveform,  U (e.g. , the  pltysics)  of  the  interfering  noise,  e.g. , 


Kx(*i,t2)ggM  = f • 


(2.8) 


Here,  p is  a "density"  function  which  measures  the  number  of  (independent) 
contributing  noise  sources  in  the  domain,  A , of  possible  sources,  which,  in 
turn,  is  determined  by  the  particular  urban  environment  we  are  considering. 
The  ^ (=  t,  9,  cp)  are  a set  of  coordinates  with  respect  to  a convenient  (at  rest) 
reference  system,  and  cP.  is  a volume,  or  surface  element  [dtdcp,  dtdGdcp],  cf. 


Sec.  4-  (5)  below.  The  9 are  a set  9 = 


W-] 


of  (time-invariant)  statistical 
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parameters  which  can  govern  the  shape  and  strength  of  the  basic  wasveform  U. 
In  view  of  (2.4),  we  see  that  (2.8)  can  be  further  anatomized: 


Ub(VM))  S 

B 

+ J Ud(VM))  Si  - 

ad 


(2.9) 


where  are  the  respective  Poisson  densities  for  the  background  (limiting 

gaussian)  process,  and  the  low-density,  transient  interference;  (it  may  be  that 
UB=  without  loss  of  generality),  and  Ag  UAd.  Accordingly,  if  we  can  de- 
termine K , we  obtain  K,_  , cf.  (2.7). 
a’  Z ’ 

Note  that  the  statistical  independence  of  a(t)  and  Z(t)  requires,  at  least, 

18 

that  their  (intensity)  spectra  do  not  noticeably  overlap.  Actually,  here  a(t) 
will  be  quite  low-frequency,  i.e.,  slowly -varying,  while  Z(t)  is  comparatively 
narrow-band,  centered  about  some  carrier  frequency  f , even  when  apertures 
and  RF  stages  of  the  receiver,  broad  compared  to  the  usual  communication 
practice,  are  used,  e.g.,  0(1,2  Mc/sec.)  vs.  0(10  kc/sec. ).  The  process  X(t) 
is,  therefore,  "narrow-band,"  and  can  be  expressed  as 


X(t)  = a(t)x(t)  cosaUgt  + a(t)y(t)  sinu^t  ~ |a(t)|  e(t)  |cos  u^t  - cp(t)J  , (2.10a) 

e ^ 0 

= E(t)  cos  [.Qt  - cp(t)]  , (2 . 10b) 

where  E(t)  = | a (t)  | e(t)  ( ^ 0)  is  the  envelope  of  X(t),  and  cp(t)  =tan  ^ [ a (t)  y(t)/ 
a(t)x(t)]  = tan  "*’[y(t)/x(t)],  is  the  phase.  From  (2.5),  we  have 


Z(t)  = x(t)  cos  u)^t  + y(t)  sin  u^t  , 


(2.10c) 


and  it  is  clear  that  since  Z(t),  and  x(t),y(t),  are  gaussian  processes,  the  phase 
cp(t)  must,  in  the  first-order  be  uniformly  distributed  over  2rr . The  envelope  e(t) 
will  likewise  have  a Rayleigh  distribution  in  the  first-order,  but  the  envelope  E(t) 
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will,  in  general,  not.  Even  if  the  bandwidth  of  Z is  comparable  to,  say,  f^ /2, 
this  model  applies,  but,  of  course,  it  is  then  physically  meaningless  to  speak 
of  X(t)  as  having  an  envelope  and  phase  in  the  customary  sense,  cf.  Eqs.  (2. 10). 

So  far,  we  have  neglected  system  noise  ( which  is  symbolically  indicated 

(N)\  ' 

in  Trq  ] • This  is  permissible  as  long  as  the  background  ambient  field  of  the 

interference  is  strong  enough  to  dominate  the  system  noise  in  the  receiver, 
which  may  often  be  the  case.  However,  if  the  background  contribution  is  exces- 
sively weak,  and  only  the  transient  interference  is  significant,  then  we  must  in- 
clude the  system  noise.  Since  this  is  gaussian  (and  stationary),  with  a spectrum 
determined  essentially  by  the  frequency  response  of  the  RF(-IF)  stages  of  the 
receiver,  we  have  again  our  Basic  Model,  but  now  only  the  impulsive  interference 
belongs  to  a Poisson  process.  Equation  (2.9)  is  accordingly  modified  to 


with  the  covariance  of  the  system  noise  becoming 


KS(T) 


*WX)  hRI(X  + T)dx  » 


’ 


(2.11) 


(2.12) 


cf.  Eq.  (3.86),  Ref.  10,  where  is  the  spectral  intensity  density  of  the  ("white") 

system  noise  initially,  and  h is  the  combined  weighting  function  (impulsive  re- 

^ 19 

sponse  of  the  RF-  and  IF-stages  of  the  receiver  ). 

We  have  now  the  beginnings  of  a quantitative  theory  of  urban  radio-noise 
interference.  Two  major  steps  follow: 

(1)  we  must  develop  the  analytical  details  of  our  models;  and 

(2)  we  must  verify,  and  modify  the  theory  where  necessary,  by  appropriate 
experim  entation. 

For  the  former  we  sketch  the  development  in  Secs.  3 and  4 below;  for  the  latter 
a program  of  experiments  is  now  underway  at  the  Transportation  Systems  Cen- 
ter (TSC) . 
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3.  AMBIENT  FIELDS  AND  WAVEFORMS: 

Here  we  shall  indicate  the  principal  steps  in  determining  the  interference 
field  incident  upon  the  receiver  and  the  resulting  noise  waveform  following  the 
receiver's  aperture.  The  present  approach  is  based  mainly  on  the  methods  of 
Ref.  7 [cf.  Sec.  IE  B). 


A.  General  Development: 


21 


We  begin  by  observing  that  the  principal  component  of  the  electro- 


.th 


magnetic  field  produced  by  a typical  j — source  S^  obeys  the  (scalar)  wave  equa- 
tion (in  rectangular  coordinates  and  for  an  assumed  homogenous,  isotropic,  and 
non-conducting  medium  here) 

1 — _ 


v2^, 


at 


- (ievi) 


(3.1) 


= 0 


elsewhere 


where  Gj^  is  the  source  function,  and  c ( = 3 • 10  m/sec.)  is  the  velocity  of  prop- 
agation; Vj  denotes  the  region  occupied  by  the  source  and  is  generally  a volume. 
The  geometry  is  shown  in  Fig.  3.1;  P(R')  is  some  point  in  space,  a vector  dis- 
tance R*  = -R'-R.  from  the  source's  origin  0„,  andhere  R.  is  the  vector  joining 

J J , ^ J a ^ 

the  origins  of  the  j—  source  and  the  receiver, with  R.  = i 'x'  + i 'y#  + i 'z#  = 

A A A ~X  ^ ^ 

i x +2  y + i z,  etc. ; d£  and  dTj  are  radiating  and  receiving  elements  in  their 
*"x  ~y  -z  — 

respective  apertures,  and  similarly  oriented  rectangular  coordinate  systems 
are  used  for  conve: 
aperture,  we  have 


are  used  for  convenience.  Letting  Aj(£,f)^  be  the  (frequency  form)  of  the  source's 


f °° 

Gi('>lJ2j). = J Ai<i- f >j  si(f>ll2j). ® » <“  = 2nf ) . 


(3.2) 


where 


JCD 

Sl(t,?|Rj)  e dt  (3.2a) 

is  the  amplitude  spectrum  of  the  j—  source  Sj^t,£  |R^  , which  in  generality  i 
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Figure  3.1.  Geometry  of  interfering  source  and  the  receiver 
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position-dependent  and  unless  otherwise  indicated,  is  a real  quantity.  An  alter- 
native form  for  At.  is 
1 


—hut  , 
e dr  , 


(3.3) 


so  that  (3.2)  is  also  written 


3 “«  3 


dr  . 


(3.4) 


Thus,  the  source's  aperture  A^  ^or  is  a linear  filter,  which  also  depends  in 
magnitude  on  location  (§).  The  spatial  Fourier  transform  of  the  aperture  is  the 


beam  pattern  <7  , e.g., 


a. 


C2,f>j  = f A (5,f)  e2ni~TdS  , 
J J y 1 ■«**  J 


(3.5) 


where  v is  a (vector)  spatial  frequency; 


A A A A A A A -1 

v(i)=i  v +i  v +i  v = k/2n  = i X + i X + i X , (3.6) 

~ — x x — y y ~z  z — ~x  x -y  y ~z  z ’ 

-1  A A /i. 

with  X^  , etc.  wavenumbers,  for  the  direction  i of  P(R),  e.g. , ^ H R/ |R  | , 
where  P(R)  is  some  point  in  space,  outside  V . Applying  (3.3),  we  see  that  the 
beam  pattern  is  alternatively  given  by 


W.L 

V — CO 


-2nifT  . 
e dT 


(3.7) 


Note  that  like  the  aperture  function  A^,  the  beam-pattern  Cl^  is  frequency- 


dependent. 


23 


We  assume,  safely,  thatP^R^,  P^R(  j are  always  in  the  far-field,  so  that 


the  solution  of  (3.1)  becomes 
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r ro  iu/t-R?/c\ 

43).^/  fi'i'(ij,f/c>f|si).e  J df  - 

1 ] v '3 


or 


or 


-ssrf  i5J"oiCi’t-T-R;/c+-i-i/0).si(T’l|£i).dT  - 

3 jy  ^ —a>  'j  N J 

„ r f oo  itu/t-R'  /c  + | • i.  /c\ 

^ sf  1 J ~“J  J 

3 •'V  •'-»  3 


(3.8a) 


(3.8b) 


df  , (3.8c) 


(3.9) 


where  ^ is  a generalized  beam  pattern: 

i7i(^t/0,f|si)j  s=s|Ai(I,f)j  si(f,Il£j).|  • 

The  spatial  frequency  v(^ )isj'f/ c = (r'/r*  \ f/c,  R'  = | R?  | , etc . In  the  rectan- 
gular  coordinates  of  the  source  we  have 


a.  a , a a 

i = i cos  cp  sin  0„  + i sin  cp„  sin  0_  + i cos  . 

— »-v  Ts  s —y  Ts  s —7.  s 


(3.9a) 


24 

Reception  is  handled  in  an  analogous  fashion.  We  have  for  the  typical 

waveform,  U.,  after  processing  by  the  receiving  aperture,  A (T],f),  with 
j R 

R{  “*  r (T|): 

"“3  “*3  ~ 


uj(t)=J  £D  J Vljt"T)s£[T^j(l)]  dT  » 


or  in  frequency  form 


(3.10a) 


uj(t)  =JV  13  J."  V3>f>  s^[flij(l)]  dt 

R 


(3.10b) 


where  S is  the  Fourier  transform  of  the  field  a?.,  (3.  8),  in  which  r.(T])  =R.-T|, 
cf.  Fig.  3.1.  As  in  (3.3),  (3.4),  (3.8b),  is  the  spatio-temporal  form  of  the 
receiving  aperture,  and  V is  its  physical  domain.  Inserting  (3.8)  into  (3.10) 
gives  us  a variety  of  equivalent  expressions  for  the  waveform  of  the  interfering 
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(N) 

noise  entering  the  time-processing  portion  of  the  receiver,  , cf.  (2.2). 
These  are,  for  this  far-field  state  (when  there  is  ignorable  mutual  coupling 
from  receiver  to  source): 


dT 


or 


4nR.  j m ^ -Rj  VC’ f ) ai\ (-R)  !/c 


',dR.\  dT' 

(3.11a) 

io )( t-R.  / c^ 

3 ' ]/  > df  , 

(3.11b) 

A A 

>5f  2J 4 

have  used  i -*  i_  , r.  = 
- ~R’  j 

R.  + i • 7]/c  in  this  far-field  approximation  along  with  the  definitions  (3.5),  (3.9). 
3 ~-R  ~ A 

[We  remember  that  here  is  the  unit  vector  directed  outward  along  R toward 
O of  the  receiver.] 

Note  that  the  operators  for  the  j — source  and  for  the  receiver,  on  com- 
paring (3.10),  (3.11)  with  (2.1),  can  be  written  now  explicitly 


rjv  >T, 


dT 


or 


U(  )f;jeto(t)df  , 


(3.12a) 


(3.12b) 


and 
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T(N){  } 1 

-M  J‘ 


or 


• 1 I I t-R.  /c  + i_  • ?/c\  ( ) 

J 4"Vv  ~J-«  V ]/  ~Ri  - J 

‘"Hal 


. dT  (3. 13a) 


Wj ( >f.j 0 


i 

l f°°  icr/t-R.  / c^ 

1 j/  >At  ■ 


2rrifiR  • j/c 


d? 


(3.13b) 
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Equations  (3.12),  (3.13)  are  the  associated  far-field  (Fraunhofer)  oper- 
ators. The  condition 

AtL2  A . ) « R.  (3.14) 

\ max/  min/.  j 

must  be  obeyed  here  [and  for  (3.8),  (3.11)],  where  L is  the  largest  physical 

xricLX 

dimension  of  the  interfering  source's  aperture  AT  , and  X-  . is  the  shortest  wave- 

I’  min 

length  of  the  significant  radiation  from  the  source.  For  example,  atf  ^100mc/sec. 

2 ^ 

and  L =1  meter,  we  have  tt  • 1 /3  « R,  or  R » 1 meter  for  far-field  opera- 
max  7 7 

tion. 


B.  Stochastic  Field  N^(t,R)  and  Random  Process  X(t): 

The  stochastic  field  representing  the  total  effect  of  all  interfering  sources, 
considered  over  all  possible  representations  of  such  sources  in  space,  time, 
and  parameter  space,  is  simply  the  sum  of  the  individual  contributions 


Nl(t,R)=j£ 


£ 


(*■  ■>), } ■ 


afL  given  by  Eqs.  (3.8) 


(3.15) 


Again,  for  the  far-field,  (3.8)  is  used,  where  now  r.=  R-R. , is  the  vector  distance 

tli  ^ ^ 

between  the  j — source  and  some  point  P(R)  in  the  field  [here  R and  R are  both 

directed  toward  the  origin,  O , cf.  Fig.  3.1] . In  operator  notation  (3.15)  be- 
comes 


Ni(t’£>  ; Im>[  } 1 ! 


(3.16a) 


E^*’ ^ E% l/r  A1  t/C’ ' 'Si< f’ "J6 

J J 


ico(t-r./c) 


df  • 
(3.16b) 


Here  ir^  is  the  imit  vector  - ^R- R^^/|R- lb  | , directed  toward  P(R).  Figure  3.2 
illustrates  the  situation  schematically,  for  a given  configuration  of  sources 
during  any  short  period  of  time. 
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Figure  3.2.  Schema  of  a typical  distribution  of  urban  noise  sources 
and  a receiver  (R) 
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As  mentioned  in  Sec.  2,  we  may  divide  the  interference  field  into  a back- 
ground, or  gaussian  field  G(t,R)  (as  the  limit  of  a Poisson  field),  and  the  added 
field  P(t,  R)  of  a small  number  of  strong,  discrete  sources,  cf.  Eq.  (2.3a).  The 
statistics  of  N(t,R)  may  then  be  determined  by  the  methods  outlined  in  Sec.  4 
following.  The  important  thing  to  observe  here  is  that  the  physics  and  structure 
of  the  source  enter  explicitly  through^  and  Sj^  above. 

Our  principal  interest,  however,  at  this  stage  of  the  study  is  the  inter- 
ference process  X(t)  in  the  receiver.  This  is  directly  die  sum 


X(t)  = {E  lk(t)j  ; j\h  given  by  Eqs.  (3.11)J  , 
which  in  operator  formalism  becomes  [ cf.  (3.12),  (3.13)  in  (2.1)] 


(3.17) 


(3.18a) 


■El, 

] 


v- 


(-v/c’f)^b} 


iUjt 

e df 


(3.18b) 


with  once  more  given  by  (3.8),  R'  ->  R^  , etc.  As  Eq.  (2.4)  indicates,  we  may 

divide  X(t)  into  two  processes,  corresponding  to  the  gaussian  background,  G(t), 

and  low-density  Poisson  process,  P(t),  of  the  few,  interfering  sources,  strong 

in  the  receiver,  cf.  Fig.  3.2.  Again,  the  physics  and  the  structure  of  the  source 

enter  explicitty  here,  through  sd.  and  Sj^. 

3 •* 


C . Generalities  and  Modifications: 

We  need  next  to  develop  the  model  somewhat  further,  and  to  include  such 
features  as  doppler  due  to  relatively  moving  sources  and  receiver,  narrow- 
bandedness, transient  signal  structures,  and  generalizations  of  the  source  mech- 
anism. We  include,  also,  an  Example. 
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(1)  Popp]  er  : 


Because  of  the  possible  relative  motion  between  a typical  interfering 
source  and  the  receiver,  a doppler  distortion  in  the  received  waveform  is  in- 
troduced; (this,  of  course,  affects  desired,  as  well  as  undesired  signals).  We 
may  calculate  this  doppler  by  a direct  adaptation  of  the  procedure  of  Sec.  4.4, 
Ref.  8.  Letting: 

t = time  at  which  leading  edge  of  the  wavefront  (from  source  j ) reaches 
the  receiver; 

t = time  at  which  leading  edge  of  the  wavefront  (from  source  j ) is  emitted 
b 

from  the  source, 

we  find  that  here,  after  modifying  the  Eqs.  (4.22a-c)  in  Ref.  8 appropriately, 


*R 


R/c  + (1~B sr)  \ 
1 + 6 


RS 


(3.19) 


The  delay  between  source  and  receiver  is,  therefore, 


V‘s‘B 


/,  r \ ' R/e-‘s(BSR+BRs) 

(R>  s)/  i+eRS 


(3.20) 


Since  t = t-T,  where  t is  the  "now"  at  the  receiver,  and  t is  the  elapsed  time 
b 

since  transmission,  we  get 


delay:  R^(t,  t)/c  = 


R/c  + (ft SR  + ^ks)  C 

1+P, 


RS 


a.  + (T-t)b.  , (3.21) 

J J 


for  the  now  time-dependent  delay  R ^/c  appearing  in  (3.11a).  Here  we  have  the 
normalized  doppler  velocities 


(3.22) 


where  v , v are  respectively  uniform  velocities  of  the  source  and  receiver  with 

respect  to  a stationary  coordinate  system  (the  same  in  which  the  velocity  c 

is  measured);  v and  v need  not  be  coplanar,  but  thej7  must  be  uniform.  For 
—v,  ~-R  2^ 

the  small  platform  and  receiver  velocities  encountered  in  practise  (3.21) 
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reduces  to: 


(small  doppler):  R.( T,  t)/c  = TQ.  + (pgR+  PRg)  ( t-  t)  ; 


tosVc(1+0hs)"Vc  • 


(3.23) 


It  is  convenient  to  define 


3+3  / 

?SrJW  +p  x . 
i + 3_  V psr  sp.y  ’ 


RS 


H = l + e = ( l - 3) 


SR  SR 
1 + 3 

-1  hRS  . , D \ 

1 ~ 3gR  (PRS+PSr) 


(3.24) 


g — ■ 3 +3  — 3 

RS  SR 


.*.  1 + 3 = p. 


for  each  j—  source  vis-a-vis  the  receiver. 

Let  us  now  "dopplerize"  our  expressions  (3.11)  and  (3. 12).  For  (3. 11a)  we 

simply  replace  R. /c  therein  by  R.(t',t)/c=T  +(t,-t)  3.,  and  the  temporal  ar- 

3/  ii  J i ^ 3 3 

gument  of  becomes  ( t - t')^1  + 3^  J - T +jRj  ’ §/c.  For  (3.  lib)  we  get 

!/*  00  A a -icoT  +i(l+3)fot 

J ^R('lR<1+B)£/C’<1  + E1)f)Cl(-Bt/C>flSI[(1  + P)f’l])e  df 

(3.25) 

and  the  operators  (3.12b),  (3.13b)  become 


1AR(  }j={l_/R("V  + P)f/C><1  + P)f)(  ’l* 


(1  + 3)  (at 


df 


and 


_tx  / , r co  i to [(1+3)  t - T ! 

Im  ' ^ si  *i'(v/°’fl(1+P)f)e  df 


(3.26a) 


(3.26b) 


With  no  relative  motion,  3 = 0 (|i  = l,  e=0),  T^-*R/c  and  (3.25),  (3.26)  reduce  as 


expected  to  our  previous  results  (3.11)-  (3 


.13).  [ 


In  the  usual  case  of  small 
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doppler,  with  comparatively  nari'ow-band  Sj^ , we  may  set  3 = 0 in  G^ 


R' 


(2)  Generalized  Sources  and  Signals: 

Wliile  the  class  j Sj(t,  %)  j.  of  sources  is  probably  adequate  to  account  for 
many  of  the  source  mechanisms  that  produce  the  urban  noise  interference  field, 
we  can  generalize  it  to  include  a much  broader  class,  by  making  time-vari- 
able, i.e.,  setting  Sj=  ^t,  t;  ^ |R^  , etc.  Here,  the  convention  is  to  regard 
"t",  the  first  variable,  as  the  ''memory"  and  "t",  the  second  variable,  as  the 
time- variation  of  the  waveform.  Thus,  the  source  function  Gj,(3.4),  now  be- 
comes 


(3.27a) 


= J Ai<I>f)jSI(£,t;S|R.)  eW,tdf  , (3.27b) 

where  now  [ in  (3.27b)]  is  the  time-variable  amplitude  spectrum  of  S^( t,  t; ...),. 
The  relations  (3. 8),  (3.9),  (3.10),  (3. 11a)  are  formally  unchanged,  except  that 
G'  is  now  time-dependent,  and  depends  on  t.  The  alternative  form  (3.11b), 
however,  is  modified,  as  is  (3.25)  when  there  is  a doppler.  The  result  for  (3.25) 
becomes 


DI(t)t**!j  *R(4f/c>f)eil,tdfJ  d2EAT<L’f'>J  V[[(i  + P)f'if-(1+P)£':£|R) 

j -00  ' Vj  -00 


2tt  if” 


X e 


ifW 


c -icu'T, 


• e 


0 


df' 


(3.28) 


27 


where  ^ is  the  bi-frequency  function  of  the  time-varying  source  S^,  viz 


(3.29) 
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jwhen  Sj  is  not  time-variable,  then  ^ f(l  + 3)  f';  f- (1  + 3)  f',  § |rJ  = Sj  [(1  + P)  f'; 
5|R]  6 [f-  (1  +P)  f'] , and  (3.28)  then  reduces  directly  to  (3.25);  with  no  doppler, 
VI  = VI(f,,f-f/;?|R)-]  The  operators  (3.26a,b)  are  readily  obtained  from  (3.28) 
by  inspection. 

To  facilitate  handling  of  the  general  class  of  transient  signals  , we 
may  replace  the  various  Fourier  integrals  over  frequency  in  (3.2)-  (3.28)  by 
suitable  contour  integrals,  e.g., 

co i+d  ( >0)  , 

Gl(t;ij2j).  AI(!,»/2"i)jSI(s/2T.i,tiS|Hj)^  , (3.30) 

28 

etc.,  where  (-ooi+d?  »i+d)  is  a Bromwich  contour. 

When  instead  of  the  distributed  sources  assumed  above  we  may  use  the 
approximation  of  a point  source,  then  we  have 


or 


SI(T’ll5}  = Si(Tl^}  6(I"0)  > 

SI(T,t;£|R)  = SI(T,t|R)  6(£-0) 


(3.31a) 

(3.31b) 


for  the  time-varying  sources  above.  The  (far-)  field  s£  , (3.8),  is  accordingly 


- a?  iu/t-R(  /c\ 

(t>£j).^J.ojAi(0>t)isi(flsj)ie  /df  > 


(3.32) 


and  is  thus  omni-directional:  <5^(0,  f)  =Aj(0,  f)^  here.  Similarly,  we  obtain  for 
the  received  wave  (3.11b) 


(pt.  source):  U 


X^(0,f).  S^f|R.)  e 


df 


(3.33) 


with  corresponding,  obvious  modifications  of  the  Eqs.  (3.12),  (3.13).  Extending 
the  point-source  assumption  to  the  cases  of  doppler  and  time-variable  sources, 
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(3.31b),  we  find  at  once  that  the  desired  waveform  (3.28)  in  the  receiver  be- 
comes 


(pt.  source): 


-iu/T 

Aj(0,f').  V[[(l+B)f';f-(l+3)f/|R]j  e °df' 

(3.34) 


(N) 

with  corresponding  forms  for  the  operators  T^j,.  ,T^,j. , cf.  (3.26a,b).  [Note 
from  (3. 24)  that  for  small  doppler,  1 +P  = g.] 

Finally,  there  remains  the  structuring  of  S (t,5  |R).  While  the  details  re- 

1 /Vv, 

main  for  verification  by  experiment,  a plausible  model  of  ignition  noise  is  (for 
point  sources  here)  the  deterministic  waveform 


S ( 1 1 R).  = A cos  u>  Yt-e.^  , 
I '-J  OJ  OjV  ]/  ’ 


to<t<to+Ts; 


= 0,  elsewhere 

(3.35) 

where  A^,  (ju^,  are  time-invariant.  The  epoch  (with  respect  to  a cycle 
of  the  carrier  fg)  indicates  a time  relationship  between  the  initiation  of  the  igni- 
tion transient  (3.35)  and  the  time  scale  of  the  observer  at  the  receiver.  The 
time-delay  between  reception  and  the  emission  of  the  transient  field  from  the 
j—  source  is  given  by  T cf.  (3.23).  This  latter  converts  a distance  /ft.  'j  into 
a time-delay  iT  j by  virtue  of  a finite  velocity  of  propagation  (c),  while  the 
former  establishes  the  coherence  or  incoherence  of  the  ignition  waveform 
with  the  observer’s  time-scale.  Generally,  we  do  not  have  to  worry  about  e : 


we 


include  it,  in  effect  in  T . [ From  the  point  of  view  of  the  ensemble,  i Sjjl , 

29  J ' ' 

Aq  ,cOq,  also  Tq  and  doppler  l + {3  ( = g),  can  be  considered  as  independent  ran- 
dom parameters,  with  some  appropriate  p.d.’s,  ultimately  to  be  established 
empirically.]  Other  plausible  waveforms  are: 


S^.Sj^Oje 


COS  <J0 


oj(t_£j)  ; 


— 0,  elsewhere  . 


(3.36) 
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For  some  other  types  of  interference  S (t,  R).  might  belong  to  a (gaussian)  noise 

1 J 

process,  or  S^.  might  effectively  persist  for  a long  time  (-00  <t<  co),  and  so  on. 
One  task  of  the  general  study  is  to  determine  the  Sp.  and  study  their  statistical 
properties. 

(3)  Narrow - band  Sources;  An  Example: 

For  practical  consideration  most  of  the  interfering  sources  are  narrow- 

band:  in  this  definition  the  bandwidth  Af^  of  the  source  emission  is  0^10  1 f^ , 

or  smaller,  where  fQ  is  the  carrier  or  central  frequency  of  the  band.  Even  in 

the  so-called  "wide-band"  situations  of  measurement,  Af  is  only  the  order  1 

4 b 

or  2 megacycles,  at  carriers  0(50  me),  in  part  because  of  the  great  dynamic 
range  required  of  the  receiver  to  measure  signals  of  such  bandwidths  without 
serious  distortion. 

The  narrow-band  (n.b.)  condition,  in  conjunction  with  far-field  operation, 
enables  us  to  simplify  our  preceeding,  general  relations  considerably.  In 

place  of  the  real  signal  sources  S^.  used  heretofore,  it  is  convenient  to  employ 
the  complex  n.b.  signal 


Id)  T 

S^tjSjR)  = S0(T,t;|)  e 


icpo(T’t;|) 


(3.37) 


where  S is  a complex  envelope;  A , cp  are  real  envelopes  and  phases.  Then 
u 0 0 

(3.29),  as  it  appears  in  (3.28),  becomes 


fa)  “ 1”  So  [(1  + fo)  > t;l]  e'2" ' [f  " (1  + W f,]  1 


dt 


(3.38) 


where  is  $ { S.l  . 

0 ^t(  Of 

Next,  if  the  apertures  are  essentialfy  frequency-insensitive  to  changes  in 
frequency  in  the  region  f +Af  , where  the  n.b.  source  has  any  significant  out- 
put,  we  can  rewrite  (3.28)  with  the  help  of  (3.38)  for  distributed  sources  Sj^  , as 
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-i<D  T 

0 Oj 

u.(t) 

1 4ttR.  R 

J 


hv°’fo)Le  df/v  ~Aife,f°) 


tooWI/c 


f “ 'iujTo 

J df'e  Vj  [a  + P)  f'5  f - (1  + P)  f0  - (1  + P)  f 0 ’ 

J — Oft 


(3.39) 


where  now  (^j)  is  given  by  (3.  38),  with  the  appropriate  arguments.  The  double 
integration  over  f and  f'  in  that  order,  gives  us  finally  (with  1 + 3 = g here): 
i|W)nft-T0/p,j 


U(t).  = 


0 


4trgR 


^r(-1r  fo/c’  fo)  T"[4  fo/°>  fo ! 


(3.39a) 

where  Cl"  S^( t-T^y/g,  t;  f 1 RU  is  another  generalized  beam  pattern, 

cf.  (3.9). 

Now,  if  the  drive  S ( ; § |R).  at  each  element  d§  of  the  source’s  aperture 

is  the  same,  i.e.,  Sj.  is  independent  of  §,  Eq.  (3.39)  takes  the  still  simpler 
J * — 

form  (with  1 + 3 = g) 


-iurT,, . 

0 0 OJ  * 

U.(t)  = — — — Ct  /-i 
r ; 4ttR.  R 
3 


^fo/c>fo)*i(Vo/c’fo)/_t 


itrt 

e df 


Jot  -iu/T 

e vr(^f',f-^f0-  ^f'|R)0  df'  > 

-oo 


(3.40) 


wherein  the  beam-patterns  are  now  factored  out.  Following  the  same  steps  used 

to  obtain  (3.39a),  we  integrate  next  over  f,  to  get 

-iu).T 
0 Oj 


Jr.  fo/°>  fo  j ^(iR.  fo/c’  f( 


l 


ot  — iu/T, 


0 


X I e ~ YI(pf/;t)()  e 


ig(uu  +u/ ] t 
^ 1 ' df' 


(3.41) 


where  Y ( ;t)Q  is  the  time-varying  system  function  of  the  source  (envelope),  e.g., 


130 


27 


yi<M»0  = f"s0(V)e-dT^1S(Vi)!  . in  (3. 40). 

J “00  ( 0] 

Carrying  out  the  integration  over  f ' gives  us  finally  the  desired  result 


(3.41a) 


U.(t)  = 


i(‘-To>) 


4ugR 


*0  (t"T0  A’ ')  ^(-R  f0  /°»  oWVoA’MI.  • (3'42> 


With  the  small  doppler  in  force  in  all  our  applications,  we  can  set  ^ = 1 every- 


where, except  in  guo^,  so  that  our  narrow-band,  complex  waveform  from  the  j 
source  becomes  finally  the  expected  result 


th 


30 


Uj(t)  " 


lgm 


,(l-To) 


4ttR 


So(t-To’t)fi’R(‘Afo/°’fo)tfi(iRfo/c’fo)|.  • ,3-43) 


If  a point-source  is  assumed,  then  Cl^  -*  <5^0,  f^  = A^O,  f^  , a real  constant, 
representing  an  omni-directional  radiator.  Also,  if  the  source  is  not  time- 
varying,  we  replace  S0^t-TQ,tjby  S0(t-T0j  in  (3.43). 


Example: 

As  a simple  example,  let  us  consider  the  ignition  transient  (3.35),  with  a 
point- source  emitter,  and  a duration  long  enough  so  that  S^.  is  narrow-band 
vis-a-vis  the  receiver's  beam-pattern  £7  . In  terms  of  the  complex  signal  (3.37) 
we  have 


/ -iv  'V' 

S = (A0  e 


A 

— jo)  e 

a o i 

S = A e J 
0j  Oj 


[tsT  ] . 


(3.44) 


Equation  (3.43)  applies  here,  and  with  (3.44),  is 

-Ai(°AH(--W°’fo) 


V(t)  * 


lAoe 


4rrR 


’ [‘o<t<to+Tsi  (3'45) 


= 0,  elsewhere 

including  spreading  loss,  doppler,  and  receiving  beam-pattern.  Other  examples 
are  readily  constructed  in  a similar  way.  Observe,  of  course,  that  the  physical 
signal  is  always  /?c  lE(t)  here. 
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D.  Summary  Remarks: 

In  Section  3,  we  have  shown  how  to  determine  the  (scalar)  interference 
field  N (t,  R)  and  more  particularly,  the  typical  waveform  U.  into  the  receiver 
following  the  receiving  aperture.  This  has  been  done  at  both  the  general  level 
and  the  specialized  level.  The  former  is  needed  for  two  reasons: 

(1)  we  will  require  a general  model,  in  order  to  handle  all  the  possible 
tj-pes  of  interference  at  some  future  stage  of  the  study; 

(2)  we  need  the  general  approach  in  order  to  establish  the  specific  approx- 
imations and  assumptions  required  for  the  special  cases. 

In  all  cases  we  include  doppler,  spreading  loss,  waveform  structure,  and  per- 
tinent geometric  features,  i.e.,  beam-patterns.  At  the  specialized  level,  we 
have  focussed  on  a typical  example,  characteristic  of  the  kind  of  interference 
we  may  expect. 


4.  STATISTICAL  MODELS: 


Here  we  shall  derive  the  pertinent  statistics  for  the  Basic  Statistical  Model 
(BSM)  described  in  Sec.  2.  We  shall  also  develop  expressions  for  the  first-order 
p.d. , covariance  function,  and  a variety  of  other  moments  of  the  equivalent  model 
(ESM),  based  on  the  relations  (2.5)-  (2.10)  presented  earlier.  Our  principal  in- 
terest here  is  in  the  received  process,  X(t),  whose  physical  structure  is  specified 
by  Eqs.  (3.17),  (3. 18)  genei'ally, and  by  Eqs.  (3.25),  (3.28),  (3.33),  (3.34), 
(3.39a),  (3.42),  (3.43),  (3.45),  in  particular. 

A.  The  (BSM)  Poisson  Process,  X(t): 

Let  us  express  X(t)  in  somewhat  more  statistical  detail  than  is  indicated 
in  Sec.  3:  we  write 


X(t)  = 


j 


(4.1) 


where  the  U.  are  the  typical  waveforms  derived  in  Sec.  3;  Hs  a set  of  coordi- 
nate variables,  and  9 are  a set  of  (time-independent)  statistical  parameters, 
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e.g.,  scale  (amplitude),  epoch,  doppler  frequency,  etc.,  cf.  the  Example, 

(3.45),  for  instance.  Then,  U is  the  basic  waveform  of  each  random  "event," 
as  indexed  by  j — the  "event"  here  is  the  generation  of  an  ignition  transient  — 
and  X.  = f t.,  0.,  cp. ) are  the  coordinates  of  the  point  in  space  where  the  event  takes 
place.  This  point  is  randomly  located,  in  a small  spatial  element  AS..,  such  that 
P ^1 1 AS  j = °g^wj^  AS.  is  the  probability  that  exactly  one  such  event  occurs  in 
AS^.  With  the  requirement  that  P^N|  AS^  « P^l  JAS..^  , N>  2,  for  sufficiently 
small  AS^,  and  the  condition  that  in  non-overlapping  regions  the  number  of  events 
associated  with  the  physical  density,  a , of  points  where  these  events  can  occur, 

O 

is  independent  for  each  region,  a Poisson  process  of  the  number  of  events  is  de- 
fined, e.g., 

-AS.  ai 


-Ad.  a / 

(N  | AS.  j = (os  AS.^N  e j S/N ! 


(N  > 0) 


(4.2) 


This,  in  turn,  defines  a Poisson  process  for  X(t).  The  key  point  here  is  that  the 
"events"  are  all  independent.  The  details  of  the  derivation  are  given  in  Secs. 
(2.1),  (3.1),  Ref.  8;  we  present  the  results  here  only: 

(1)  It  is  found  that  the  n*~  order  characteristic  function  (c.f. ) for  X(t),  con- 


sidered at  the  time  t 


t , is 
n’ 


£ P(i><(exp  i £^(VM)|-1 

l X 1 I 


th 


n>  1 


(4.3) 


The  corresponding  n—  order  p.d.,  W , is  the  n-fold  Fourier  transform  of  (4.3), 
Here  the  process  "density"  p is  a non-negative  function,  and  the  operator 
p - p(A)(  ) dX  is  a counting  functional,  which  adds  up  all  the  contributions  of 

the  individual  sources,  without  regard  to  their  magnitude;  A is  the  region  in 
which  the  sources  lie,  for  example,  as  indicated  in  Fig.  3.2. 

The  first-order  c.f.  is,  directly: 
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with  higher-order  c.f.  's  similarly  obtained  from  (4.3). 

The  chief  difficulty  with  (4.3)  and  its  transform  is  their  analytic  com- 
plexity. Although  this  is  an  "exact"  model  of  the  urban  radio  noise  environment 
under  the  basic  assumptions  presented  earlier  (in  Sec.  2),  it  is  not  practically 
amenable  to  closed-form  results  in  the  distributions,  which  severely  limits  its 
analytic  usefulness  whenever  these  are  directly  required  (i.e.,  in  system  de- 
sign, etc.).  However,  in  the  case  of  the  moments,  and  particularly  those  of 
the  first-  and  second-order,  closed  forms  are  readily  obtained,  and  in  this  re- 
spect our  model  is  practically  useful,  particularly  in  conjunction  with  the 
equivalent  model  (ESM)  (Sec.  2).  Finally,  note  that  in  general,  because  of 
the  distributed  nature  and  secular  variation  of  the  sources  (with  time) 

X(t)  is  not  strictly  a stationary  process.  (Sec.  4-(8)  below,  however.) 

(2)  As  we  pointed  out  in  Sec.  2,  we  can  express  X(t)  as  the  sum  of  two  indepen- 
dent Poisson  processes,  X^+X^,  with  densities  and  . The  former 
represents  a high-density  process,  which  is  asymptotically  gaussian.  Re- 
writing (4.3)  since  p = p + p here,  we  have 

B JD 


(4.  5a) 


(4.  5b) 


Furthermore,  expanding  in  an  Edgeworth  series  gives  us 

G 
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(F)^exp[1&^2+ir<V]>  <4-6> 

where  § = , §n  ] is  a row  vector,  (X^)  = (X^),  . . . , (X^>  a (row)  vector 

of  mean  values,  and 


*b  - [<*(*0  %)>  - ■ <xfe>]B  ’ l^-1-  •••»“]  > (4-6a» 

is  an  (nxn)  covariance  matrix  of  the  process,  X . The  process  X remains 

JD  P 

Poisson,  with  a low-density. 


(3)  The  moments  of  X(t)  =X  (t)  +X^(t)  [ = G(t)  + P(t),  cf.  (2.3b)]  are  readily  found 
from  (4.3),  or  (4.5),  in  the  usual  way  by  differentiating  the  characteristic  func- 
tion. We  have  specifically, 


<X(t)>  = 


5-.  = ° 


(4.7a) 


and 
K 


{x(‘l>t2) 


J Pb^bCi-’5))  ~+/  PD(il{UD(‘l;^i))  (‘‘■7b) 

B ""  D ~~ 


(BSM) 


*1*2 12 ' vh  vti  F\ 


° 


j pw(u(54.£)'-,(‘2;M))e®=Kx(ti>t2)B+Kx(t i>‘2)d  > <4-8> 


cf.  (2.9),  for  the  covariance  function.  For  the  mean  intensity,  we  have 


dX 


(4.9) 
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and  we  proceed  similarly  for  the  higher  moments  of  all  orders  and  degrees. 
[We  shall  need  to  examine  the  higher-order,  higher-degree  moments  in  our  la- 
ter studies,  when  we  seek  to  develop  the  equivalent  model  (ESM)  more  fully.] 


(4)  The  behavior  of  X(t),  in  the  BSM,  for  large  amplitudes  is  not  at  all  easy  to 
determine  from  (4.3)  or  its  Fourier  transform.  However,  for  small  ampli- 
tudes [ i.  e. , values  of  X small  compared  to  the  occasional  large  values  of  X 
when  the  impulsive  part,  X , dominates],  which  occur  most  of  the  time  and 
which  are  attributable  to  the  general  background  portion,  X^,  we  can  show 
the  expected  gaussian  character. 

Here  we  shall  consider  only  the  first-order  density  (n  = l)  with  zero  mean, 
as  an  illustration:  from  (4.4)  and  (4.5)  we  can  write 


2 2 / 2 2 / E[(r2m)  +(r2J 


(h) 


G+P 


(4.10) 


where 


Wb.D  =1(B  D) 


(4.11) 


(B,D) 


P(B,D)(^ 


<u®’dX 


d\  . 


Now  the  summation  in  (4. 11)  can  be  replaced  by  the  finite  series,  with  remainder, 
of  the  expansion  of  exp^-§^  , since  (4.11)  can  be  alternatively  written 


(M 


2 t 2 2 \ / 

"5l(aB  + GD  )/2 


exp 


G+P 


^ ... 

4!  -4id?4  FG 


05 


1_ 


; 0 < 0 


= e(h) 


< 1 . 

(4.12) 


Now  lim 


oo  corresponds  to  I X I -*  0,  and  here  we  have 


lim 

lql-“ 


4 


lim  exp 

?ll  -*00 


<r4 


4! 


->1  , 
(4.13) 
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(with  an  Hermitian  polynomial),  so  that 


2 2 / 
"?iaG+P/ 


lim 


\h 


(F0, 


= lim 


G+P 


2 2 2 
, o^  ^ = o „ + a 
’ G+P  G P 


(4.14) 


which  is  the  expected  gaussian  behavior  at  small  values  of  X^=  X^  + 


(5)  The  process  density  p now  needs  to  be  established  from  the  controlling  geom- 
etry. Here  we  shall  assume  first  that  all  the  interfering  sources  lie  on  a 
single  surface  (not  necessarily  flat),  which  reflects  the  fact  that  vehicular 
traffic  moves  on  the  streets  of  the  urban  region  in  question.  The  receiver, 
however,  may  be  on  or  above  this  surface.  The  process  density  p is  ac- 
cordingly specified  by 


ps(^  dX  = as0  ds  , (4.15) 

and  since  dX  = dtdcp,  (ct  = R),  cf.  Fig.  4.1,  and  X = [ t,0( t,  cp),  cp]  are  the  coordi- 
nates of  a point  in  the  surface  A that  is  available  to  the  receiver  at  O . Thus, 
the  shaded  region  in  Fig.  4. 1 is  that  portion  of  the  overall  urban  region  within 

the  beam  Cl  of  the  receiver.  Accordingly,  from  (4. 15)  we  can  write 
R 

= CTs[t>0R(t’  ‘fit)  ’ %]  6s/' it%  • <4' 16) 

Following  the  analysis  of  Sec.  3.1,  Ref.  8,  with  L = 0 therein,  we  find  that 

pS(i)=as(t’eR’%)tc2  ’ °tahR’  (4-17) 

with 

00SeR  = l,R/R  = hR/Ct  ' <4'18) 

[if  the  receiver  is  in  the  same  surface  as  the  sources  (and  that  surface  is  flat, 
cf.  Fig.  4.1),  then  tt/2,  (hp=  0).]  Here  cr^  is  the  physical  density  of  the 
radio  noise  sources,  i.  e. , the  number  per  unit  area,  a quantity  which  needs  to 
be  determined  for  different  urban  regions  and  different  times  of  day,  week,  etc., 
as  part  of  the  experimental  portion  of  the  study. 
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Figure  4.1.  Geometry  of  an  interfering  source  (at  Pj)  and  a receiver  at 
height  (hj^  > 0),  for  an  urban  region  A containing  independent 
radio  noise  sources.  ( Exaggerated  in  the  vertical  direction.) 
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Fields 

(6)  The  Poisson  Noise/ N (t,  R)  =G(t,  R)  +P(t,  R),  cf.  (2.3a)  are  determined  in 

th  ^ 

similar  fashion.  The  n_  order  c. f.  ofN(t,  R),  cf.  (3.16),  analogous  to 
(4.  3),  becomes  here 


F 

n 


; i C 

Tr  n’~n 


f / 

n 

\ 

exp 

p(A)  exp/i 
\ “ \ 

E c/(vs.;xi) 

-l\  dX 

,4=i 

(4.19) 


where  now  Splays  the  role  of  U eai’lier,  with  R?  therein  replaced  by  r.  = R-  R.  , 

^3  J 

for  each  £ = 1, ..  ,,n  , cf.  (3.8),  and  A is  the  domain  of  all  the  sources,  not 

O O 

just  those  "viewed"  by  the  receiver,  as  is  the  case  in  (1)  - (5). 

We  may  proceed  as  in  (2)  - (5)  above.  For  instance,  the  lower-order  mo- 
ments of  this  ambient  noise  field  are 


p(X)(s£(  t,  R;  X,0))  dX  , 
/***  'vv'  'v^' 


Sj-ij-li-S,).  <4-20> 


for  the  mean,  and  for  the  covariance  function  of  the  field,  we  have 


KNI(ti’5i;t2>E2)  =/A  pW^(ti>5ri,£-)^(‘2.52o.g))e 


VS2;M)>  ^ . <4-21) 


with  R'j-»r^  = R^-R  , etc.  in  The  mean  intensity  at  P(R)  is  likewise 


N (t,  R) 

J. 


M, 


p(A)<  s£(t,  R;  X,0)  ) dX 
0 


(4.22) 


and  so  on  for  the  other  moments.  We  can  parallel  the  anafysis  of  (2)  - (4)  for 
X(t),  here  for  N , in  similar  fashion;  the  details  are  left  to  a later  study.  For 
the  surface  distribution  of  ambient  sources  sketched  in  Fig.  4.1,  the  process 
density  p is  given  by  (4. 17),  (4.18). 


(7)  Some  Special  Results;  The  Narrow-band  Process; 

Let  us  assume  narrow-band  sources,  small  doppler,  a non-time-varying 
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signal,  and  far- field  reception,  so  that  (3.43)  applies  for  U . Furthermore,  we 
have 


(4.23) 


since  the  epoch  e.  [of.  (3.45)  and  comments  following  (3.35)]  is  assumed  ran- 
domly and  uniformly  distributed  (0,2tt).  In  fact,  all  odd,  first-order  moments 

of  U.  vanish,  e.g.,  /u2m  =0,  not  unexpectedly  here.  However,  the  covari- 

j \ <t)  / 

ance  function  K ft  , t^  j is  non-vanishing.  For  narrow-band  processes  here  we 
have,32  generally  [for  real  X(t) ] : 


= Hx(ti*'a) oos  V + btfVb) si 


R&- 


(ViAx)° 


'V 


sin  (ju^t 


T $t2-‘l 


B A 

* 1UJnT ) 
. \*  0 
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( x\ 
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where  B 


X 


is  defined  by 


(v  l2) 1 


Bx(tl,t2)‘*(RX+IAx) 


i J 


with 


RX  2 BX  5 


Kr~~~  ^ B 


2 X * 


(4.24a) 

(4.24b) 

(4.24c) 


(4.25) 


Specifically,  for  the  present  Poisson  cases  we  have,  now  for  the  complex  U of 
Sec.  3C-3  above, 


B, 


(ti,t2)  = e 0 J p<y><(u(y^)u(vM)*)  & . 


(4.26) 


Now  using  (3.43),  (4.17),  and  X = ^X,  } djt  = dX  dcp^,  with  X = R/c  hen 

get  for  (4.26) 
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/ \ 1 f i _ ,2  , _ ,2  °s{X’?R) 

bx(Wt  =-7) J^r!  Kl  — x— 

v ' (4n) 


° ®0 (*l‘  x ~)  ®0 (v T " x ~J 


0=^Vn 


with  i in  (7  , £7_,  cf.  (3. 43)  given  by 
-*R  R X 


dX  dcpR  ’ 
(4.27) 


A fA  A A "I 

R = - i cos  m sin  0_,  + i sin  cp_  sin  0_  + i cos  ©„ 
~R  [y-x  YR  R ^-y  YR  R ~z  RJ 


e = cos 

ri 


-i 


(hR/cx) 


(4.28) 


Here  the  averages  indicated  by  ( ) are  over  doppler  (e),  cf.  (3.24),  over  "car- 

0 

/W* 

rier"  frequency  f which  also  appears  in  the  beam  patterns,  cf.  (3.43)  again, 
over  epoch  (e)  and  amplitude  , contained  implicitly  in  the  complex  envelope 
[cf.  (3.45)  in  a particular  instance].  As  is  evident  from  (4.24)-  (4.27),  X(tj 
is  not  in  general  a wide-sense  stationary  process,  since  ^ Kx(lt2_tl=T  l) 

here,  unless  certain  conditions,  readily  met  in  practice  are  obeyed  [vide  4- (8)  below] . 

With  a uniform  beam  in  the  receiver  and  point  sources  j^-*  Aj^O,  f^jl,  both 
insensitive  to  frequency  over  a band  ±Af  about  f , the  average  carrier  frequency,  Eq. 
(4.27)  may  be  simplified  to 

R 1 , . /„  - \ , 2 f CTS(X^r) 


B 


X 


. M°’fo)'2/CTs’ 

(4tt)  V 7 JA 


x / e 


’ 10J  £T 
0 A 


/ *• 


dA  dcpR  , (4.29) 


with  proper  attention  to  the  domain  A implied  by  the  nonvanishing  of  | G 


R1 


If  U ^ U , cf.  (2.9),  and  certainly  p ^ p , then  B becomes 

AD  JD  13  U X 


Bx(tl,t2)B  4 


(4.30) 
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where  each  B in  (4.  30)  is  a term  of  the  form  (4.27);  a = a + a , as  well, 
X.  13  JD 

with  cr  » o In  this  way,  by  an  anatomization  of  the  waveform  we  arrive  at 
a corresponding  anatomization  of  the  statistics.  From  an  experimental  point 
of  view  the  covariance  function  is  also  an  important  descriptor  of  the  noise 
process. 

Finally,  we  mention  that  there  is  a spectral  statistic  associated  with 
here.  This  is  the  cointensity  spectrum  [Ref.  8,  Sec.  8.5],  defined  by  the 
double  Fourier  transform 


“Vr1^ 


dt„  dt. 


(4.31) 

This  quantity  measures  the  amount  of  statistical  correlation  between  components 
of  X(t)  at  frequency  f and  frequency  f . If  X(t)  is  stationary  (as  it  can  be  ap- 
proximated  in  many  important  cases  [vide  4- (8)  following], then  (4.31)  reduces  to 


yx(fl’f2)=yx(fl)6(f2'fl)  > <4'3: 

where  57  (f)  is  the  familiar  (auto-)  spectral  (intensity)  of  the  (now  stationary) 

A 

process  X(t),  The  auto-  and  cointensity  spectra  are  also  quantities  of  experi- 
mental interest  in  building  workable  models  of  this  class  of  radio  noise. 


(8)  Macrostationarity:  The  interfering  radio  noise  process  X(t)  is  not  strictly 
or  even  wide-sense  stationary,  because  of  (a), the  finite  duration  of  the 
emitted  signals,  and  (b)tyhe  secular  variation  of  source  density  (p)  and  signal 
level  with  geography  (A)  and  time  of  day,  week,  etc.  However,  in  most  in- 
stances X(t)  can  be  regarded  as  effectively  stationary  in  the  wide  sense  at 
least,  from  the  viewpoint  of  the  usual  observer,  whose  receiving  system  is 
influenced  by  this  noise.  This  is  because  (lb the  number  of  such  sources  in 
A is  large  enough  that  there  is  always  at  least  one  or  more  emissions  going  on 
at  any  given  time,  and  (2), the  average  emission,  S^,  has  a duration  long  com- 
pared to  the  time  it  takes  to  propagate  over  the  maximum  effective  distance  A 

max 

between  source  and  receiver  in  the  domain  of  possible  sources,  A . Beyond 
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this  distance  we  say  that  the  spreading  loss  ( ~X-2  in  energy)  is  so  great 

that  the  source  projects  negligible  interfering  energy  into  the  receiver. 

Thus,  if  the  maximum  effective  dimension  is  0(30  km,  i.  e. , X -30km/ 

(c  = 3*10  km/sec.  ) = 0. 1 msec. , or  less,  a dimension  very  probably  on  the 

large  side  in  practice,  a 1 msec,  emission  will  persist  for  ten  times  the 

time  it  takes  the  wavefront  to  cover  X . The  result  is  what  we  call  here 

max 

a macro  stationary  process,  i.  e. , one  that  is  not  truely  stationary  in  the 
analytical  sense  but  which  can  be  treated  as  such  for  time  periods  during 
which  the  secular  changes  in  the  source  conditions  are  small.  Thus,  in  a 
period  of  heavy  traffic,  say  4- 6p.m.  of  a city  weekday,  we  may  expect 
that  the  secular  time  scale  is  0(VS  hr.  - 1 hr. ),  much  longer  than  a typical 
signal  period,  and  certainly  much  longer  than  any  one  (non-zero)  signal 
duration  and  X 

max 

Mathematically,  the  conditions  for  macrostationarity  may  be  established 
as  follows:  we  start  with  (4.27),  for  the  usual  narrow-band  signals,  and  expand 

A A , A A 

the  signal  envelope  S0,  viz.  S ( t— X;  6)  = S (t;  9)  - XS  (t;0),  so  that  a typical  con- 
V U ^ U ^ U /v*' 

A 

dition  for  neglecting  during  the  time  X is 

0 max 


r a 2\  1 

Vp 

r A 2\1 

S' 

4-> 

o 
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2 » X 

max 

(vt:e> ) 

ya 


(4.33) 


this  defines  our  characteristic,  interfering  signal  of  "long  duration"  here.  Ac- 
cordingly, we  can  now  replace  sVt^  2~^;0)  So(tl  2’®)  in  which  in 

(4.24c)  in  turn,  allows  us  to  write  for  the  covariance  of  the  process 


10)  T 


> ‘2=tl  + T> 


(4.34) 


^(‘i’V7)  e 0 (so(‘i/)X(tt;£))( 

/i 

with  = ^exp^io>QT£^)c  the  characteristic  function  of  the  doppler,  e,  where 

it  is  assumed  that  f is  known, ^ and  where 

CA  = ~ / I^rI2  l^fl2  as(X,CpR)  dX  dcpR/X  ’ (4.34a) 

a (real)  constant  (>  0),  depending  only  on  the  geometry  of  the  receiver  vis-a-vis 
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the  interfering  sources.  Next,  if  the  epochs  of  initiation  of  the  are  uniformly 
distributed  in  time,  over  a period  comparable  to  or  larger  than  the  mean  dura- 
tion of  a typical  interfering  source,  a situation  entirely  reasonable  in  practice, 
we  see  that 


the  (complex)  covariance  of  the  source  envelope;  then  0 = (e,  0'),  0' ^amplitude, 
durations,  etc.  Thus,  we  have  for  (4. 34) 


Kx(tl’tl+T)  “ CA  ^e|Fl(‘‘V)i 


V A 

iu)QTj  e Kg  (t) 
e 0 


(4.36) 


which  is  a function  of  t only,  and  hence,  X,  in  the  above  sense,  is  "macrosta- 
tionary. " With  ignorable  doppler,  (4.  36)  reduces  to  the  still  simpler  form 


(no  doppler): 


Kq  (t) 
S0 


(4.37) 


For  a typical  interfering  signal  like  (3.35),  we  get  at  once 


Kx(WT)iA0CA<(cos'"0(1'E)T) 


aocacosV  > 


(4.38) 


this  last  for  ignorable  doppler.  The  covariances  for  other,  more  complex  inter- 
fering signals  are  readity  found  in  similar  fashion. 

Finally,  this  property  of  macrostationarity  is  very  important  for  experi- 
mental measurement,  for  it  allows  us  to  construct  an  experimental  data  ensemble 
directly  from  a continuous  record  of  the  output  of  a single  receiver  subject  to 
this  mixed-process  (Gauss  and  Poisson)  interference.  The  procedure  is  to  divide 
the  record  up  into  equal  intervals  long  compared  to  the  reciprocal  of  the  band- 
width of  the  interference  (as  observed  in  the  receiver),  but  short  compared  to 
the  (much  longer)  period  during  which  secular  variations  in  the  traffic  intensity, 
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and  hence  ignition  noise  level  and  character,  begin  to  manifest  themselves.  Be- 
fore we  calculate  from  these  ensemble  data  the  various  desired  statistics  (mo- 
ments, p.d.'s,  etc.),  it  is  necessary  to  establish  the  validity  of  the  data,  e.g., 
we  must  verify  that  we  have  a valid  ensemble  — one  where  the  member  func- 
tions are  independent  and  constitute  a homogeneous  set.  For  these  purposes  vari- 
ous statistical  tests  are  available;  illustrative  applications  have  been  discussed, 

34 

for  example,  by  Middleton  in  connection  with  ocean  reverberation,  and  by 
35 

Plemons,  in  an  experimental  study  of  the  scattering  of  acoustic  waves  from 
lake  surfaces.  Applications  of  these  concepts  to  our  present  problem  is  im- 
mediate and  direct. 


B . The  Equivalent  Model  for  X(t): 

As  we  noted  in  Sec.  2,  the  equivalent  statistical  model  (ESM)  is  postulated 
to  have  the  form 


X(t)  = a(t)  Z(t)  , (4.39) 

cf.  (2.5),  where  Z(t)  is  gaussian({Z)  =0  ) , and  a(t)  is  slowly-varying  vis-a-vis 
Z(t),  which  is  narrow-band  about  some  carrier  frequency  f^,  determined  by  the 
receiver.  The  process  a(t)  is  also  zero-mean,  i.  e. , (a)  =0.  Both  a and  Z are 
statistically  independent. 

First,  let  us  show  that  for  small  amplitudes  of  X,  X is  gaussian,  at  least 
in  the  first-order,  as  it  must  be  to  be  consistent  with  the  Basic  Model  [cf.  A, 
(4)  above].  Since  Z(t)  is  normal,  we  can  write  the  first-order  p.d.  of  X as 


Letb  = l/a,  [-a><a<co,  and.’. 


-»  < b < co] ; then  (4.  40)  becomes 


(4.  40) 


V.Z 

\ f3/y,TT(T  n / \ ^/y.riry  / 

Since  X is  small,  we  expand  (4.  41)  and  get  (with  K a normalizing  constant' 


(4.41) 
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lW.Z  = K ([<lb r‘>  - lfclx2/24  + • ■ ])b  = 

/<lbl><l' 


2/2A: 


,-lv-l 


(4.42  a) 
(4.42b) 


which  shows  the  expected  normality  for  small  amplitudes. 

For  large  amplitudes  ( | X j — • «)  it  is  known  empirically  that  W^(X)  falls  off 
much  more  slowly  than  if  X(t)  were  gaussian  only.  This  is  because  of  the  more 
or  less  discrete^large  impulsive  transients,  that  greatly  exceed  the  background, 
and  which  tend  to  take  on  the  amplitude  distributions  of  the  impulses  themselves. 
To  generate  a class  of  first-order  p.  d.  's  in  a that  yields  the  proper  behavior  at 
large  amplitudes  in  X,  we  follow  the  procedure  suggested  by  Hall  [Ref.  5,  pp. 
13-15]  for  his  model  of  atmospheric  noise,  except  that  here  the  impulsive 
events  are  all  independent,  unlike  the  situation  occurring  in  the  lightning  dis- 
charges which  generate  atmospheric  noise.  Accordingly,  we  consider  the  re- 

-1  2 
ciprocal  process  b(t)  = a(t)  , and  require  it  to  have  a first-order  x -p.d.  of  v 

degrees  of  freedom,  with  scale  parameter  ct^,  viz: 


w1(b) 


(v/2) 


v/2 


v-1 


T(v/2)  cr 


2/n  2 
-Vb/2o 

b 


( -co  < b < <»)  ; 

cr,  > 0 ; v>  0 . 
b 


(4.43) 


Then  one  readily  shows  that  for  a = l/b, 
v/2 


Wx(a)  = 


iv/ 21 


r(v/2) 


/o  2 2 

-v/2ct,  a 

■ — v — 1 b 

a e 


(— oo<a<°=)  , 


and 


lim  W (a)  = lim  | a | V 1 , ( v > 0)  , 


1 . 

i' 7 

II  v1  1 

> 

1 a 1 -*  co 

8 

T 

ctS 

(4.  44) 


(4.44a) 


which  is  the  desired  form  of  large-amplitude  dependence  for  X,  since  |a|-»  « im- 
plies |X |->  co  and  from  (4.41),  we  have  ( remembering  that  |dZ/dX|  = l/a) 
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W,(X).  , = lim  ----- — ~~=  lim  

'al'"CD  |X|-«  |X|V+  |X  | — co  |x| 


v+1  * 


(4.44b) 


Our  next  step  is  to  obtain  W^(X)  explicitly.  For  this  we  use  (4.44)  in  (4.41), 
2f  / 2 2/ 2 


Wl(X)  = j 

J — r 


,-v-2  Av 
a I — 7=7f  da  ; 


= jy/2)v/.2  -v 
AV  T(v/2)  b ’ 

and  get  finally 


wi«“  r(v/2) 


r[(v-n)/2]  (y7 

y/tT 


/ 2 2 / 2\ 

(x  + v°2/ ab) 


(•J  +l)/2 


Letting 


(4.45) 


, v > 0 . (4.46) 


°0  HV^  az/CTb  mid  V + 1 = y (>  1)  , (°Z>  °) 


(4.47) 


we  can  express  W^(X)  in  the  somev/hat  more  compact  form 


(4.48) 


Note  that  (4.44b)  is  indeed  the  behavior  as  |x|-»<»,  and,  in  particular,  that 


= r[(v-n)/2]  v 

v T(y/2)Jr\  0 


(4.  48  a) 


The  characteristic  function  (c.f. ) for  W^(X)  is 
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(4.49a) 
(4.49b) 


v = y-1 
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where  K , is  a modified  Bessel  function  of  the  second  kind.  The  various  mo- 

^ l 37 

ments  <X2^),  when  they  exist,  are  found  most  easily  from 


YW 


X2kdX 


2k 


y/2 


(-l)kr(k  + l/2)(a0) 

[(3  - v)/2)  vAT  5 


y > 2k+l 


< 2k+l 


(4.  50) 


In  a similar  fashion  we  can  readily  obtain  the  first-order  statistics  of  the 
envelope  E = | a(t)  | e(t)  [cf.  Sec.  2,  Eqs.  (2.10)].  Wehave 


W1(E)  - (w^E 


wx(E|a)  = 


Ee 


r,2/o  2 2 

-E/2a  crz 


a cr 


(4.  51) 


which  from  (4.  44)  for  W^(0)  is  found  to  be 


wi(E)  = cye/(e2  + <to 


where 


y/2 


, 0 < E < co  , 


=2r(y/2)a^  ^(  Y/2-1)  . 


(4.52) 
(4.  52a) 


As  alreadj'  noted  in  Sec.  2,  cf.  remarks  following  Eqs.  (2.10),  the  first-order 
p.d.  of  the  phase,  cp,  is  uniform  on  2tt.  The  various  first-order  moments  of 
the  envelope  are  also  easity  obtained.  We  have  directly 
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W, 


Ek+1  dE  a0  F(k/2  + y/2  " k/2  " X) 


o 4,2  2\ 
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, y > k+2 
, y ^ k+2 


(4.  53) 
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Finally,  note  that  for  small  amplitudes 


Wi(E)  - = 2^E/<E2)^e“E  ^ , E>0  , 

<E2)  = 2/C 


(4.53a) 


which  shows  the  expected  Rayleigh  form  for  the  envelope  [analogous  to  the 
gaussianp.d.  for  the  instantaneous  amplitude,  cf.  (4.42)]. 

The  parameters  and  v are  to  be  determined  empirically.  We  expect 
that  they  will  show  secular  variations  with  the  time  of  day,  week,  etc.  They 
may  also  vary  with  the  urban  region  chosen  for  study.  They  will  certainly  de- 
pend on  the  bandwidth  and  frequency  allocation  of  the  receiver  used.  In  any  case, 
they  are  a part  of  the  experimental  data  to  be  investigated  in  the  overall  study. 

The  higher-order  statistics  of  X = aZ  need  to  be  studied  theoretically  in 
conjunction  with  the  basic  model  (BSM)  developed  here  [Secs.  3,  4A] . As  noted 
in  Sec.  2,  for  example,  and  with  the  aid  of  (2.  7),  (2.  9),  and  the  results  of  Sec. 
4A-7  above,  we  have 


^azfV l2)  = Ka(*l> l2)  Kz(‘l’ l2)  = 'S/T h)  <“•  54> 

where  K^.  is  given  explicitly  by  (4.27)  or  (4.29).  Moreover,  because  X(t) 

may  be  expected  to  be  macro  stationary  [Sec.  4A-8],  we  may  replace  (4.  54)  by 

Kx(tl’tl+T)  ~ Ka(T)  Kz(T)  5 (4.54a) 

with  considerable  simplification  of  the  higher-order  statistics  of  X. 

Because  of  the  independence  of  the  basic  interference  " events"  in  the  urban 

5 

model,  in  contradistinction  to  the  atmospheric  situation,  we  may  expect  that 
a(t)  may  be  quite  slowly  varying  over  reasonable  observation  intervals,  t^+  T^ 
so  that  during  such  periods  T,  a(t)  is  effectively  constant.  If  this  is  the  case  in 
conjunction  with  macrostationarity,  calculation  of  the  higher-order  statistics  of 
the  equivalent  statistical  model  (ESM)  is  greatly  simplified.  For  example,  we 
have 
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KX=az(tl,t2)  = <^2^)  KZ(T)  = ICK(T)BSM  5 tEq*  (2*9)and  Sec*  4A-3,8],  (4.55) 
which  determines  K and,  therefore,  aid  orders  of  p.d.  's  for  Z,  since  Z(t)  is 

/j 

postulated  to  be  zero-mean  and  gaussian.  From  (4.50)  we  have 


(Y>3)  . 


(4.56) 


The  relative  invariance  of  a (t ) with  time  is  accordingly  a very  important  question, 
to  be  tested  experimentally. 

Finally,  even  though  it  may  be  established  empirically  that  y is  3 (or  less), 

38 

and  consequently  (a2)  and  (X2)  ->  <»,  our  ESM  may  still  be  used.  Physically, 
of  course,  (X2)  must  be  finite.  What  probably  occurs  in  reality  is  that  not  one, 
but  two  laws  for  a(t)  are  required  at  large  amplitudes:  one  with  a value  y(  < 3), 
say,  for  large  amplitudes,  up  to  some  very  high  level,  at  and  above  which  there 
is  not  onty  possible  saturation  at  the  front-end  of  the  receiver,  but  also  a modi- 
fication of  y,  now  larger  than  3,  to  insure  the  necessary  finiteness  of  (X2). 


5.  CONCLUDING  REMARKS: 

In  this  first  paper  (Part  I),  we  have  laid  the  foundations  for  a workable 
analytical  model  of  urban  radio-noise  interference.  Although  attention  here  has 
been  directed  primarily  to  the  environment  of  vehicular  mobile  land-transporta- 
tion, the  general  results  and  approaches  used  are  available  for  other  classes  of 
noise  mechanisms,  both  civilian  and  military.  In  so  doing,  we  have  also  raised 
various  questions  that  can  only  be  resolved  by  experiment  and/or  simulation. 

Our  analytical  model  includes  the  physics  of  propagation,  coupling  to  the  medium 
by  means  of  the  source  and  receiving  apertures  (beam  patterns),  source  and 
platform  motion  (doppler),  the  general  geometry  of  sources  and  receiver,  and 
general  waveforms  for  die  interference,  both  broad-  and  narrow-band.  The 
basic  statistical  model  (BSM)  is  Poisson  in  the  location  in  space  of  possible 
sources,  leading  to  a Poisson  process  for  the  received  waveforms  [which  is 
shown  to  be  essentially  stationary  (macrostationary)  for  the  type  of  emissions 
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considered  here] . The  fundamental  assumption  of  independent  transient  sources 

appears  entirely  reasonable.  By  suitable  modification  of  Hall’s  atmospheric  mod- 
5 

el,  we  have  constructed  an  equivalent  statistical  model  (ESM)  that  has  the  great 
advantage  of  analytic  tractibility  in  its  probability  distributions,  unlike  the  BSM, 
for  which  generally  no  closed-form,  or  even  simple  series  developments  seem  pos- 
sible. (This,  of  course,  does  not  preclude  simulation  of  the  BSM  and  its  distri- 
butions, if  need  be.)  On  the  other  hand,  our  BSM  provides  us  with  the  explicit, 
detailed  forms  of  the  various  moments  of  the  received  process,  e.g. , the  mean 
(usually  zero),  the  average  intensity,  (X^),  the  covariance  function, 

t0x(**)>]1  and  whatever  higher-order  and  higher  degree  moments 
desired.  By  means  of  these  moments  we  related  the  non-empirical  por- 
tion of  the  equivalent  model  (ESM)  to  the  detailed  physics  and  geometry  of  the 
underlying  interfering  sources  and  their  radiation.  The  "equivalent"  model  re- 
quires empirical  calibration  [for  example,  two  parameters  (Y,  °q)>  cf-  (4.48) 
are  being  measured  at  TSC] . Other  measurements  [ such  as  testing  the  invariance 
of  the  scale  process,  a(t)]  and  analytical  extensions  [higher-order  p.d.'s]  also 
mentioned  in  preceeding  Secs.  2 -4,  are  currently  under  investigation. 


I;  <x( 

may  be 
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GLOSSARY  OF  PRINCIPAL  SYMBOLS: 


A.  Aq 
a(t) 

signal  amplitude 
scaling  process 

AI(l’f) 

aperture  weighting  (source) 

AR(l’f) 

a\'aRiCl\ 

timeform  of  aperture  weighting 
aperture  weighting  (receiver) 
beam  patterns 

B.  BSM 

Basic  Statistical  Model 

PRS’  8SR 

P 

normalized  doppler  speeds 
a combination  of  doppler  speeds 

C.  c 

velocity  of  propagation 

D.  6(  ) 

a della  function 

E.  ESM 

Equivalent  Statistical  Model 

e(t) 

envelope  of  Z(t) 

1 

a receiver  coordinate 

e 

sum  of  dopplers 

A 

t 

epoch:  observer  vs.  initial  source  emission 

F f f 

’0 

F.,F 
1’  n 

frequencies 

characteristic  functions 

G.  G(t,R) 

gauss  noise  field 

G(t) 

a gauss  process 

GI 

source  function 

Y 

parameter  of  a p.  d. 
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II. 

hRI 

* 

weighting  function  of  RF-IF  receiver 
(including  receiving  aperture) 

hR 

A A A A A A 

A, 

height  of  receiver  above  source  plane 

I. 

i„;i  j1  »i  >*  j 

R x x y y z’ 

i' 

Z 

: unit  vectors 

J. 

K. 

BX  ’ Si  ’ ICN  ’ Ka  ’ 
5 

Av> 

Kz 

: covariance  functions 

: a source  coordinate 

L. 

5c 

"sine"  component  of 

field  of  a noise  source 

A 

domain  of  interfering  sources 

X 

a (time)  variable 

X 

a set  of  coordinates 

L 

max 

max.  dimension  of  source  aperture 

M. 

M- 

: 

1 + a doppler  speed 

N. 

N^R) 

interference  field 

V 

vector  wave  number 

V 

a statistical  parameter 

o. 

: 

angular  frequencies 

P. 

Poisson  noise  field 

P(t) 

a Poisson  process 

<P0i«P 

phases 

^R 

azimuthal  angle 

Q. 
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R.  Rx 

p(M 

r,  r',  R,  R',  R 

"cosine"  component  of  I<^ 
Poisson  source  density 
distance  ( and  vectors) 

S'  °B>aD-aS 

°0 

A 

so 

Sj(t,R) 

Sj(f,5) 

physical  densities  of  sources 
parameter  of  a p.d. 
complex  signal  envelope 
interfering  noise  source 
amplitude  spectrum  of  source 

T.  Tg 

To 

,-ar  ’I'm  >1r0 
0 

signal  duration 

time  delay  on  path  of  propagation 
transformations 

set  of  parameters  of  interfering  signal 

T 

correlation  time,  t -t 

Zj  _L 

"’VVVr 

times 

u.  u. 

J 

typical  received  waveform  from  j ~ 1 source 

V V V 

P R 

volume  of  source  and  receiver  aperture 

V . V 

~R’~S 

velocities 

W.  Vx 

WPW! 

auto-  or  cointensity  spectrum 

p.d.  's 

X.  X(t) 

the  interference  process  in  the  receiver 

Y-  YI 

% 

Z.  Z(t) 

( time- varying)  system  function  of  source 
bi-frequency  function  of  source 
a gaussian  random  process 
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LIST  OF  FIGURE  CAPTIONS:  (DOT:  DM,  Paper  #1) 


Fig.  3.1 

Geometry  of  interfering  source  and  the  receiver 

Fig.  3.2 

Schema  of  a typical  distribution  of  urban  noise  sources  and  a 
receiver  (R) 

Fig.  4.1 

Geometry  of  an  interfering  source  ^at  P^  and  a receiver  at  height 
^h^>_  0^  , for  an  urban  region  A containing  independent  radio  noise 
sources.  ( Exaggerated  in  the  vertical  direction. ) 
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